Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,642,358 T→A 6.1% intergenic (+27/+22) dtpB → / ← rsmJ dipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,642,3580TA 6.1% 692.9 / 25.7 196intergenic (+27/+22)dtpB/rsmJdipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (89/95);  new base A (7/5);  total (96/100)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 199 total aligned reads displayed.

AGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAA  >  NC_000913/3642090‑3642630
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
agcaTTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCg                                                                                                                                                                                                                                                                             <  2:241904/274‑1 (MQ=255)
    gttGATTGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGTGCCGCTTGAGACGTTTCCCGTCTATACCAACGGGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAcc                                                                                                                                                                                                                                                                        <  2:272628/274‑1 (MQ=255)
                    gtcGGGGTGTGTGTACTGACGCAGGCTGTCGCGTTCTTGGTGGGTGTCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTGGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCGGGTCACGATGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                        <  4:132989/274‑1 (MQ=255)
                                                                                  tACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTTGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTCCACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTttgcgtggct                                                                                                                                                                                          >  4:293849/1‑274 (MQ=255)
                                                                                                     aTTACCGATCCGCTTGAGACGTTGCCCGTCTATCCCAACGTGTTTGGTAAGCTGGGTCTGGTCACGCTGGGCGTTGAAGTAGTGATGCTTTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgc                                                                                                                                                                        <  2:293522/274‑1 (MQ=255)
                                                                                                            aTCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGACTGCGACGCCGGAAAGCCATTAATTATTATTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACAGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGt                                                                                                                                                                 >  2:200650/1‑274 (MQ=255)
                                                                                                                gCTTGAGACGTGTCCCGTCTATACCAACGCTTTTGGTATGATTGGTCTGTGCACGCTGGGCTTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATAAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                            <  2:55247/275‑1 (MQ=255)
                                                                                                                         gTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGtt                                                                                                                                                                                                                                                                                 <  7:635640/149‑1 (MQ=255)
                                                                                                                              ccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGc                                                                                                                                                                                                                                                                            >  6:865375/1‑149 (MQ=255)
                                                                                                                              ccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGc                                                                                                                                                                                                                                                                            >  8:1452257/1‑149 (MQ=255)
                                                                                                                                 tCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAc                                                                                                                                                                                                                                                                         <  5:152927/149‑1 (MQ=255)
                                                                                                                                   tataCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTAccc                                                                                                                                                                                                                                                                       <  6:396069/149‑1 (MQ=255)
                                                                                                                                     taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTg                                                                                                                                                                                                                                                                     >  6:617755/1‑149 (MQ=255)
                                                                                                                                     taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTg                                                                                                                                                                                                                                                                     <  8:442560/149‑1 (MQ=255)
                                                                                                                                      aCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgt                                                                                                                                                                                                                                                                    >  8:356689/1‑149 (MQ=255)
                                                                                                                                        cAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                                                                                                                                                  >  5:217065/1‑149 (MQ=255)
                                                                                                                                        cAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                                                                                                                                                  >  8:666817/1‑149 (MQ=255)
                                                                                                                                          aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTATCCTgtttta                                                                                                                                                                                                                                                                <  8:546961/149‑1 (MQ=255)
                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttag                                                                                                                                                                                                                                                               >  8:625836/1‑149 (MQ=255)
                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttag                                                                                                                                                                                                                                                               <  5:1423365/149‑1 (MQ=255)
                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACCTTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc                                                                                                                                  >  2:124130/1‑274 (MQ=255)
                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTTGCGCGTAGTCCGGTCGCTTTACCACCCCCCGGTTGGTCGc                                                                                                                                  >  2:33667/1‑274 (MQ=255)
                                                                                                                                               tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCt                                                                                                                                                                                                                                                                      >  4:180322/1‑138 (MQ=255)
                                                                                                                                               tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttt                                                                                                                                                                                                                                                           <  8:1162604/149‑1 (MQ=255)
                                                                                                                                                ttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttta                                                                                                                                                                                                                                                          <  8:843417/149‑1 (MQ=255)
                                                                                                                                                  ggTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                        <  7:821709/149‑1 (MQ=255)
                                                                                                                                                   gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                                                                                                                                                  <  2:180322/138‑1 (MQ=255)
                                                                                                                                                   gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTacac                                                                                                                                                                                                                                                       >  6:334885/1‑149 (MQ=255)
                                                                                                                                                    tAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                                                                                                                                                      >  7:619488/1‑149 (MQ=255)
                                                                                                                                                    tAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                                                                                                                                                      <  6:1328621/149‑1 (MQ=255)
                                                                                                                                                      aGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTgg                                                                                                                                                                                                                                                    <  6:564198/149‑1 (MQ=255)
                                                                                                                                                       gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtt                                                                                                                                                                                                                                                                   >  6:1434523/1‑133 (MQ=255)
                                                                                                                                                            gtctggtcACGCTGGGCGTTGCAGTAGTGATGGTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACGCTGGCGTAcc                                                                                                                                                                                                                                              >  8:932029/1‑149 (MQ=255)
                                                                                                                                                            gtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACTCTGTTTTAGTTTTACACTGGCGTACGTGCATAGATATCAAAGCGATGCACTTTATTGACTACCGCGTTTGGCGTGGCGACATTCGGCAGCGGTGGCGCGTAGTCCGGGAGCTTTAGCACCACGCGTTTGGTCGCACGCAAGTGTGCAgtgt                                                                                                                 >  1:42201/1‑271 (MQ=255)
                                                                                                                                                             tctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                                                                                                                                             >  8:828100/1‑149 (MQ=255)
                                                                                                                                                             tctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                                                                                                                                             <  6:351884/149‑1 (MQ=255)
                                                                                                                                                               tggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAACCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGc                                                                                                                                                                                                                                           >  5:707882/1‑149 (MQ=255)
                                                                                                                                                                 gtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                         >  6:931860/1‑149 (MQ=255)
                                                                                                                                                                  tcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATa                                                                                                                                                                                                                                        >  7:1287499/1‑149 (MQ=255)
                                                                                                                                                                   cACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCTATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAg                                                                                                                                                                                                                                       >  8:330507/1‑149 (MQ=255)
                                                                                                                                                                    aCGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                                                                                                                                                      <  5:852393/149‑1 (MQ=255)
                                                                                                                                                                    aCGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTCAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                                                                                                                                                      >  7:774971/1‑149 (MQ=255)
                                                                                                                                                                      gCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgt                                                                                                                                                                                                                                                                    <  7:510403/117‑1 (MQ=255)
                                                                                                                                                                      gCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGata                                                                                                                                                                                                                                    >  5:536071/1‑149 (MQ=255)
                                                                                                                                                                       gtGGGGGTTGCAGTAGGGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCTGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGGAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCAt                                                                                                     <  4:124130/274‑1 (MQ=255)
                                                                                                                                                                       cTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGatat                                                                                                                                                                                                                                   <  7:831562/149‑1 (MQ=255)
                                                                                                                                                                          ggctgtGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGACGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                            <  2:63914/212‑1 (MQ=255)
                                                                                                                                                                           gCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaa                                                                                                                                                                                                                                >  4:255542/1‑148 (MQ=255)
                                                                                                                                                                           gCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCAGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaaa                                                                                                                                                                                                                               <  6:500393/149‑1 (MQ=255)
                                                                                                                                                                           gCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTCAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaaa                                                                                                                                                                                                                               <  8:774971/149‑1 (MQ=255)
                                                                                                                                                                            cGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAg                                                                                                                                                                                                                              <  7:1087394/149‑1 (MQ=255)
                                                                                                                                                                             gTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGc                                                                                                                                                                                                                             >  5:1098133/1‑149 (MQ=255)
                                                                                                                                                                                gCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGAGAGCAAGGGTGCGGGCTCCAGGAATCCATCGGCAGCaa                                                                                            >  1:37283/1‑275 (MQ=255)
                                                                                                                                                                                  agtagtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGc                                                                                                                                                                                                                        >  5:1222126/1‑149 (MQ=255)
                                                                                                                                                                                  agtagtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGc                                                                                                                                                                                                                        >  8:1215746/1‑149 (MQ=255)
                                                                                                                                                                                    tagtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGccc                                                                                                                                                                                                                      <  6:239974/149‑1 (MQ=255)
                                                                                                                                                                                    tagtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGccc                                                                                                                                                                                                                      <  7:1282674/149‑1 (MQ=255)
                                                                                                                                                                                      gtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCtt                                                                                                                                                                                                                    >  7:232164/1‑149 (MQ=255)
                                                                                                                                                                                         aTGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCCCGCGTTTGGTCGCCAGAACGCGTGCAGGCTCCCGTAATCCATCGGCATCAAAATCCGGTc                                                                                   >  4:327998/1‑275 (MQ=255)
                                                                                                                                                                                                    ggTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTGTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGa                                                                         <  2:175056/274‑1 (MQ=255)
                                                                                                                                                                                                              cTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                                                                                                                                            <  7:1107234/149‑1 (MQ=255)
                                                                                                                                                                                                              cTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                                                                                                                                            >  7:1373392/1‑149 (MQ=255)
                                                                                                                                                                                                               tGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtggc                                                                                                                                                                                           >  7:334749/1‑149 (MQ=255)
                                                                                                                                                                                                               tGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtggc                                                                                                                                                                                           <  5:373327/149‑1 (MQ=255)
                                                                                                                                                                                                                 aaaCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGt                                                                                                                                                                                                                                                                                  <  2:140457/60‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATCGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGAc                                                                                                                                                                                        >  7:1202536/1‑149 (MQ=255)
                                                                                                                                                                                                                  aaCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGAc                                                                                                                                                                                        <  5:617755/149‑1 (MQ=255)
                                                                                                                                                                                                                   aCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACa                                                                                                                                                                                       >  7:143429/1‑149 (MQ=255)
                                                                                                                                                                                                                   aCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACa                                                                                                                                                                                       >  5:299704/1‑149 (MQ=255)
                                                                                                                                                                                                                    cGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAt                                                                                                                                                                                      >  6:1061731/1‑149 (MQ=255)
                                                                                                                                                                                                                     gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                                                                                                                                                     >  8:1402064/1‑149 (MQ=255)
                                                                                                                                                                                                                     gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                                                                                                                                                     >  8:1098986/1‑149 (MQ=255)
                                                                                                                                                                                                                     gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                                                                                                                                                     >  7:1054727/1‑149 (MQ=255)
                                                                                                                                                                                                                      cATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtt                                                                                                                                                                                                                                                             <  5:421674/76‑1 (MQ=255)
                                                                                                                                                                                                                      cATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                                    >  7:941189/1‑149 (MQ=255)
                                                                                                                                                                                                                        tGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGc                                                                                                                                                                                  >  7:396191/1‑149 (MQ=255)
                                                                                                                                                                                                                         gATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGcc                                                                                                                                                                                 >  5:6895/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAg                                                                                                                                                                               >  5:787168/1‑149 (MQ=255)
                                                                                                                                                                                                                               gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                                                                                                                                                           >  6:774616/1‑149 (MQ=255)
                                                                                                                                                                                                                               gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                                                                                                                                                           >  7:101399/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aaaGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTa                                                                                                                                                                   <  8:206147/149‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaGCCATTAACTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTa                                                                                                                                                                   <  5:664971/149‑1 (MQ=255)
                                                                                                                                                                                                                                             attaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCgg                                                                                                                                                             >  8:299757/1‑149 (MQ=255)
                                                                                                                                                                                                                                              ttaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                            >  5:1445728/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  ttattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCt                                                                                                                                                        >  5:386057/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  ttattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTc                                                                                                                <  4:21732/189‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTa                                                                                                                                                     <  7:42235/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                                                                                                                                                    <  7:496219/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                                                                                                                                                    >  6:1145759/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      tCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCCTCAAGATCCGGTCCCACCAGCGACTGAAAGACACGc                                                            >  1:2247/1‑237 (MQ=255)
                                                                                                                                                                                                                                                        ttGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTacca                                                                                                                                                  >  5:189325/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aactaaaaCAGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTaccacc                                                                                                                                                >  8:1011358/5‑149 (MQ=255)
                                                                                                                                                                                                                                                           cAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTAccacca                                                                                                                                               >  5:650720/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGt                                                                                                                                          >  7:1026673/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcgc                                                                                                                                                                      >  6:100650/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                   aGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcgc                                                                                                                                                                      <  5:100650/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc                                                                                                                                    <  8:1287499/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCa                                                                                                                                <  6:787168/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGCCACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTCTGCGGGAACATCGGGTCGAGATaaa  >  1:12447/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                            ttATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGc                                                                                                                              >  7:282385/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
AGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAA  >  NC_000913/3642090‑3642630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑