Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,724 A→C 7.9% N105T (AAC→ACC)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,7240AC 7.9% 225.4 / 9.1 89N105T (AAC→ACC) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (42/39);  new base C (2/5);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

GAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869575‑869872
                                                                                                                                                     |                                                                                                                                                    
gAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCaa                                                                                                                                                      >  2:815859/1‑150 (MQ=255)
 aaaTCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTAAACAGTGAAATTGCGGGAAGGCATTAAATTCAAGGATGTCACCGATTTCaaa                                                                                                                                                     >  1:50602/1‑149 (MQ=255)
  aaTCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcc                                                                                                                                                    >  2:351182/1‑149 (MQ=255)
   aTCGTTTTACCAGGGGCTGTTCTGTCTGGATAAATAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCTGGATGGCACCGATTTCAAcgc                                                                                                                                                   <  1:116729/150‑1 (MQ=255)
    tCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgcc                                                                                                                                                  <  1:1294384/150‑1 (MQ=255)
      gTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgccgc                                                                                                                                                >  2:1485714/1‑150 (MQ=255)
      gTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGAGGAGAGTTATACCGTTTCAGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgccgc                                                                                                                                                >  1:494318/1‑150 (MQ=255)
           aCCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGcggcgg                                                                                                                                           <  1:537296/150‑1 (MQ=255)
             cAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGGGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTg                                                                                                                                         <  2:457631/150‑1 (MQ=255)
             aaGGGGCTGTTCGGTCTGGATAAAGAGATGAAAATGAAAAACGTGCTGGCGTAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAATTTCCAGGATGGCACCGATTTCAACGCCGAGGCGGTg                                                                                                                                         <  2:266810/149‑1 (MQ=255)
                gggCTGTTCGGGCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGaaa                                                                                                                                      <  1:809986/150‑1 (MQ=255)
                 gtcggTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCGAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaat                                                                                                                                     <  2:603199/146‑2 (MQ=255)
                 ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAg                                                                                                                                     <  1:928998/150‑1 (MQ=255)
                  gCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGc                                                                                                                                    <  1:1410334/150‑1 (MQ=255)
                  gCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGc                                                                                                                                    >  1:715983/1‑150 (MQ=255)
                   cTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTAAACGCCGCGGCGTTGAAAGCg                                                                                                                                   >  1:624802/1‑150 (MQ=255)
                   cTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACAGTTTCGGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCg                                                                                                                                   <  2:222407/150‑1 (MQ=255)
                   cTGTGCGGGCTGGATAAAGAGAGGAAACTGAAATTCGTGCGGGCGGAGAGTTATTCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGTAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGGTGAAATCg                                                                                                                                   <  2:286883/150‑1 (MQ=255)
                     gtgcGGTCTGGATAAAGCGATGAAACAGAAAAACTTTCTGGCGGATAGTAATACCGATTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGACGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGGTGAAAGCGaa                                                                                                                                 <  1:850516/147‑1 (MQ=255)
                     gctCGGTCTGGATAACGAGATGCAACTGAAAAACGTGCTGGAGGAGAGTTATACCGGTTCCGATGACGGCAGTACTAACACCGTGACCTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGaa                                                                                                                                 <  1:269738/148‑1 (MQ=255)
                     gTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGaa                                                                                                                                 >  1:1165670/1‑150 (MQ=255)
                      ttCGGTCTGGATAAAGAGATGAAACTGAATAACGGGCTGGCGGAGATGTATACCGTTTGCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAAt                                                                                                                                <  2:822265/150‑1 (MQ=255)
                      gtCGGTCTGGATAAAGAGAGGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGCGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCggtg                                                                                                                                <  1:1300795/149‑4 (MQ=255)
                          gTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTgg                                                                                                                            <  2:379652/150‑1 (MQ=255)
                          gTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTgg                                                                                                                            <  1:1021084/150‑1 (MQ=255)
                          gTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTgg                                                                                                                            <  1:1220098/150‑1 (MQ=255)
                                  aaaGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCCGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGcc                                                                                                                    >  2:1105691/1‑150 (MQ=255)
                                  aaaGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGcc                                                                                                                    >  1:578958/1‑150 (MQ=255)
                                  aaaGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGcc                                                                                                                    >  2:1124799/1‑150 (MQ=255)
                                    agatATCAAACTGAAAAACGTGCTGGCGGAGAGCTATCCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAg                                                                                                                  <  1:1354077/150‑1 (MQ=255)
                                    agagATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAg                                                                                                                  >  1:9339/1‑150 (MQ=255)
                                     gagaTGAAACTGAAAAACGTGCTGGCGTAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGCGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGc                                                                                                                 <  1:323656/150‑1 (MQ=255)
                                          gAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATcc                                                                                                            <  1:1195638/150‑1 (MQ=255)
                                           aaaCTGAAAAACGTGCTGGCGGAGGGTTATACCGTTTCCGATGAGGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCCGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCg                                                                                                           >  2:1570844/1‑150 (MQ=255)
                                              cTGAAAAACGAGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGAGGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCg                                                                                                        <  1:1064168/150‑1 (MQ=255)
                                                     aCGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACCCCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAatcatc                                                                                                 >  1:1046701/1‑150 (MQ=255)
                                                       gTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCtt                                                                                               <  1:714029/150‑1 (MQ=255)
                                                        tGCTGGCGGAGAGTTATACCGTTTCCTATGACGGCATTACTTACACCGTGAAATTGCTGGAAGGCATTAAATTCCAGGATGGCACCTATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTa                                                                                              <  1:823/150‑1 (MQ=255)
                                                        tGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTa                                                                                              >  1:1207935/1‑150 (MQ=255)
                                                           gggCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGATTCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAAc                                                                                           <  1:1625260/149‑1 (MQ=255)
                                                               ggAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCta                                                                                       <  2:729849/150‑1 (MQ=255)
                                                                   agTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAc                                                                                   <  1:556090/150‑1 (MQ=255)
                                                                      tataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTg                                                                                >  2:71292/1‑150 (MQ=255)
                                                                      tataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCAGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATACCCTg                                                                                >  1:92711/1‑150 (MQ=255)
                                                                       ataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               >  1:727012/1‑150 (MQ=255)
                                                                       ataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               >  2:1058682/1‑150 (MQ=255)
                                                                        taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAACAGATATAACCTGtt                                                                              >  2:185505/1‑149 (MQ=255)
                                                                         aCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtat                                                                             >  2:1364680/1‑150 (MQ=255)
                                                                          ccGTTTCCGCTGACGGCATATCTTACACCGTGATATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtata                                                                            <  1:1428867/150‑1 (MQ=255)
                                                                          ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtata                                                                            >  1:567529/1‑150 (MQ=255)
                                                                                ccGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAata                                                                      <  2:9339/150‑1 (MQ=255)
                                                                                 cGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAatat                                                                     <  2:1184902/150‑1 (MQ=255)
                                                                                 cGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAatat                                                                     >  1:333416/1‑150 (MQ=255)
                                                                                 cGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACAGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAatat                                                                     >  1:996878/1‑150 (MQ=255)
                                                                                   aTGACTGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGGATCATCTTAAACGCTATAACCTGTATAAGAATATTg                                                                   <  2:494318/150‑1 (MQ=255)
                                                                                    tGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGc                                                                  >  1:920835/1‑150 (MQ=255)
                                                                                       ctgCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaa                                                               <  2:576529/148‑1 (MQ=255)
                                                                                       cGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaa                                                               >  1:1666771/1‑150 (MQ=255)
                                                                                         gCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa                                                             >  2:1434602/1‑150 (MQ=255)
                                                                                          cATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAAc                                                            >  2:118034/1‑150 (MQ=255)
                                                                                             tacttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACggc                                                         >  1:166972/1‑149 (MQ=255)
                                                                                             tacttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGa                                                         >  1:1447677/1‑150 (MQ=255)
                                                                                             tacttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGa                                                         <  1:1407015/150‑1 (MQ=255)
                                                                                              acttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGaa                                                        >  1:1162344/1‑150 (MQ=255)
                                                                                              acttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGaa                                                        >  1:471259/1‑150 (MQ=255)
                                                                                              acttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAAAGCTATAACCTGTATAAGAATATTGCTAAAACGGaa                                                        >  2:1505456/1‑150 (MQ=255)
                                                                                                 tacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGCTCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcg                                                     >  2:1252625/1‑150 (MQ=255)
                                                                                                 tacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcg                                                     >  1:522212/1‑150 (MQ=255)
                                                                                                    acCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatc                                                  <  1:1485714/150‑1 (MQ=255)
                                                                                                       gTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatcgat                                               >  2:1715789/1‑150 (MQ=255)
                                                                                                       gTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatcgat                                               >  1:654215/1‑150 (MQ=255)
                                                                                                        tGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCgatcgatc                                              <  1:1313790/150‑1 (MQ=255)
                                                                                                          aaaTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCg                                            <  2:652528/150‑1 (MQ=255)
                                                                                                            aTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGACGCGGCGGAGAAAGCGAATCTGGACCGGGCCAGCTATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAAAATTACTAAAACGGAAGCGATCGATCCGAc                                          >  2:1144657/1‑150 (MQ=255)
                                                                                                                cGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCTGACAACg                                      <  1:71292/150‑1 (MQ=255)
                                                                                                                     aGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTaaa                                 >  1:1420507/1‑150 (MQ=255)
                                                                                                                         aTTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAtt                             <  2:1296219/150‑1 (MQ=255)
                                                                                                                         aTTAAATTACAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGTGAATATAGAACGGGCCAGCCATCCGGCGAATCATCTTAAACGATATTACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCAGACAACGGTAAAGAtt                             >  1:1097924/1‑150 (MQ=255)
                                                                                                                              aTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCt                        <  1:828146/150‑1 (MQ=255)
                                                                                                                                tCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCa                      <  2:1447677/150‑1 (MQ=255)
                                                                                                                                   aGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAAc                   <  1:1364680/150‑1 (MQ=255)
                                                                                                                                        ggCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCg              <  1:624429/150‑1 (MQ=255)
                                                                                                                                        ggCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCg              >  2:407467/1‑150 (MQ=255)
                                                                                                                                         gCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGCCAACGGTAAAGATTACCCTCAAACAGCCGt             >  2:1370905/1‑150 (MQ=255)
                                                                                                                                            tcgctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAACCGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGGtct          >  1:1063617/5‑150 (MQ=255)
                                                                                                                                             cGATTTCATCGCCGCGGCGGTGAAAGCGAATCTGGACCTGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTtctc         <  2:558914/150‑1 (MQ=255)
                                                                                                                                              gctttCACCGCCGCGGCGTTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGGAAAGATTACCCTCAAACAGCCGTTCTCa        >  1:180031/3‑150 (MQ=255)
                                                                                                                                              gATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCa        >  1:828659/1‑150 (MQ=255)
                                                                                                                                                 ttCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCg     <  2:578958/150‑1 (MQ=255)
                                                                                                                                                   cAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGCTCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGtt   >  1:694153/1‑150 (MQ=255)
                                                                                                                                                   cAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGtt   <  2:727012/150‑1 (MQ=255)
                                                                                                                                                    aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGttt  <  2:764469/150‑1 (MQ=255)
                                                                                                                                                    aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAACCGCAAACAGCCGTTCTCAGCGttt  >  2:1191313/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869575‑869872

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