Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 541,524:1 +A 6.2% coding (960/1302 nt) ybbY → putative uracil/xanthine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913541,5241.A 6.2% 250.6 / 12.2 65coding (960/1302 nt)ybbYputative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base . (27/34);  new base A (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGACCTTCCATTAACGCTTCGTCCTCTGCTCAGTAACGGCTTA  >  NC_000913/541242‑541803
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                        
gCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTACTTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCTCTTACACGTGGGCTCCTGttatttat                                                                                                                                                                                                                                                                                    >  2:672881/1‑289 (MQ=255)
  aGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGTGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATg                                                                                                                                                                                                                                                                                   >  1:316284/1‑288 (MQ=255)
   gTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGc                                                                                                                                                                                                                                                                                 <  2:362038/289‑1 (MQ=255)
      gTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAg                                                                                                                                                                                                                                                                               <  2:504543/288‑1 (MQ=255)
          ccGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTa                                                                                                                                                                                                                                                                          >  1:327038/1‑289 (MQ=255)
                       gCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCtgg                                                                                                                                                                                                                                                              <  2:327038/288‑1 (MQ=255)
                            gCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGTGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGgcg                                                                                                                                                                                                                                                         <  2:316284/288‑1 (MQ=255)
                              tttGTGGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTTGTAAATATCAGCAATACCTACGGTGCGATTCTGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCt                                                                                                                                                                                                                                                       <  1:746469/288‑1 (MQ=255)
                               ttGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTgg                                                                                                                                                                                                                                                     >  1:290410/1‑289 (MQ=255)
                                                tCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGa                                                                                                                                                                                                                                    <  1:257742/289‑1 (MQ=255)
                                                         ggTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGTCAGCGTTATTTTCCTGCTGGTGGCGCTGGTTCCTGCACTCCAGCGACTGTTTTg                                                                                                                                                                                                                            >  2:470638/1‑288 (MQ=255)
                                                            gATTACAGGTCTGGTAAATATCATCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGt                                                                                                                                                                                                                        <  2:686405/289‑1 (MQ=255)
                                                             aTTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGt                                                                                                                                                                                                                        <  1:534345/288‑1 (MQ=255)
                                                                       tGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTAATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTa                                                                                                                                                                                                               >  1:667252/1‑288 (MQ=255)
                                                                         gTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCg                                                                                                                                                                                                           >  1:209817/1‑289 (MQ=255)
                                                                         gTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCg                                                                                                                                                                                                           >  1:196367/1‑289 (MQ=255)
                                                                           aaaTATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGt                                                                                                                                                                                                          <  2:424758/288‑1 (MQ=255)
                                                                           aaaTATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGt                                                                                                                                                                                                          <  1:164115/288‑1 (MQ=255)
                                                                             atatCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGag                                                                                                                                                                                                       <  2:290410/289‑1 (MQ=255)
                                                                                   gCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGt                                                                                                                                                                                                                        <  1:329663/266‑1 (MQ=255)
                                                                                   gCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGt                                                                                                                                                                                                                        >  2:329663/1‑266 (MQ=255)
                                                                                   gCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGc                                                                                                                                                                                                 >  1:143948/1‑289 (MQ=255)
                                                                                    cAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTTCTGGTGGCGCTGGTTCCTTCCTTCACGCGAATGTTTTGCAGTATCCCTTTACCCGTGagtagt                                                                                                                                                                                                   >  2:708967/1‑286 (MQ=255)
                                                                                     aaTACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTg                                                                                                                                                                                                  <  1:708967/286‑1 (MQ=255)
                                                                                          cTACGGAGCGTGTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACCATGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATg                                                                                                                                                                                          <  1:420522/289‑1 (MQ=255)
                                                                                            aCGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCa                                                                                                                                                                                                                          >  2:637327/1‑255 (MQ=255)
                                                                                            aCGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCa                                                                                                                                                                                                                          <  1:637327/255‑1 (MQ=255)
                                                                                              ggTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTgg                                                                                                                                                                                      >  2:622265/1‑289 (MQ=255)
                                                                                              ggTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTgg                                                                                                                                                                                      >  2:202501/1‑289 (MQ=255)
                                                                                                tGCGATTCGGGGCACAGATGTTGTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGt                                                                                                                                                                                     <  1:672881/288‑1 (MQ=255)
                                                                                                 gCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTAATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGt                                                                                                                                                                                     <  2:667252/288‑1 (MQ=255)
                                                                                                  cGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGtt                                                                                                                                                                                    >  1:80693/1‑287 (MQ=255)
                                                                                                  cGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGGTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGtt                                                                                                                                                                                    >  1:355012/1‑287 (MQ=255)
                                                                                                  cGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCCGGTGGCGCTGGTTCCTGCACTCACGCGCATGGTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCCgggt                                                                                                                                                                                    >  2:553063/1‑285 (MQ=255)
                                                                                                       cGGGGCACGGATGGTTTTTAACCGAAGAAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAGCCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCTATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTTACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATc                                                                                                                                                                             <  2:231436/289‑1 (MQ=255)
                                                                                                                         tATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTcc                                                                                                                                                            >  2:628459/1‑288 (MQ=255)
                                                                                                                         tATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTcc                                                                                                                                                            >  2:392712/1‑288 (MQ=255)
                                                                                                                            ccGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCCGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGCCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCcgc                                                                                                                                                          >  1:29266/1‑287 (MQ=255)
                                                                                                                             cgcgtcacgggtaaGGGAATACGCGTTATCGTGGTATCTTTGTGGGGACCGGATTCAGTCCGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCt                                                                                                                                                       <  2:265343/276‑1 (MQ=255)
                                                                                                                                agcagGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGgt                                                                                                                                                    >  1:243232/1‑289 (MQ=255)
                                                                                                                                     ggCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTAATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCtt                                                                                                                                                                                                                 <  1:374020/224‑1 (MQ=255)
                                                                                                                                     ggCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTAATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCtt                                                                                                                                                                                                                 >  2:374020/1‑224 (MQ=255)
                                                                                                                                               aTACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAAt                                                                                                                                     <  1:628459/289‑1 (MQ=255)
                                                                                                                                                        aTCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCt                                                                                                                                                       <  2:578481/262‑1 (MQ=255)
                                                                                                                                                        aTCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCt                                                                                                                                                       >  1:578481/1‑262 (MQ=255)
                                                                                                                                                             cgtAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCt                                                                                                                        <  2:196367/288‑1 (MQ=255)
                                                                                                                                                              gtAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCg                                                                                                                      <  1:470638/289‑1 (MQ=255)
                                                                                                                                                               tAGCTTTGTGGCGACCGGATTTATGAAGCTGAGAACCGTACAGCTGGCGGTAATTCCATTTTCACCGTTTGTGTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCg                                                                                                                      <  2:355012/288‑1 (MQ=255)
                                                                                                                                                                            aCCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGa                                                                                                         >  2:300151/1‑288 (MQ=255)
                                                                                                                                                                                 aTTTATGGAGCTGATAACCGTACCGCTGGCGGTAATTCCATTTGTACCGTTTGTTTCTTCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCTGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGGGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGAc                                                                                                        <  2:209817/284‑1 (MQ=255)
                                                                                                                                                                                  tttATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCa                                                                                                   <  1:392712/288‑1 (MQ=255)
                                                                                                                                                                                                   cGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGc                                                                                 <  1:622265/289‑1 (MQ=255)
                                                                                                                                                                                                    gTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCa                                                                                <  2:143948/289‑1 (MQ=255)
                                                                                                                                                                                                     tACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTAATCGACTCGCATTGCCGTTATTTGTCGGCat                                                                               >  1:631743/1‑289 (MQ=255)
                                                                                                                                                                                                             gCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCACTCTTTACCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTaa                                                                       >  2:472210/1‑289 (MQ=255)
                                                                                                                                                                                                                    ttCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCAtt                                                                >  2:238435/1‑289 (MQ=255)
                                                                                                                                                                                                                               ccGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTgtgt                                                     <  1:273951/289‑1 (MQ=255)
                                                                                                                                                                                                                                      tttCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCTACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCa                                              <  1:553063/289‑1 (MQ=255)
                                                                                                                                                                                                                                            ccATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGACCt                                        <  2:243232/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        taACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGACCTTCCATTAACGCt                            <  1:472210/289‑1 (MQ=255)
                                                                                                                                                                                                                                                               gACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCtgt                                                       >  1:595402/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                 cTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGACCTTCCATTAACGCTTCGTCct                     <  1:300151/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGACCTTCCATTAACGCTTCGTCct                     <  1:238435/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCt                                                                                                                                                                                                                                                       >  2:584054/1‑58 (MQ=255)
                                                                                                                                                                                                                                                                    gCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCt                                                                                                                                                                                                                                                       <  1:584054/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gcgTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGCCCTTCCCTTAACGCTTCGTCCCTGTCTCAGTAACGGCttt  >  2:673894/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                        
GCAGTGGTTCCCGCTCGGCAGTGGCGGTGCTTTGTCGCCGGGAATTATTCTGACGGCGGTGATTACAGGTCTGGTAAATATCAGCAATACCTACGGTGCGATTCGGGGCACGGATGTTTTTTATCCGCAGCAGGGCGCAGGGAATACGCGTTATCGTCGTAGCTTTGTGGCGACCGGATTTATGACGCTGATAACCGTACCGCTGGCGGTAATTCCATTTTCACCGTTTGTTTCATCCATTGGTTTATTAACCCAGACTGGCGATTACACGCGGCGTTCGTTTATTTATGGCAGCGTTATTTGCCTGCTGGTGGCGCTGGTTCCTGCACTCACGCGACTGTTTTGCAGTATCCCTTTACCCGTGAGTAGTGCGGTCATGCTGGTTTCTTATCTGCCTTTACTCTTTTCCGCGCTGGTGTTTAGCCAGCAAATAACGTTTACCGCTCGCAATATTTATCGACTCGCATTGCCGTTATTTGTCGGCATATTTTTAATGGCATTACCGCCTGTGTATCTGCAAGACCTTCCATTAACGCTTCGTCCTCTGCTCAGTAACGGCTTA  >  NC_000913/541242‑541803

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