Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,724 A→C 6.6% 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 6.6% 254.0 / 7.1 87N105T (AAC→ACC) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base C (3/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869577‑869872
                                                                                                                                                   |                                                                                                                                                    
aaTCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGACTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcg                                                                                                                                                    <  1:516518/150‑1 (MQ=255)
        tACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCgcggcg                                                                                                                                            >  2:1059395/1‑150 (MQ=255)
        tACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCgcggcg                                                                                                                                            <  2:895842/150‑1 (MQ=255)
         aCCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGcggcgg                                                                                                                                           >  2:332875/1‑150 (MQ=255)
         aCCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGcggcgg                                                                                                                                           <  2:55387/150‑1 (MQ=255)
         aCCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGcggcgg                                                                                                                                           >  1:1207076/1‑150 (MQ=255)
          ccAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCgcggcgtt                                                                                                                                          >  2:1524159/1‑148 (MQ=255)
          ccAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCgcggcgtt                                                                                                                                          >  2:778793/1‑148 (MQ=255)
          ccAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGt                                                                                                                                          >  1:573514/1‑150 (MQ=255)
           cAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGGCGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTg                                                                                                                                         <  1:619598/150‑1 (MQ=255)
           cAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTg                                                                                                                                         >  1:1148461/1‑150 (MQ=255)
           cAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTg                                                                                                                                         <  2:1262459/150‑1 (MQ=255)
             ggggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaa                                                                                                                                       <  1:778793/150‑1 (MQ=255)
              gggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGaaa                                                                                                                                      >  1:775599/1‑150 (MQ=255)
               ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAg                                                                                                                                     <  2:1420803/150‑1 (MQ=255)
               ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAg                                                                                                                                     <  2:1411413/150‑1 (MQ=255)
               ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAg                                                                                                                                     >  1:495530/1‑150 (MQ=255)
                gCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGc                                                                                                                                    <  1:1303444/150‑1 (MQ=255)
                   gTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGCTTTCACCGCCGCGGCGTTGAAATCggt                                                                                                                                 <  1:46436/150‑3 (MQ=255)
                   gTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGaa                                                                                                                                 >  1:1314867/1‑150 (MQ=255)
                   gTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGaa                                                                                                                                 <  1:500619/150‑1 (MQ=255)
                   gTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGaa                                                                                                                                 >  2:1121834/1‑150 (MQ=255)
                    ttCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCggtg                                                                                                                                <  1:609916/150‑4 (MQ=255)
                     tCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATc                                                                                                                               <  2:1218184/150‑1 (MQ=255)
                       ggTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTg                                                                                                                             >  2:1312944/1‑150 (MQ=255)
                              aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCggg                                                                                                                      <  2:735561/150‑1 (MQ=255)
                              aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCggg                                                                                                                      >  2:369277/1‑150 (MQ=255)
                               tAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGc                                                                                                                     >  1:1197611/1‑150 (MQ=255)
                                  agagATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGACCAg                                                                                                                  >  1:769192/1‑150 (MQ=255)
                                    agaTGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGATAGCGAATCTGGACCGGGCCAGCg                                                                                                                >  1:92744/1‑150 (MQ=255)
                                       tGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGCCCGGTCCAgattcg                                                                                                             <  1:549256/150‑6 (MQ=255)
                                       tGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATc                                                                                                             >  2:223433/1‑150 (MQ=255)
                                            cTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGTCATTAAATACCAGGATGGCACCGATTTCAACGCCGCGGCGTTGAAAGCGAATGTGGACCGGGCCAGCGATCCGGCg                                                                                                        <  2:1438277/150‑1 (MQ=255)
                                            cTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCg                                                                                                        >  1:1194651/1‑150 (MQ=255)
                                             tGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGa                                                                                                       >  1:912941/1‑150 (MQ=255)
                                              gAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGaa                                                                                                      <  1:223433/150‑1 (MQ=255)
                                              gAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGaa                                                                                                      <  1:367148/150‑1 (MQ=255)
                                                aaaaCGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAatc                                                                                                    >  2:799153/1‑150 (MQ=255)
                                                aaaaCGTGCTGGCGGAGAGTTATACCGATTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAatc                                                                                                    >  2:1464958/1‑150 (MQ=255)
                                                    cGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCt                                                                                                <  1:1437203/150‑1 (MQ=255)
                                                          ggCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACg                                                                                          >  1:898924/1‑150 (MQ=255)
                                                          ggCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACg                                                                                          >  1:610426/1‑150 (MQ=255)
                                                             ggAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCta                                                                                       <  1:43812/150‑1 (MQ=255)
                                                                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATc                                                                                                             <  2:1432105/123‑1 (MQ=255)
                                                                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATc                                                                                                             >  1:1432105/1‑123 (MQ=255)
                                                                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAcc                                                                                  <  1:388349/150‑1 (MQ=255)
                                                                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAAcc                                                                                  >  2:149838/1‑150 (MQ=255)
                                                                  gTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGCCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATATTAAACGCTATAAcc                                                                                  >  1:1160639/1‑150 (MQ=255)
                                                                      taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCTGTTAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGta                                                                              <  1:1102907/150‑1 (MQ=255)
                                                                          gTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAg                                                                          <  2:1387018/150‑1 (MQ=255)
                                                                          gTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAg                                                                          >  1:604694/1‑150 (MQ=255)
                                                                             tCCGCTGAGGTCATTACTTACAGCTTGAAATTGCTGTATGGCATTAAATTACAGGATGGAACCTATTTCAACGCCGCGTCGGTGAAAGCGAATCTGGACCGGCCCAGCGTTCCGGGGAATCATCTTACACGCTATAACCTGTATAAGAat                                                                       <  2:579336/150‑1 (MQ=255)
                                                                              ccGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGATAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAata                                                                      <  2:92744/150‑1 (MQ=255)
                                                                               cGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAatat                                                                     >  2:1114769/1‑150 (MQ=255)
                                                                                    aCGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCCCCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGGGAATCATCTTAAACGGTCTACCCTGTATAAGAATATTGCtt                                                                >  2:1203536/1‑149 (MQ=255)
                                                                                    aCGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTa                                                                <  1:1059395/150‑1 (MQ=255)
                                                                                    aCGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTa                                                                <  1:799153/150‑1 (MQ=255)
                                                                                     cGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaa                                                               <  1:1121834/150‑1 (MQ=255)
                                                                                         aTTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACg                                                           <  1:882021/150‑1 (MQ=255)
                                                                                           tacttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGa                                                         >  2:654075/1‑150 (MQ=255)
                                                                                             cttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAg                                                       <  1:365227/150‑1 (MQ=255)
                                                                                                 cacCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgat                                                   >  1:1149465/1‑150 (MQ=255)
                                                                                                  acCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatc                                                  <  2:1233541/150‑1 (MQ=255)
                                                                                                   ccGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatcg                                                 >  1:640040/1‑150 (MQ=255)
                                                                                                        aaaTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGAGAAAGCGAATCTGGACCGGGACAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATACg                                            >  1:625340/1‑150 (MQ=255)
                                                                                                          aTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGAc                                          <  2:561852/150‑1 (MQ=255)
                                                                                                               gggAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACgg                                     >  2:795411/1‑150 (MQ=255)
                                                                                                                ggAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGt                                    <  2:1389307/150‑1 (MQ=255)
                                                                                                                    ggCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAg                                <  1:369277/150‑1 (MQ=255)
                                                                                                                     gCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAgc                               >  2:1381991/1‑149 (MQ=255)
                                                                                                                     gCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGa                               >  1:1056990/1‑150 (MQ=255)
                                                                                                                     gCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCAGACAACGGTAAAGa                               >  1:1259652/1‑150 (MQ=255)
                                                                                                                      cATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAt                              >  2:1381912/1‑150 (MQ=255)
                                                                                                                       aTTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAtt                             >  2:1221106/1‑150 (MQ=255)
                                                                                                                          aaaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAcc                          <  2:1314867/150‑1 (MQ=255)
                                                                                                                           aaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAccc                         <  1:1470699/150‑1 (MQ=255)
                                                                                                                           aaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAccc                         <  1:1761/150‑1 (MQ=255)
                                                                                                                            aTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCt                        >  2:977079/1‑150 (MQ=255)
                                                                                                                             ttCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTc                       <  2:763950/150‑1 (MQ=255)
                                                                                                                             ttCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTc                       >  2:682291/1‑150 (MQ=255)
                                                                                                                               ccAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCaa                     <  1:1464958/150‑1 (MQ=255)
                                                                                                                               ccAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGACCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCaa                     <  2:769192/150‑1 (MQ=255)
                                                                                                                                      ggCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCg              <  1:1008099/150‑1 (MQ=255)
                                                                                                                                        cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGtt            >  2:196057/1‑150 (MQ=255)
                                                                                                                                           cGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTtctc         <  1:1564013/150‑1 (MQ=255)
                                                                                                                                             aTTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAATAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAg       <  2:912941/150‑1 (MQ=255)
                                                                                                                                             aTTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTGTCAg       <  2:495530/150‑1 (MQ=255)
                                                                                                                                                  aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGttt  <  2:1428975/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869577‑869872

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