Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,141,631 C→G M3I (ATG→ATC)  unknown ← Minor pilin of type IV secretion complex, VirB5

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,141,6310CG98.4% 170.7 / ‑3.5 61M3I (ATG→ATC) unknownMinor pilin of type IV secretion complex, VirB5
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (31/29);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

CCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACGTTGTTTATAAATCGGAATCCAGTCATTAGCGGCTTCGCCGACTTCTTCGATAATTTGATGTAGGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCGGTTTCGAGCCCTTTACGCACTAAGAAAGTACGGGAATCCGGTTGTAATTCACGGATAATTTCAAA  >  1/1141386‑1141871
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ccAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGAtttt                                                                                                                                                                                                                                             >  1:828561/1‑251 (MQ=255)
         gCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTGTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGt                                                                                                                                                                                                                                     <  1:592353/250‑1 (MQ=255)
                   cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCtt                                                                                                                                                                                                                              <  1:142079/247‑1 (MQ=255)
                   cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                                                                                                                                                                                                                          <  1:725014/251‑1 (MQ=255)
                   cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                                                                                                                                                                                                                          <  1:691574/251‑1 (MQ=255)
                   cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                                                                                                                                                                                                                          <  1:294045/251‑1 (MQ=255)
                    aCGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATaca                                                                                                                                                                                                                         <  1:337627/251‑1 (MQ=255)
                     cGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATacac                                                                                                                                                                                                                        <  1:58906/251‑1 (MQ=255)
                        aTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACaa                                                                                                                                                                                                                      >  1:342920/1‑250 (MQ=255)
                        aTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAAc                                                                                                                                                                                                                     <  1:799367/251‑1 (MQ=255)
                          tccggtctgtGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACgt                                                                                                                                                                                                                   <  1:364165/247‑1 (MQ=255)
                                ctgtGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATacac                                                                                                                                                                                                                        >  1:336927/1‑240 (MQ=255)
                                     aGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATaaa                                                                                                                                                                                                         <  1:567908/250‑1 (MQ=255)
                                      ggCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTGGCTTGAACTAATGCCGGTATGGACAAGACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACgt                                                                                                                                                                                                                   <  1:679189/239‑1 (MQ=255)
                                                   gCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAAt                                                                                                                                                                                                        >  1:237594/1‑237 (MQ=255)
                                                                               gagatacggtttctgTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCaataat                                                                                                                                                              >  1:176362/15‑251 (MQ=255)
                                                                                    acggtttctgTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGAt                                                                                                                                                         <  1:483430/242‑1 (MQ=255)
                                                                                          tctgTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTaaga                                                                                                                                                   <  1:51409/248‑4 (MQ=255)
                                                                                               ccattGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCaat                                                                                                                                                                                             >  1:329920/5‑202 (MQ=255)
                                                                                                    gttttgcttctTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTAAAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGtt                                                                                                                                         >  1:737734/9‑251 (MQ=255)
                                                                                                             ctTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACt                                                                                                                                <  1:644079/251‑1 (MQ=255)
                                                                                                               tttGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTgg                                                                                                                              <  1:598101/251‑1 (MQ=255)
                                                                                                                       ctgttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGa                                                                                                                      >  1:728130/26‑251 (MQ=255)
                                                                                                                       ctgttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGa                                                                                                                      >  1:135463/26‑251 (MQ=255)
                                                                                                                         gttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGtt                                                                                                                                         >  1:336742/24‑230 (MQ=255)
                                                                                                                         gttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTgg                                                                                                                              <  1:133564/218‑1 (MQ=255)
                                                                                                                         gttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGt                                                                                                                             >  1:84945/24‑242 (MQ=255)
                                                                                                                         gttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAAt                                                                                                                    >  1:20359/24‑251 (MQ=255)
                                                                                                                         gttgaacagtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAAt                                                                                                                    >  1:34751/24‑251 (MQ=255)
                                                                                                                                 gtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACCAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGcc                                                                                                                                                                             >  1:689314/16‑186 (MQ=255)
                                                                                                                                 gtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAg                                                                                                                                                                                       >  1:267871/16‑176 (MQ=255)
                                                                                                                                 gtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTg                                                                                                                                             >  1:430051/16‑218 (MQ=255)
                                                                                                                                 gtttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGt                                                                                                            >  1:580475/16‑251 (MQ=255)
                                                                                                                                   ttgtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACgt                                                                                                                                                                                                                   <  1:594438/133‑1 (MQ=255)
                                                                                                                                     gtaacatactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGtttt                                                                                                         >  1:773688/12‑250 (MQ=255)
                                                                                                                                          atactgaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGc                                                                                                    >  1:571004/7‑250 (MQ=255)
                                                                                                                                               gaTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGt                                                                                                                                          <  1:471331/206‑1 (MQ=255)
                                                                                                                                                aTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAAc                                                                                             <  1:443052/251‑1 (MQ=255)
                                                                                                                                                aTTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAAc                                                                                             <  1:655629/251‑1 (MQ=255)
                                                                                                                                                               gcagcagcaccATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACttcttcaat                                                                                                                                                                 >  1:172880/11‑165 (MQ=255)
                                                                                                                                                                       accATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAg                                                                      <  1:279456/249‑1 (MQ=255)
                                                                                                                                                                        ccATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAAc                                                                                             <  1:208441/226‑1 (MQ=255)
                                                                                                                                                                        ccATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACtt                                                                             <  1:262707/242‑1 (MQ=255)
                                                                                                                                                                         cATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGc                                                                    <  1:372885/251‑1 (MQ=255)
                                                                                                                                                                            caactaCCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAgg                                                                     >  1:280069/1‑247 (MQ=255)
                                                                                                                                                                               ctacCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGttt                                                               >  1:149322/3‑250 (MQ=255)
                                                                                                                                                                                tacCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCg                                                             <  1:842939/250‑1 (MQ=255)
                                                                                                                                                                                  ccGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATcc                                                                                                                                     >  1:874890/1‑177 (MQ=255)
                                                                                                                                                                                  ccGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAg                                                           >  1:56737/1‑251 (MQ=255)
                                                                                                                                                                                        tACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCttt                                                     >  1:897662/1‑251 (MQ=255)
                                                                                                                                                                                             ctagCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGc                                                 <  1:329409/247‑1 (MQ=255)
                                                                                                                                                                                               agCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCa                                                >  1:72966/2‑249 (MQ=255)
                                                                                                                                                                                                             ccGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCACTAAGAAAGTACggg                                 <  1:564945/250‑1 (MQ=255)
                                                                                                                                                                                                                gtaTAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAAc                                                                                             >  1:371586/3‑187 (MQ=255)
                                                                                                                                                                                                                gtaTAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGc                                                                    >  1:267966/3‑212 (MQ=255)
                                                                                                                                                                                                                    agacaagaCAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCACTAAGAAAGTACGGGAATCCGGt                         >  1:187684/8‑251 (MQ=255)
                                                                                                                                                                                                                              gCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAg                                                                                                             <  1:691146/157‑1 (MQ=255)
                                                                                                                                                                                                                              gCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCACTAAGAAAGTa                                     >  1:35442/1‑229 (MQ=255)
                                                                                                                                                                                                                               cAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCACTAAGAAAGtgcgtaaagggctcgaaac                    <  1:750249/245‑19 (MQ=255)
                                                                                                                                                                                                                               cAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCACTAAGAAAGTACGGGAATCCGGTTGTAAtt                  <  1:451251/247‑1 (MQ=255)
                                                                                                                                                                                                                                           gTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTAAGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCTGTTTCGAGCCCTTTACGCACTAAGAAAGTACGGGAATCCGGTTGTAATTCACGGATAATTTCaaa  >  1:184391/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
CCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACGTTGTTTATAAATCGGAATCCAGTCATTAGCGGCTTCGCCGACTTCTTCGATAATTTGATGTAGGATTTCTGCGTTATCCTTACTGGTTGACAGAATTTTTAAGTTTTCACGCCCTAAACCGCTTAAATCTAACTTAGCTAAGGCGGTTTCGAGCCCTTTACGCACTAAGAAAGTACGGGAATCCGGTTGTAATTCACGGATAATTTCAAA  >  1/1141386‑1141871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑