Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,470,108 A→C 5.8% intergenic (‑255/+37) chiA ← / ← tufA periplasmic endochitinase/translation elongation factor EF‑Tu 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,470,1080AC 5.8% 1440.2 / 59.1 417intergenic (‑255/+37)chiA/tufAperiplasmic endochitinase/translation elongation factor EF‑Tu 1
Reads supporting (aligned to +/- strand):  ref base A (202/189);  new base C (9/16);  total (211/206)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 422 total aligned reads displayed.

AGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTT  >  NC_000913/3469960‑3470257
                                                                                                                                                    |                                                                                                                                                      
aGAAAAAAGCCTGACACATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGaa                                                                                                                                                       <  1:6316270/150‑1 (MQ=255)
aGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGaa                                                                                                                                                       <  2:4047255/150‑1 (MQ=255)
  aaaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgc                                                                                                                                                     >  1:5034068/1‑150 (MQ=255)
    aaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgcgc                                                                                                                                                   <  2:841069/150‑1 (MQ=255)
        gCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtt                                                                                                                                               >  1:2636685/1‑150 (MQ=255)
         ccTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCttt                                                                                                                                              >  1:3098290/1‑150 (MQ=255)
          cTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt                                                                                                                                             <  2:1002900/150‑1 (MQ=255)
          cTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTGAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCtttt                                                                                                                                             <  2:1152645/150‑1 (MQ=255)
              cAAATGACGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCGAAAAGGGCGCTTCGGCGCCCtttttgca                                                                                                                                         <  1:2814710/150‑5 (MQ=255)
              cAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaat                                                                                                                                         <  2:103487/150‑1 (MQ=255)
                aaTGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaattc                                                                                                                                       >  2:7144921/1‑150 (MQ=255)
                      gCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACt                                                                                                                                 <  2:2305381/150‑1 (MQ=255)
                       cATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTa                                                                                                                                >  1:2297715/1‑150 (MQ=255)
                          tttAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCTGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAAtt                                                                                                                             >  2:1962672/1‑150 (MQ=255)
                           ttAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGTATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGCAGCGCCCTTTTCAATTCACAAGTAAttg                                                                                                                            >  2:5063943/1‑149 (MQ=255)
                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTATATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgt                                                                                                                       <  1:6409175/150‑1 (MQ=255)
                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgt                                                                                                                       >  2:3737081/1‑150 (MQ=255)
                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgt                                                                                                                       <  2:3572827/150‑1 (MQ=255)
                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgt                                                                                                                       >  2:1217068/1‑150 (MQ=255)
                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgt                                                                                                                       <  1:242277/150‑1 (MQ=255)
                                 aCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAATGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgtgt                                                                                                                     >  1:5586230/1‑150 (MQ=255)
                                 aCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgtgt                                                                                                                     >  1:2958288/1‑150 (MQ=255)
                                 aCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACgtgt                                                                                                                     <  1:1286020/150‑1 (MQ=255)
                                   aGAAACATTCATATTTAGAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTaa                                                                                                                   >  2:1348701/1‑150 (MQ=255)
                                   aGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTaa                                                                                                                   >  1:3625872/1‑150 (MQ=255)
                                     aaaCATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAAtt                                                                                                                 >  2:5288695/1‑150 (MQ=255)
                                        cATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGc                                                                                                              >  2:6869057/1‑150 (MQ=255)
                                            cATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAg                                                                                                          >  2:5523092/1‑150 (MQ=255)
                                              tatTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTACACGT‑TAATTAgttttgaa                                                                                                        <  1:1224720/150‑8 (MQ=255)
                                                tttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACt                                                                                                      <  1:595774/150‑1 (MQ=255)
                                                tttAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACt                                                                                                      >  1:3911143/1‑150 (MQ=255)
                                                  tAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACttt                                                                                                    >  2:6804161/1‑150 (MQ=255)
                                                    aaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAg                                                                                                  >  1:6415176/1‑150 (MQ=255)
                                                     aaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGc                                                                                                 <  2:1734183/150‑1 (MQ=255)
                                                       tGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGcaa                                                                                               <  2:2672317/150‑1 (MQ=255)
                                                        gTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGcaac                                                                                              <  1:3359352/150‑1 (MQ=255)
                                                         ttAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCAAGAACTTAAGcaaca                                                                                             >  2:5875821/1‑150 (MQ=255)
                                                                gaattgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCcg                                                                                      >  1:5738218/1‑150 (MQ=255)
                                                                 aattgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCcgc                                                                                     >  1:538441/1‑150 (MQ=255)
                                                                  attgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCcgcg                                                                                    >  2:2951204/1‑150 (MQ=255)
                                                                             aaaTATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTATTTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACg                                                                         <  1:4758722/150‑1 (MQ=255)
                                                                                  tGAATAATTTATTCGTTCTGACATTACGAATATGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgc                                                                    <  1:4219753/150‑1 (MQ=255)
                                                                                  tGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgc                                                                    <  2:5757113/150‑1 (MQ=255)
                                                                                   gAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgcc                                                                   >  1:2345619/1‑150 (MQ=255)
                                                                                    aataatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCt                                                                  >  2:6474730/1‑150 (MQ=255)
                                                                                       aatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCa                                                               >  1:3213297/1‑150 (MQ=255)
                                                                                        atTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCAc                                                              >  1:2879295/1‑150 (MQ=255)
                                                                                         tttATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACg                                                             >  1:330061/1‑150 (MQ=255)
                                                                                          ttATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACgg                                                            >  2:4232088/1‑150 (MQ=255)
                                                                                            aTTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGAt                                                          <  1:1030403/150‑1 (MQ=255)
                                                                                            aTTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTACACGGAt                                                          >  1:4827617/1‑150 (MQ=255)
                                                                                             ttCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGAtt                                                         >  1:2048910/1‑150 (MQ=255)
                                                                                              tCGTTCTTACATTACTATTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTg                                                        <  2:6685830/150‑1 (MQ=255)
                                                                                              tCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTg                                                        <  2:6197031/150‑1 (MQ=255)
                                                                                              tCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTg                                                        <  2:4769872/150‑1 (MQ=255)
                                                                                                gTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCg                                                      <  1:4455911/150‑1 (MQ=255)
                                                                                                    tGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAAc                                                  >  1:1783492/1‑150 (MQ=255)
                                                                                                         gTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGa                                             >  2:658731/1‑150 (MQ=255)
                                                                                                            cGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcg                                          >  1:806424/1‑150 (MQ=255)
                                                                                                             gAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGACGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCCGAACTTTAGCAACAACGCCCGCGCCAACGGTCCGGCCGCCTTCACGGATTGCGAAACGCAGACcgt                                         >  1:3696407/1‑150 (MQ=255)
                                                                                                             gAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgt                                         >  1:7046548/1‑150 (MQ=255)
                                                                                                                tAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtcgt                                      >  2:3652148/1‑150 (MQ=255)
                                                                                                                 aaGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCgtcgtc                                     >  2:5201905/1‑150 (MQ=255)
                                                                                                                  aGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTcc                                    >  2:6043422/1‑150 (MQ=255)
                                                                                                                    atatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCAt                                  >  2:6968312/1‑150 (MQ=255)
                                                                                                                    atatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCAt                                  >  1:6678472/1‑150 (MQ=255)
                                                                                                                    atatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCAt                                  >  2:1916786/1‑150 (MQ=255)
                                                                                                                     tatGCCGTCAACACATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                                 <  1:4571165/150‑1 (MQ=255)
                                                                                                                     tatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                                 <  1:5070831/150‑1 (MQ=255)
                                                                                                                     tatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                                 <  2:4167283/150‑1 (MQ=255)
                                                                                                                      atGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcg                                >  1:3505539/1‑150 (MQ=255)
                                                                                                                       tGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcgc                               >  2:4841783/1‑150 (MQ=255)
                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCCTCGCGAGCg                          >  1:5839224/1‑150 (MQ=255)
                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCg                          <  1:1569546/150‑1 (MQ=255)
                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCg                          <  1:4096813/150‑1 (MQ=255)
                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCAGCGCGATCg                          >  2:1315616/1‑150 (MQ=255)
                                                                                                                             caacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCgg                         <  2:4543810/150‑1 (MQ=255)
                                                                                                                             caacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCgg                         <  2:1876059/150‑1 (MQ=255)
                                                                                                                              aacaaATGCAAAAAGGGCGCCGAAGCGCCCTGTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCACCGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGCTCggg                        >  2:627070/1‑150 (MQ=255)
                                                                                                                                caaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTg                      <  2:4333744/150‑1 (MQ=255)
                                                                                                                                caaAGGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTg                      <  1:6828361/150‑1 (MQ=255)
                                                                                                                                 aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTgg                     >  2:1624610/1‑150 (MQ=255)
                                                                                                                                  aaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGa                    <  2:2456557/150‑1 (MQ=255)
                                                                                                                                   aTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGAt                   <  2:5598809/150‑1 (MQ=255)
                                                                                                                                      cAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAg                >  2:7036400/1‑150 (MQ=255)
                                                                                                                                       aaaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAgg               <  2:6319400/150‑1 (MQ=255)
                                                                                                                                        aaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAggg              >  1:6240646/1‑150 (MQ=255)
                                                                                                                                         aaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGACACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGt             >  2:383131/1‑150 (MQ=18)
                                                                                                                                         aaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGt             >  1:1532378/1‑150 (MQ=255)
                                                                                                                                         aaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACGTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCAGCGCGATCGGGTGGATCAGGGt             >  2:5690451/1‑150 (MQ=18)
                                                                                                                                           aGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaa           >  1:5316186/1‑150 (MQ=255)
                                                                                                                                            gggCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaac          <  1:4988976/150‑1 (MQ=255)
                                                                                                                                             ggCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaaca         >  2:502579/1‑150 (MQ=255)
                                                                                                                                             ggCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaaca         >  2:4680760/1‑150 (MQ=255)
                                                                                                                                               cgcCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaacaac       <  2:3138442/150‑1 (MQ=255)
                                                                                                                                               cgcCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaacaac       <  1:4367347/150‑1 (MQ=255)
                                                                                                                                                gcCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAAcc      >  2:335955/1‑150 (MQ=255)
                                                                                                                                                  cGAAGCGCCCTTTTCAATTCAAAACTAATTA‑CCGTGTAATTAGCCCAGAACTTTAGCACCCACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCAt    <  2:4902253/150‑1 (MQ=255)
                                                                                                                                                    aaGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCAttt  <  2:1492409/150‑1 (MQ=255)
                                                                                                                                                    aaGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCAttt  <  1:6908808/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTA‑ACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTT  >  NC_000913/3469960‑3470257

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