Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,710,703 G→A 5.5% intergenic (‑228/+96) eptB ← / ← yhjX Kdo2‑lipid A phosphoethanolamine 7''‑transferase/putative pyruvate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,710,7030GA 5.5% 230.2 / 2.0 72intergenic (‑228/+96)eptB/yhjXKdo2‑lipid A phosphoethanolamine 7''‑transferase/putative pyruvate transporter
Reads supporting (aligned to +/- strand):  ref base G (40/28);  new base A (2/2);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATCG  >  NC_000913/3710560‑3710851
                                                                                                                                               |                                                                                                                                                    
taTAAGTAAAATCTACACGCTCGAATAAGTAGACCACAGCAAAGCGAAAAAGTGTGGGGAAAGAGTTTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAtt                                                                                                                                                <  1:1309479/149‑1 (MQ=255)
    aGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGcc                                                                                                                                            <  1:373539/150‑1 (MQ=255)
     gTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCAAGAAGCACAAAATTACAAAAATGCAAGTATGAAAAGAGACAATAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAAATGGcac                                                                                                                                           >  2:1301256/1‑148 (MQ=255)
      tAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCa                                                                                                                                          >  1:595969/1‑150 (MQ=255)
       aaaaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAATCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaa                                                                                                                                         >  2:1032056/1‑150 (MQ=255)
        aaaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCaaa                                                                                                                                        >  2:414409/1‑150 (MQ=255)
         aaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttt                                                                                                                                       <  2:367437/150‑5 (MQ=255)
         aaTCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAATTAACTGGccctttt                                                                                                                                       <  2:667038/150‑5 (MQ=255)
             tACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTa                                                                                                                                   <  1:648722/150‑1 (MQ=255)
              acacGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                                                                  <  2:1412967/150‑1 (MQ=255)
              acacGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTAc                                                                                                                                  <  1:610254/150‑1 (MQ=255)
                          aaGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGttt                                                                                                                      <  1:1173349/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                     <  1:844671/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                     <  1:1619328/150‑1 (MQ=255)
                           aGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTta                                                                                                                     >  2:181649/1‑150 (MQ=255)
                              aGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATAc                                                                                                                  <  1:391378/150‑1 (MQ=255)
                                     cgCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAACAGAGACTTTAGCGGTACAAAAAGCAGTAAAAAGGGCCAGGTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGt                                                                                                           <  2:307579/149‑1 (MQ=255)
                                           gggAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGc                                                                                                     <  2:656068/150‑1 (MQ=255)
                                               aaaaGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGttt                                                                                                 <  1:181649/150‑1 (MQ=255)
                                                aaaGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTg                                                                                                <  1:414409/150‑1 (MQ=255)
                                                     gtgGGGAAAGAGTGTTCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGca                                                                                           <  1:1032056/150‑1 (MQ=255)
                                                     gtgGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGca                                                                                           <  1:955202/150‑1 (MQ=255)
                                                      tgGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGtttcgccagc                                                                                          >  2:877504/1‑143 (MQ=255)
                                                      tgGGGAAAGAGTGTGAATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGAATTTAGCGTTACAAAAAGCATTAAAAAGGGCAAGTTAATTGTCCAAAATTTTACGAGTCATTGTTTATACATTTCGTTTATGAGtttcgacagc                                                                                          >  2:871896/1‑143 (MQ=255)
                                                           aaaGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGa                                                                                     <  2:46747/150‑1 (MQ=255)
                                                               agTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCt                                                                                 >  2:686248/1‑150 (MQ=255)
                                                                    gCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGc                                                                            >  2:995480/1‑150 (MQ=255)
                                                                         aaGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGt                                                                       <  1:1250542/150‑1 (MQ=255)
                                                                           gCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAg                                                                     <  1:1432600/150‑1 (MQ=255)
                                                                           gCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAg                                                                     >  1:13381/1‑150 (MQ=255)
                                                                              caAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTc                                                                  >  1:84942/1‑150 (MQ=255)
                                                                                  tttACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCgg                                                              >  2:983467/1‑150 (MQ=255)
                                                                                   ttACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGc                                                             <  1:941972/150‑1 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                          >  2:608373/1‑150 (MQ=255)
                                                                                      cAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGa                                                          >  2:1241266/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaaa                                                  >  2:750133/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaaa                                                  >  1:223967/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTaaa                                                  >  2:1323008/1‑150 (MQ=255)
                                                                                              cAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGACAGTTAATTGGCCCAAATTTTACGATTCATTGTTTATACATATCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCAAGTAGCTCACGGCTGAGCATTaaa                                                  >  1:181970/1‑150 (MQ=255)
                                                                                                 gTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAggg                                               <  1:1586922/150‑1 (MQ=255)
                                                                                                    tGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGc                                            >  2:46819/1‑150 (MQ=255)
                                                                                                    tGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGc                                            >  1:1196383/1‑150 (MQ=255)
                                                                                                      aaaaGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCa                                          >  2:634228/1‑150 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                      <  1:184488/150‑1 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                      <  1:849282/150‑1 (MQ=255)
                                                                                                          gagaCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATgcg                                      >  2:172068/1‑150 (MQ=255)
                                                                                                             aCTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcca                                   >  1:1148565/1‑149 (MQ=255)
                                                                                                               tttAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCa                                 <  2:1432459/150‑1 (MQ=255)
                                                                                                                 tAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACg                               >  1:157752/1‑150 (MQ=255)
                                                                                                                      ttACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACa                          >  1:709288/1‑150 (MQ=255)
                                                                                                                       tACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACAt                         >  1:549413/1‑150 (MQ=255)
                                                                                                                        aCAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACAtt                        >  1:597429/1‑150 (MQ=255)
                                                                                                                         cAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACAttt                       >  1:149364/1‑150 (MQ=255)
                                                                                                                            aaaGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCt                    <  2:1196383/150‑1 (MQ=255)
                                                                                                                              aGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGc                  >  2:906852/1‑150 (MQ=255)
                                                                                                                               gCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGct                 >  2:935411/1‑150 (MQ=255)
                                                                                                                               gCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGct                 >  1:1526199/1‑150 (MQ=255)
                                                                                                                                  gTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTg              >  2:1481983/1‑150 (MQ=255)
                                                                                                                                   aaatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTgg             >  1:665676/7‑150 (MQ=255)
                                                                                                                                    aatttgggCCAATTAACTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGc            >  2:2221/6‑150 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcc                                    <  1:1082464/126‑1 (MQ=255)
                                                                                                                                    aaaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGcc                                    >  2:1082464/1‑126 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATTAGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt           >  2:219679/1‑150 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGTAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt           <  1:952782/150‑1 (MQ=255)
                                                                                                                                     aaaaGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCt           <  1:686248/150‑1 (MQ=255)
                                                                                                                                        aGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGAc        <  2:149364/150‑1 (MQ=255)
                                                                                                                                          ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGa      >  1:1329206/1‑150 (MQ=255)
                                                                                                                                          ggCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGa      >  1:1093880/1‑150 (MQ=255)
                                                                                                                                            ccAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAt    >  1:1524407/1‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg  >  1:1252202/1‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg  >  2:84887/1‑150 (MQ=255)
                                                                                                                                              aGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAAtcg  >  1:253104/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GATAAGTAAAATCTACACGCGCGAATAAGTAGACCACAGCAAAGGGAAAAAGTGTGGGGAAAGAGTGTGCATGAAGCACAAATTTACAAAAATGCAAGTTTGAAAAGAGACTTTAGCGTTACAAAAAGCAGTAAAAAGGGCCAGTTAATTGGCCCAAATTTTACGAGTCATTGTTTATACATTTCGTTGATGCGTTTGATGCACACGGAAGCTGAAGCCCAGTAGCTCGCGGCTGAGCATTAAAGGGAGCCATGCGCCTCACGCAACATTTTCTGCTCTGGCTGACGAATCG  >  NC_000913/3710560‑3710851

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: