Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 928,207 T→A 7.8% Q330L (CAA→CTA)  cydC ← glutathione/cysteine ABC transporter export permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913928,2070TA 7.8% 250.7 / 8.8 77Q330L (CAA→CTA) cydCglutathione/cysteine ABC transporter export permease/ATPase
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base A (3/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAG  >  NC_000913/928059‑928352
                                                                                                                                                    |                                                                                                                                                 
tttGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGttt                                                                                                                                                   >  2:301842/1‑149 (MQ=255)
tttGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGttt                                                                                                                                                   <  1:419567/149‑1 (MQ=255)
  tGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTgg                                                                                                                                                 <  3:65270/149‑1 (MQ=255)
        aTCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAAGCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCa                                                                                                                                           <  2:415154/149‑1 (MQ=255)
         tCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTgg                                                                                                                                                 >  2:402201/1‑142 (MQ=255)
         tCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTgg                                                                                                                                                 <  1:402201/142‑1 (MQ=255)
            ggTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGaa                                                                                                                                       <  4:137226/149‑1 (MQ=255)
                 gCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCGCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTg                                                                                                                                  <  1:395520/149‑1 (MQ=255)
                  cccGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGa                                                                                                                                 >  3:396724/1‑149 (MQ=255)
                           cGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGtt                                                                                                                        <  3:408890/149‑1 (MQ=255)
                            gCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGttt                                                                                                                       >  1:6475/1‑149 (MQ=255)
                             cTATATGTTCCCCGGCGTTTACCTGAAGATAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGtttt                                                                                                                      <  2:72242/149‑1 (MQ=255)
                              tataTGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAAATTGCTCCGGATGAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACTCGAGTTTGGGTATCAGGAAAGGAGACCTCCGGttttt                                                                                                                     >  4:411259/1‑149 (MQ=255)
                                tatGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGa                                                                                                                   <  1:470251/149‑1 (MQ=255)
                                tatGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGa                                                                                                                   <  3:563047/149‑1 (MQ=255)
                                    ttCCCCGGCGTTTACCTTAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGATATCGGCAACACGTGTTTGGGTATCAGTAAAGGTCACCTCCGGTTTTTGATCCg                                                                                                               <  4:348862/149‑1 (MQ=255)
                                       cccGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTa                                                                                                            <  2:6475/149‑1 (MQ=255)
                                        ccGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCATAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                                                                                           <  4:371413/149‑1 (MQ=255)
                                         cGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAg                                                                                                          >  2:475762/1‑149 (MQ=255)
                                         cGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAg                                                                                                          >  2:80220/1‑149 (MQ=255)
                                             gTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCag                                                                                                      >  3:27351/1‑149 (MQ=255)
                                              tttACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCaga                                                                                                     <  4:318693/149‑1 (MQ=255)
                                              tttACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCaga                                                                                                     <  1:470001/149‑1 (MQ=255)
                                                 aacTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAt                                                                                                  <  2:37522/147‑1 (MQ=255)
                                                     gAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAtacgt                                                                                              >  3:140408/1‑149 (MQ=255)
                                                     gAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAt                                                                                                  >  1:326641/1‑145 (MQ=255)
                                                     gAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAt                                                                                                  <  2:326641/145‑1 (MQ=255)
                                                       agagAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAtacgtac                                                                                            >  2:489094/1‑149 (MQ=255)
                                                       agagAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGAACCGTTAAGTCAGAGAtacgtac                                                                                            >  4:466570/1‑149 (MQ=255)
                                                             aTCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCaga                                                                                      <  1:190165/149‑1 (MQ=255)
                                                             aTCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCaga                                                                                      >  2:234344/1‑149 (MQ=255)
                                                                 cTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTg                                                                                                                                                  <  1:44352/85‑1 (MQ=255)
                                                                 cTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTg                                                                                                                                                  >  2:44352/1‑85 (MQ=255)
                                                                   ttAAGTGCCTGTTGAGATTTCTCCGGATAAGAGAACTGAACAAACCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAAt                                                                                <  2:308172/149‑1 (MQ=255)
                                                                        tGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACtt                                                                           <  1:471231/149‑1 (MQ=255)
                                                                         gCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCAGTACCGTCAGCGAAACACGATCTGCAACACGAGTCTGGTTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGACAGAGATACGTACGGCTGAGGCAATGACTTg                                                                          >  3:350318/1‑149 (MQ=255)
                                                                          ccTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc                                                                         >  4:67851/1‑149 (MQ=255)
                                                                           cTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGcc                                                                        <  3:466570/149‑1 (MQ=255)
                                                                              ttGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCa                                                                     >  4:417778/1‑149 (MQ=255)
                                                                                        tCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGaaa                                                           <  3:426755/149‑1 (MQ=255)
                                                                                         ccGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAAt                                                          <  4:396724/149‑1 (MQ=255)
                                                                                         ccGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAAt                                                          <  1:85255/149‑1 (MQ=255)
                                                                                            gATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCa                                                       <  1:301842/149‑1 (MQ=255)
                                                                                                 gTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                  >  4:387192/1‑149 (MQ=255)
                                                                                                 gTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                  >  4:54036/1‑149 (MQ=255)
                                                                                                 gTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                  >  1:406343/1‑149 (MQ=255)
                                                                                                 gTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                  >  4:270923/1‑149 (MQ=255)
                                                                                                 gTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                  >  1:104296/1‑149 (MQ=255)
                                                                                                 gTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                  >  2:66655/1‑149 (MQ=255)
                                                                                                  tgaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGtt                                                 >  2:501312/1‑149 (MQ=255)
                                                                                                      ctgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg                                             >  2:256400/1‑149 (MQ=255)
                                                                                                            aTCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCa                                       >  3:133421/1‑149 (MQ=255)
                                                                                                              cccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGt                                     >  2:482966/1‑149 (MQ=255)
                                                                                                               ccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTg                                    >  1:529312/1‑149 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                  >  1:235474/1‑149 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTAGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                   >  4:338465/1‑148 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTAGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                   >  4:114000/1‑148 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTAGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                   <  3:114000/148‑1 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTAGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                   >  4:415171/1‑148 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTAGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                   <  3:415171/148‑1 (MQ=255)
                                                                                                                 gTAACGTCAGCGAAACGCGATCGGCAACACGAGTTAGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                   <  3:338465/148‑1 (MQ=255)
                                                                                                                      gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCaaa                             <  2:93465/149‑1 (MQ=255)
                                                                                                                         aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAcgc                          >  3:537943/1‑149 (MQ=255)
                                                                                                                          gCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAcgcg                         <  4:133421/149‑1 (MQ=255)
                                                                                                                          gCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAcgcg                         >  3:146101/1‑149 (MQ=255)
                                                                                                                                cgcgATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAAc                   <  2:452962/149‑1 (MQ=255)
                                                                                                                                   gATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTAc                                               >  2:236672/1‑118 (MQ=255)
                                                                                                                                   gATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTAc                                               <  1:236672/118‑1 (MQ=255)
                                                                                                                                   gATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAAcgcg                <  3:73491/149‑1 (MQ=255)
                                                                                                                                   gATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAAcgcg                >  3:240543/1‑149 (MQ=255)
                                                                                                                                    aTCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACgcgc               >  1:122511/1‑149 (MQ=255)
                                                                                                                                     tCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCa              <  3:67851/149‑1 (MQ=255)
                                                                                                                                      cGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCag             <  4:537943/149‑1 (MQ=255)
                                                                                                                                          aaCACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCagaaga         <  4:140408/149‑1 (MQ=255)
                                                                                                                                           acacGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGAc        <  1:482966/149‑1 (MQ=255)
                                                                                                                                               gAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAAc    <  1:234344/149‑1 (MQ=255)
                                                                                                                                                 gTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAg  <  1:80220/149‑1 (MQ=255)
                                                                                                                                                 gTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAg  <  1:475762/149‑1 (MQ=255)
                                                                                                                                                 gTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAg  >  4:88210/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTTGCCGCATCCGGTTCGCCCGAGAATCGCTATATGTTCCCCGGCGTTTACCTGAAGAGAAATCCCTTTAAGTGCCTGTTGAGATTGCTCCGGATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAG  >  NC_000913/928059‑928352

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