Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 975,291 A→T 6.1% intergenic (+137/+32) ycfJ → / ← suhB flagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947975,2910AT 6.1% 204.6 / 4.3 66intergenic (+137/+32)ycfJ/suhBflagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base T (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

TCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTT  >  NC_002947/975148‑975436
                                                                                                                                               |                                                                                                                                                 
tCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTg                                                                                                                                             >  2:81068/1‑150 (MQ=255)
  aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCCCAAGAGGTGCt                                                                                                                                           >  2:188367/1‑150 (MQ=255)
   gTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCtt                                                                                                                                          <  1:536689/150‑1 (MQ=255)
          gCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTtttgttt                                                                                                                                   <  1:868697/150‑1 (MQ=255)
             gtgtCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGcc                                                                                                                                <  1:646632/150‑1 (MQ=255)
             gtgtCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGcc                                                                                                                                <  2:44431/150‑1 (MQ=255)
             gtgtCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGc                                                                                                                                 <  1:753839/149‑1 (MQ=255)
             gtgtCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGc                                                                                                                                 >  2:753839/1‑149 (MQ=255)
             gcgtCTGTGCGGGCCTTTTCGCGGGTAAACCCGAGCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGGCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGcc                                                                                                                                <  2:155879/148‑1 (MQ=255)
              cgtcgctgCCGTCCTCTTCGCGGGTAAACCCGATCCTACAGGGGTCCTGCGACCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAACGCCCCTCCAGAGGTGCTTTTTGTTTGCCg                                                                                                                               <  2:125771/144‑1 (MQ=37)
                tCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCtt                                                                                                                                          >  1:416866/1‑137 (MQ=255)
                tCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCtt                                                                                                                                          <  2:416866/137‑1 (MQ=255)
                                  cgGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCtt                                                                                                           <  1:408214/150‑1 (MQ=255)
                                     gTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCAAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCAt                                                                                                        >  1:675349/1‑150 (MQ=255)
                                         aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGcc                                                                                                                                <  1:835069/122‑1 (MQ=255)
                                         aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGcc                                                                                                                                >  2:835069/1‑122 (MQ=255)
                                          cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGc                                                                                                                                            >  1:175683/1‑109 (MQ=255)
                                          cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGc                                                                                                                                            <  2:175683/109‑1 (MQ=255)
                                          cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTa                                                                                                                 >  2:94669/1‑136 (MQ=255)
                                          cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTa                                                                                                                 <  1:94669/136‑1 (MQ=255)
                                          cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCttt                                                                                                   >  2:463179/1‑150 (MQ=255)
                                           ccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTc                                                                                                  >  2:566317/1‑150 (MQ=255)
                                           ccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCAGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCTTGAGGTGCttt                                                                                                                                         >  1:60090/1‑111 (MQ=38)
                                           ccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCAGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCTTGAGGTGCttt                                                                                                                                         <  2:60090/111‑1 (MQ=38)
                                            cGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt                                                                                                 <  2:102471/150‑1 (MQ=255)
                                            cGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt                                                                                                 <  1:1125818/150‑1 (MQ=255)
                                             gATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                <  2:99785/150‑1 (MQ=255)
                                                   aCAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGt                                                                                          <  2:420052/150‑1 (MQ=255)
                                                       gggCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCgg                                                                                      >  2:189776/1‑150 (MQ=255)
                                                                   ccATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagt                                                                          >  1:740858/1‑150 (MQ=255)
                                                                     aTGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagtca                                                                        <  1:463179/150‑1 (MQ=255)
                                                                       gTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAACAGGGGCTTTTTGTTTGCCGGCAATCTGCGAATAGCACAACACTCTTGTTGGCAGGTGCGGCTGGATGGcagtcagt                                                                      >  1:804907/1‑150 (MQ=38)
                                                                           gAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCtt                                                                                                                                          <  2:418886/78‑1 (MQ=255)
                                                                           gAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCtt                                                                                                                                          >  1:418886/1‑78 (MQ=255)
                                                                               gtggttACCCGGGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCtt                                                              <  2:125047/146‑1 (MQ=255)
                                                                               gggTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACAAGAAGGTGAAGAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCtt                                                              <  2:870165/150‑1 (MQ=255)
                                                                                      cccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                >  2:149722/1‑109 (MQ=255)
                                                                                      cccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                <  1:149722/109‑1 (MQ=255)
                                                                                      cccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCAc                                                       <  2:1084206/150‑1 (MQ=255)
                                                                                       ccGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACt                                                      >  2:81498/1‑150 (MQ=255)
                                                                                        cgcgAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACtt                                                     <  2:7226/150‑1 (MQ=255)
                                                                                          cgAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGa                                                   <  2:107861/150‑1 (MQ=255)
                                                                                           gAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGAt                                                  <  1:81068/150‑1 (MQ=255)
                                                                                           gAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTCCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGAt                                                  <  1:88222/150‑1 (MQ=255)
                                                                                             agagGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGt                                                <  2:898291/150‑1 (MQ=255)
                                                                                                  tCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGTGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCt                                           >  2:258418/1‑150 (MQ=255)
                                                                                                      tACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCaa                                       >  1:667588/1‑150 (MQ=255)
                                                                                                       aCAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAAc                                      <  1:81498/150‑1 (MQ=255)
                                                                                                         aGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGCCCGGCTGGATGGCAGTACTTATCTCCTTTAAGCACGTGGTGTTGACTGcgacga                                    >  2:290464/1‑150 (MQ=255)
                                                                                                            ccACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAtgt                                 >  1:604837/1‑150 (MQ=255)
                                                                                                              aCGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt                                   <  1:751553/146‑1 (MQ=255)
                                                                                                              aCGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAt                                   >  2:751553/1‑146 (MQ=255)
                                                                                                                  aaaaCCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCtt                           >  1:941793/1‑150 (MQ=255)
                                                                                                                    aaCCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTg                         <  2:931195/150‑1 (MQ=255)
                                                                                                                     accaccAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGt                        <  2:675349/150‑1 (MQ=255)
                                                                                                                      ccaccaGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATATGGccgtcgag                       <  2:667588/150‑7 (MQ=255)
                                                                                                                         ccaGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAg                    <  1:566317/150‑1 (MQ=255)
                                                                                                                           aGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGa                  >  2:646435/1‑150 (MQ=255)
                                                                                                                               acaAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTc              >  2:980739/1‑150 (MQ=255)
                                                                                                                                     aaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAg                                                                                                                <  2:391238/46‑1 (MQ=21)
                                                                                                                                     aaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAg                                                                                                                >  1:391238/1‑46 (MQ=21)
                                                                                                                                           cTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGaa                                                           <  1:280755/93‑1 (MQ=255)
                                                                                                                                           cTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGaa                                                           >  2:280755/1‑93 (MQ=255)
                                                                                                                                           cTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGtt  >  1:330427/1‑150 (MQ=255)
                                                                                                                                             cAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATgg                                                                              >  2:1082381/1‑72 (MQ=255)
                                                                                                                                             cAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATgg                                                                              <  1:1082381/72‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTT  >  NC_002947/975148‑975436

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