Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 975,291 A→T 10.2% 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,2910AT10.2% 257.5 / 18.5 88intergenic (+137/+32)ycfJ/suhBflagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base T (6/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

GGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAA  >  NC_002947/975145‑975439
                                                                                                                                                  |                                                                                                                                                    
gggTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAagag                                                                                                                                                   <  1:754168/150‑1 (MQ=255)
  gTCAGTGATTTGCGGTGTTAGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGt                                                                                                                                                 <  2:938964/150‑1 (MQ=255)
   tCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTg                                                                                                                                                >  1:198381/1‑150 (MQ=255)
   tCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTg                                                                                                                                                >  2:1235323/1‑150 (MQ=255)
     aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCt                                                                                                                                              >  2:184157/1‑150 (MQ=255)
           ttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTtttgt                                                                                                                                        <  1:184157/150‑1 (MQ=255)
               gtggcctgtaccgacctcttcgcgggtaaacccgctcctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggccggCACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGc                                                                                                                                    >  2:1117331‑M2/91‑150 (MQ=255)
               ggTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCTTGAGGTGCttttttgt                                                                                                                                       <  1:902672/147‑4 (MQ=255)
               ggTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCTTGAGGTGCttttttgt                                                                                                                                       >  2:902672/1‑144 (MQ=255)
               ggTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGc                                                                                                                                    <  1:1235323/150‑1 (MQ=255)
                gtgtCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGcc                                                                                                                                   >  1:1082780/1‑150 (MQ=255)
                 tgtCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTCCCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCg                                                                                                                                  <  2:982423/150‑1 (MQ=255)
                           GGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgcg                                                                                                                        <  1:190164‑M1/150‑1 (MQ=255)
                              ctctTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCAAGCAAAACCACAAGGCAAAACAAAACACCTCAAGAGGCGCTTTTTGTTTTCCGGCAATCTGCGCtt                                                                                                                     >  2:1156073/1‑150 (MQ=38)
                               tcttcgcgggtaaacccgctcctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggtcggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTa                                                                                                                    >  1:1145137‑M2/75‑150 (MQ=255)
                                 ttCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAaga                                                                                                                                                    >  1:1138922/1‑116 (MQ=255)
                                 ttCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAaga                                                                                                                                                    <  2:1138922/116‑1 (MQ=255)
                                        gtaaacccgctcctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggccggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCAt                                                                                                           >  1:1395680‑M2/66‑150 (MQ=255)
                                        gTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCAt                                                                                                           <  1:123018/150‑1 (MQ=255)
                                        gTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCAt                                                                                                           <  1:912300/150‑1 (MQ=255)
                                         tAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGAGGTTCGAGCTGCATg                                                                                                          >  2:7891/1‑150 (MQ=39)
                                          aaaCCGGATACTACAGGGGCAAGGCGAGCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGc                                                                                                         <  2:109751/150‑1 (MQ=255)
                                            aCCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCtt                                                                                                       >  1:262771/1‑150 (MQ=255)
                                             cccGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCttt                                                                                                      >  1:100414/1‑150 (MQ=255)
                                              ccgctcctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggccggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCttttt                                                                                                                                          >  1:1307518‑M2/60‑113 (MQ=255)
                                              ccgctcctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggccggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCttttt                                                                                                                                          <  2:1307518‑M2/54‑1 (MQ=255)
                                               cGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCt                                                                                                    <  1:1395357/150‑1 (MQ=255)
                                                gctcctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggtcggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                   >  2:1299706‑M2/58‑150 (MQ=255)
                                                 aTCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAagag                                                                                                                                                   <  2:617519/101‑1 (MQ=255)
                                                 aTCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAagag                                                                                                                                                   >  1:617519/1‑101 (MQ=255)
                                                   cctacatggatcgcattgcccctgtaggatcgggtttacccgcgaagaggtcggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCa                                                                                                <  2:1145137‑M2/96‑1 (MQ=255)
                                                     tACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAgg                                                                                              >  2:17949/1‑150 (MQ=255)
                                                      acatggatcgcattgcccctgtaggagcgggtttacccgcgaagaggtcggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                   >  2:555934‑M2/52‑144 (MQ=255)
                                                      acatggatcgcattgcccctgtaggagcgggtttacccgcgaagaggtcggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTg                                                                                                   <  1:555934‑M2/93‑1 (MQ=255)
                                                        aGGGGCAATGCGATCCATGTAGGAGCGGGCTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGc                                                                                           >  1:520820/1‑150 (MQ=255)
                                                            gCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCt                                                                                       <  2:315664/150‑1 (MQ=255)
                                                                tGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGCATACGCGAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGAt                                                                                   <  1:1156073/150‑1 (MQ=255)
                                                                cattgcccctgtaggatcgggtttacccgcgaagaggccggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGAt                                                                                   <  2:1395680‑M2/109‑1 (MQ=255)
                                                                 gCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCCGTAATGGCCACGCAAAACCACCAGGCACAAAACAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATg                                                                                  <  1:673103/150‑1 (MQ=255)
                                                                  cGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATgg                                                                                 <  1:1348417/150‑1 (MQ=255)
                                                                  cGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATgg                                                                                 <  2:485726/150‑1 (MQ=255)
                                                                    aTCCATGTAGGAGCGGGTTTACCAGCGAAGAGGTCAGTACAAGCCACGCAAAACCACCAGGCACAAAAACGCACCTCAAGAGGGGCTTGTTGTCTGCCGTCAATCGGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGca                                                                               <  1:441686/150‑1 (MQ=255)
                                                                      ccATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCGTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagt                                                                             >  1:449786/1‑150 (MQ=255)
                                                                       cATGTATGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagtc                                                                            <  1:649597/150‑1 (MQ=255)
                                                                           tAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTa                                                                        <  1:63683/150‑1 (MQ=255)
                                                                             ggAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTAcc                                                                      >  1:432577/1‑150 (MQ=255)
                                                                                gCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCgg                                                                                                                                 <  2:260181/88‑1 (MQ=255)
                                                                                gCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCAACAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCgg                                                                                                                                 >  1:260181/1‑88 (MQ=255)
                                                                                  gggTTTACTCTCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCtt                                                                 <  1:1136995/150‑1 (MQ=255)
                                                                                     tttACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATTCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGaa                                                              <  2:199508/150‑1 (MQ=255)
                                                                                     tttACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGaa                                                              <  2:1082780/150‑1 (MQ=255)
                                                                                      ttACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAg                                                             <  2:520820/150‑1 (MQ=255)
                                                                                           cgcgAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACtt                                                        <  2:198381/150‑1 (MQ=255)
                                                                                              gAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGAt                                                     >  1:999697/1‑150 (MQ=255)
                                                                                              gAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGAt                                                     <  2:432577/150‑1 (MQ=255)
                                                                                              gAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACATCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGAt                                                     >  1:968757/1‑150 (MQ=255)
                                                                                                agagGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATgg                                                                                 >  1:178591/1‑120 (MQ=255)
                                                                                                agagGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATgg                                                                                 <  2:178591/120‑1 (MQ=255)
                                                                                                    gTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc                                               <  2:190164/150‑1 (MQ=255)
                                                                                                    gTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc                                               >  2:781503/1‑150 (MQ=255)
                                                                                                    gTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGATTTTTGTTTGCCGGCAATCATTTCGTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGcc                                               >  1:174234/1‑150 (MQ=255)
                                                                                                      cggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAg                                                                                                                   >  1:733749‑M2/4‑80 (MQ=255)
                                                                                                      cggTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAg                                                                                                                   <  2:733749‑M2/77‑1 (MQ=255)
                                                                                                        gTACAGGCAACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTTCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTTAAGCACTTCACGTTGCCTGCa                                           >  2:508409‑M2/2‑150 (MQ=255)
                                                                                                         tACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCtt                                                                 >  1:173375/1‑127 (MQ=255)
                                                                                                         tACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCtt                                                                 <  2:173375/127‑1 (MQ=255)
                                                                                                            aGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCCAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGa                                       >  1:1152603/1‑150 (MQ=255)
                                                                                                             ggCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgcgc                                                                                                                       <  1:750417/69‑1 (MQ=255)
                                                                                                             ggCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTgcgc                                                                                                                       >  2:750417/1‑69 (MQ=255)
                                                                                                              tctaggCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAtg                                     >  1:49421/6‑150 (MQ=255)
                                                                                                              gCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCttt                                                                                                      >  1:1042217/1‑85 (MQ=255)
                                                                                                              gCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCttt                                                                                                      <  2:1042217/85‑1 (MQ=255)
                                                                                                              gCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAtg                                     >  2:1370377/1‑150 (MQ=255)
                                                                                                              gCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAtg                                     >  2:12415/1‑150 (MQ=255)
                                                                                                               ccACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGAtgt                                    >  2:1387490/1‑150 (MQ=255)
                                                                                                                cACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGc                                                >  1:1097389/1‑137 (MQ=255)
                                                                                                                cACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGc                                                <  2:1097389/137‑1 (MQ=255)
                                                                                                                       aaCCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTg                            >  1:1014527/1‑150 (MQ=255)
                                                                                                                           accaGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGa                        <  1:1299706/150‑1 (MQ=255)
                                                                                                                            ccaGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAg                       <  1:1117331/150‑1 (MQ=255)
                                                                                                                                gCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAg                   <  2:1136066/150‑1 (MQ=255)
                                                                                                                                    aaaaaaGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGt               <  2:174234/150‑1 (MQ=255)
                                                                                                                                     aaaaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagtcagt                                                                         <  1:915374/91‑1 (MQ=255)
                                                                                                                                     aaaaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGcagtcagt                                                                         >  2:915374/1‑91 (MQ=255)
                                                                                                                                     aaaaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACtt                                                        <  2:458412/108‑1 (MQ=255)
                                                                                                                                     aaaaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACtt                                                        >  1:458412/1‑108 (MQ=255)
                                                                                                                                      aaaaGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTgg             <  2:968757/150‑1 (MQ=255)
                                                                                                                                       aaaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGc            >  1:522012/1‑150 (MQ=255)
                                                                                                                                        aaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGCTGCATCTCAGTACCGCCGTGAAGCACTTGATGTTGAATACAACGCTCTGGCCTCTTTCGAGGAAGTCGTGGcc           >  1:386378/1‑150 (MQ=255)
                                                                                                                                        aaGCACCTCTTGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTAACGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGGAGTCGTGGcc           >  2:266438/1‑150 (MQ=255)
                                                                                                                                            aCCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTTCAACGATGTGGCCTTTGTCGa                        >  2:270811/1‑133 (MQ=255)
                                                                                                                                            aCCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGa                        <  1:270811/133‑1 (MQ=255)
                                                                                                                                               tCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTg    >  1:669876/1‑150 (MQ=255)
                                                                                                                                                 aaGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGaa  <  1:12415/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCGCTTAGCGCTTCATGCTTTCTGGCAGGTGCGGCTGGATGGCAGTCAGTACCGCCTTGAAGCACTTGATGTTGCCTGCAACGATGTGGCCTTTGTCGAGGAAGTCGTGGCCACCGTTGAA  >  NC_002947/975145‑975439

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: