Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,915,978 C→T 6.9% L308L (CTG→TTG) ‡ barA → sensory histidine kinase BarA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,915,9780CT 6.9% 268.3 / 6.7 87L308L (CTG→TTG) ‡barAsensory histidine kinase BarA
Reads supporting (aligned to +/- strand):  ref base C (36/45);  new base T (4/2);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

GACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCA  >  NC_000913/2915831‑2916124
                                                                                                                                                   |                                                                                                                                                  
gACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTg                                                                                                                                                  >  1:866850/1‑150 (MQ=255)
   cAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAAt                                                                                                                                               <  1:255153/150‑1 (MQ=255)
   cAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAAt                                                                                                                                               >  2:1390959/1‑150 (MQ=255)
     ggCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACt                                                                                                                                                   <  2:1312522/144‑1 (MQ=255)
     ggCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACt                                                                                                                                                   >  1:1312522/1‑144 (MQ=255)
         aCGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGtt                                                                                                                                         <  1:1270082/150‑1 (MQ=255)
          cGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGgtgt                                                                                                                                          >  1:1717379/1‑148 (MQ=255)
          cGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGgtgt                                                                                                                                          <  2:1717379/148‑1 (MQ=255)
          cGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTa                                                                                                                                        <  1:1048450/150‑1 (MQ=255)
           gTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTAt                                                                                                                                       <  2:1407521/150‑1 (MQ=255)
           gTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTAt                                                                                                                                       <  1:1065224/150‑1 (MQ=255)
            tCCGATCTGCGTGAAACGCTGGAGCAGTTTGAAATTCAGAACGTTGAGTTTGTTCTGGCGAAAAAGCGCGCGCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTGACACGAGCTGCGTACACCACTGATTGGTGTTAtt                                                                                                                                      <  2:399791/150‑1 (MQ=255)
            tCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTAtt                                                                                                                                      >  1:298369/1‑150 (MQ=255)
            tCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTAtt                                                                                                                                      >  1:1806248/1‑150 (MQ=255)
             ccGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCGCTGAATGGTGTTATTg                                                                                                                                     <  2:552501/150‑1 (MQ=255)
             ccGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTg                                                                                                                                     <  1:1840382/150‑1 (MQ=255)
                  cTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCttt                                                                                                                                <  2:1611433/150‑1 (MQ=255)
                   tGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCATTCAGTGgtgtacgcagc                                                                                                                                   >  2:92296/1‑139 (MQ=255)
                   tGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCATTCAGTGgtgtacgcagc                                                                                                                                   <  1:92296/146‑8 (MQ=255)
                   tGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTa                                                                                                                               <  1:793489/150‑1 (MQ=255)
                       tGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCg                                                                                                                           >  2:846990/1‑150 (MQ=255)
                       tGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCg                                                                                                                           >  2:1528619/1‑150 (MQ=255)
                       tGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCg                                                                                                                           >  2:386818/1‑150 (MQ=255)
                         aaaCGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGcc                                                                                                                         <  1:93715/150‑1 (MQ=255)
                              cTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACg                                                                                                                    <  2:1545669/150‑1 (MQ=255)
                                ggAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCt                                                                                                                  <  2:455849/150‑1 (MQ=255)
                                ggAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCt                                                                                                                  >  2:1460552/1‑150 (MQ=255)
                                 gAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTg                                                                                                                 <  2:161642/150‑1 (MQ=255)
                                   gCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGaa                                                                                                               <  1:12323/150‑1 (MQ=255)
                                     aGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGaaaa                                                                                                             >  2:1618775/1‑150 (MQ=255)
                                        tGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAg                                                                                                          <  1:1208416/150‑1 (MQ=255)
                                          gAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGaa                                                                                                        >  2:1294630/1‑150 (MQ=255)
                                          gAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGaa                                                                                                        >  1:1430872/1‑150 (MQ=255)
                                          gAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGaa                                                                                                        <  1:1390959/150‑1 (MQ=255)
                                           aaaTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCATTCAGTGgtgtacg                                                                                                                                       <  2:1812361/118‑4 (MQ=255)
                                           aaaTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCATTCAGTGgtgtacg                                                                                                                                       >  1:1812361/1‑115 (MQ=255)
                                           aaaTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAAt                                                                                                       <  1:1248643/150‑1 (MQ=255)
                                             aTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTa                                                                                                     <  1:1528619/150‑1 (MQ=255)
                                             aTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTa                                                                                                     <  1:1742555/150‑1 (MQ=255)
                                               tCAGAACGTTGAGTTTGATCTGTCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTTAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAac                                                                                                   <  1:1075985/150‑1 (MQ=255)
                                                cAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAaca                                                                                                  <  2:779458/150‑1 (MQ=255)
                                                 aGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAacac                                                                                                 <  2:1143031/150‑1 (MQ=255)
                                                  gAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCTAGTTTCTGGCAAATATGTCACACGAGCTGCGGACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAACCAGAATTAACAcc                                                                                                <  2:1761736/150‑1 (MQ=255)
                                                   aaCGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCa                                                                                               <  2:1161952/150‑1 (MQ=255)
                                                    aCGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCaa                                                                                              <  1:14470/150‑1 (MQ=255)
                                                     cGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAAc                                                                                             >  2:1077946/1‑150 (MQ=255)
                                                         gAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAg                                                                                         <  2:779036/150‑1 (MQ=255)
                                                              aGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGa                                                                                    >  1:472245/1‑150 (MQ=255)
                                                                aTCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATc                                                                                  <  2:1416814/150‑1 (MQ=255)
                                                                  cTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTAACCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCAc                                                                                >  2:1264124/1‑150 (MQ=255)
                                                                     gCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTg                                                                             <  1:1460552/150‑1 (MQ=255)
                                                                     gCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTg                                                                             <  1:1294630/150‑1 (MQ=255)
                                                                      cGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTGTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGa                                                                            >  1:1394590/1‑150 (MQ=255)
                                                                       gAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGaa                                                                           >  2:54973/1‑150 (MQ=255)
                                                                             gcgcgcCCAGGAAGCGGCGCGTATTAAATGCGAGTTTCTGGCAAATATGTCAGACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCTCCTGACGCTGAAAACAGAATTAACACCAACGCATCGCGATCACCTGATTACGAt                                                                     <  2:95577/150‑1 (MQ=255)
                                                                                      ggAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTc                                                            >  1:288780/1‑150 (MQ=255)
                                                                                       gAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCg                                                           <  1:1726853/150‑1 (MQ=255)
                                                                                       gAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCg                                                           >  2:718896/1‑150 (MQ=255)
                                                                                         aGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGc                                                         <  2:866850/150‑1 (MQ=255)
                                                                                                cgTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACAcc                                                                                                <  1:1816233/104‑1 (MQ=255)
                                                                                                cgTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACAcc                                                                                                >  2:1816233/1‑104 (MQ=255)
                                                                                                    ttAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGc                                              <  1:846990/150‑1 (MQ=255)
                                                                                                     aaaaTCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAAAAAtttcat                                             >  2:1803585/2‑147 (MQ=255)
                                                                                                      aaaTCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTg                                            <  1:274303/150‑1 (MQ=255)
                                                                                                       aaTCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTgg                                           >  1:733304/1‑150 (MQ=255)
                                                                                                         tCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCa                                         >  2:883279/1‑150 (MQ=255)
                                                                                                          ccGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCaa                                        >  1:1927403/1‑150 (MQ=255)
                                                                                                              gTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAattat                                    >  1:1681464/1‑150 (MQ=255)
                                                                                                               tttCTGTCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAattatt                                   <  1:834145/150‑1 (MQ=255)
                                                                                                                   tGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatg                               <  1:1618775/150‑1 (MQ=255)
                                                                                                                   tGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatg                               >  2:697781/1‑150 (MQ=255)
                                                                                                                   tGGCAAATATGTCACACGAGCTGCGTACACAACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatg                               >  1:743004/1‑150 (MQ=255)
                                                                                                                     gCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatgat                             >  2:147535/1‑150 (MQ=255)
                                                                                                                        aaTATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTAGTAATGATGtt                          >  1:360140/1‑150 (MQ=255)
                                                                                                                             gTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGa                     >  1:1945932/1‑150 (MQ=255)
                                                                                                                              tCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAg                                                                                                          <  2:1384283/64‑1 (MQ=255)
                                                                                                                              tCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAg                                                                                                          >  1:1384283/1‑64 (MQ=255)
                                                                                                                               cacaCGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACt                   >  2:1497443/1‑150 (MQ=255)
                                                                                                                                   cGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTtctc               <  2:1712248/150‑1 (MQ=255)
                                                                                                                                   cGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTtctc               <  2:1719253/150‑1 (MQ=255)
                                                                                                                                    gAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCg              >  1:35450/1‑150 (MQ=255)
                                                                                                                                    gAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCg              <  1:1957200/150‑1 (MQ=255)
                                                                                                                                     aGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGa             >  2:562201/1‑150 (MQ=255)
                                                                                                                                      gccaataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGaa            >  1:1955201/8‑150 (MQ=255)
                                                                                                                                      gccaataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGaa            >  1:121587/8‑150 (MQ=255)
                                                                                                                                        tGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAAc          >  1:1877958/1‑150 (MQ=255)
                                                                                                                                           gTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTgg       <  2:1193320/150‑1 (MQ=255)
                                                                                                                                              caccacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAg    >  2:1266137/1‑150 (MQ=255)
                                                                                                                                               accacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGc   <  1:555428/150‑1 (MQ=255)
                                                                                                                                                ccacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCa  >  1:934772/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GACCAGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCA  >  NC_000913/2915831‑2916124

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