Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 975,173 G→A 7.1% intergenic (+19/+150) 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,1730GA 7.1% 228.4 / 4.9 80intergenic (+19/+150)ycfJ/suhBflagella biosyntheis regulator/bifunctional inositol monophosphatase/glycerol‑2‑phosphatase
Reads supporting (aligned to +/- strand):  ref base G (40/34);  new base A (4/2);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTACAACGTGACCTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCG  >  NC_002947/975034‑975321
                                                                                                                                           |                                                                                                                                                    
cTACAACGTGACCTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTcgcg                                                                                                                                            >  2:1136995/1‑150 (MQ=255)
 tACAACGTGACCTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCgg                                                                                                                                           <  1:886561/150‑1 (MQ=255)
    aaCGTGACCTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTa                                                                                                                                        >  2:1348417/1‑150 (MQ=255)
     aCGTGACCTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTaa                                                                                                                                       >  1:1010843/1‑150 (MQ=255)
            cTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGAt                                                                                                                                >  2:63683/1‑150 (MQ=255)
                     cGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAggg                                                                                                                       <  2:1361644/150‑1 (MQ=255)
                           atccagccgcacctgccagaaagcatgaagcgctaagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATg                                                                                                                 <  2:1395680‑M1/41‑1 (MQ=255)
                             aGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCg                                                                                                               >  1:109751/1‑150 (MQ=255)
                             aGGTAGGCAAGGTGAAGATGGATAGCGAGCAGGGATCGACCATAACGCTGGAAAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGAAAGTGATTTCCGGAGTCTATGCCGACATCTTCGCGAGTAAACCAGATCCTCCAGAGGCAATGCg                                                                                                               >  2:649597/1‑150 (MQ=255)
                                     aaGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCgg                                                                                                                                           >  2:551869/1‑114 (MQ=255)
                                     aaGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCgg                                                                                                                                           <  1:551869/114‑1 (MQ=255)
                                        tgccagaaagcatgaagcgctaagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGACCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAgg                                                                                                    <  2:1145137‑M1/54‑1 (MQ=255)
                                          gAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCCCGTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGtagtag                                                                                                  >  2:871418/1‑150 (MQ=255)
                                           cagaaagcatgaagcgctaagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGACCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCAATGCGATCCATGt                                                                                                       <  1:555934‑M1/51‑1 (MQ=255)
                                           cagaaagcatgaagcgctaagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGACCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCAATGCGATCCATGt                                                                                                       >  2:555934‑M1/94‑144 (MQ=255)
                                           cagaaagcatgaagcgctaagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGACCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGc                                                                                                 >  2:1299706‑M1/94‑150 (MQ=255)
                                            aGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTaa                                                                                                                                       >  1:606601/1‑111 (MQ=255)
                                            aGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTaa                                                                                                                                       <  2:606601/111‑1 (MQ=255)
                                               tGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGa                                                                                                                                 >  2:1223574/1‑114 (MQ=255)
                                               tGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGa                                                                                                                                 <  1:1223574/114‑1 (MQ=255)
                                                   tCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGCTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATACGATCCATCTAGGCGCGGGTTTAc                                                                                         >  1:315664/1‑150 (MQ=255)
                                                   atgaagcgctaagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTAc                                                                                         >  1:1395680‑M1/86‑150 (MQ=255)
                                                       gAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAAcc                                                                                                                                    <  2:873979/103‑1 (MQ=255)
                                                       gAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAAcc                                                                                                                                    >  1:873979/1‑103 (MQ=255)
                                                           cGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTCAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGCag                                                                                 >  2:441686/1‑150 (MQ=255)
                                                           cGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGaa                                                                                  <  2:859203/149‑1 (MQ=255)
                                                           cGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGaa                                                                                  >  1:859203/1‑149 (MQ=255)
                                                            taagcgcagattgccggcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGACCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAaga                                                                                >  1:1145137‑M1/77‑150 (MQ=255)
                                                            gggATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAaga                                                                                <  2:756374/150‑1 (MQ=255)
                                                                tCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTc                                                                            >  1:199508/1‑150 (MQ=255)
                                                                  gACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCgg                                                                                                                                           >  2:542715/1‑85 (MQ=255)
                                                                  gACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCgg                                                                                                                                           <  1:542715/85‑1 (MQ=255)
                                                                  gACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCgg                                                                          <  1:177427/150‑1 (MQ=255)
                                                                     cATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGttt                                                                                           >  2:401028/1‑130 (MQ=255)
                                                                     cATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGttt                                                                                           <  1:401028/130‑1 (MQ=255)
                                                                     cATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTAc                                                                       >  2:1395357/1‑150 (MQ=255)
                                                                     cATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTAGGTAc                                                                       >  1:1136066/1‑150 (MQ=255)
                                                                       ttCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGaa                                                                                  <  1:817766/137‑1 (MQ=255)
                                                                       ttCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGaa                                                                                  >  2:817766/1‑137 (MQ=255)
                                                                       ttCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAg                                                                     <  2:1010843/150‑1 (MQ=255)
                                                                       ttCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGGGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAg                                                                     <  2:1298193/150‑1 (MQ=255)
                                                                        tCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAgg                                                                    <  2:900461/150‑1 (MQ=255)
                                                                            gcaaacaaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtgCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCAc                                                                >  2:1117331‑M1/61‑150 (MQ=255)
                                                                               gACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCa                                                             >  2:912300/1‑150 (MQ=255)
                                                                                aCAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCaa                                                            >  2:961613/1‑150 (MQ=255)
                                                                                  aaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCgg                                                                                               >  1:1307518‑M1/55‑113 (MQ=255)
                                                                                  aaaaagcacctcttgaggtgcttttttgtgcctggtggttttgcgtggcctgtaCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCgg                                                                                               <  2:1307518‑M1/59‑1 (MQ=255)
                                                                                  aaGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCaaaa                                                          >  1:485726/1‑150 (MQ=255)
                                                                                   aGAAGGCCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTTCGATGCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAac                                                         <  1:87958/150‑1 (MQ=38)
                                                                                    gAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAacc                                                        <  1:961613/150‑1 (MQ=255)
                                                                                          cAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAgg                                                  >  2:344013/1‑150 (MQ=255)
                                                                                          cAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGCGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAgg                                                  <  1:871418/150‑1 (MQ=255)
                                                                                           aaGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAaccac                                                      <  1:928536/145‑1 (MQ=255)
                                                                                           aaGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAaccac                                                      >  2:928536/1‑145 (MQ=255)
                                                                                           aaGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGc                                                 <  1:344013/150‑1 (MQ=255)
                                                                                           aaGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGATGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGc                                                 >  2:1173626/1‑150 (MQ=255)
                                                                                            aGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGca                                                <  1:862112/150‑1 (MQ=255)
                                                                                             gCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGcac                                               >  2:673103/1‑150 (MQ=255)
                                                                                               tGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATACATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACaa                                             >  2:754168/1‑150 (MQ=255)
                                                                                                       aGCGAAGAGCGTACGTGATATGCCGTCTCTGAGCAGGTCTCTACTTCTGCAAATCCTATCCAGCAGGGGCAAGGCGATCCAGGTAGAAGCGGGTTTACCCGCGACGACGTGGGGACATGCCACGCAAAACCACCAGGCGCAAAAAAGCAc                                     >  2:1085973/1‑150 (MQ=255)
                                                                                                        gCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAgg                                                  <  1:55573/136‑1 (MQ=255)
                                                                                                        gCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAgg                                                  >  2:55573/1‑136 (MQ=255)
                                                                                                               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)
                                                                                                                    aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAgg                                                                    >  1:494715/1‑106 (MQ=255)
                                                                                                                    aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAgg                                                                    <  2:494715/106‑1 (MQ=255)
                                                                                                                    aGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCt                        >  2:184157/1‑150 (MQ=255)
                                                                                                                     gTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAg                                                                     <  1:33228/104‑1 (MQ=255)
                                                                                                                     gTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAg                                                                     >  2:33228/1‑104 (MQ=255)
                                                                                                                          ttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGaa                                                                                  >  1:419692/1‑86 (MQ=255)
                                                                                                                          ttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGaa                                                                                  <  2:419692/86‑1 (MQ=255)
                                                                                                                          ttGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTtttgt                  <  1:184157/150‑1 (MQ=255)
                                                                                                                              ggTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCTTGAGGTGCttttttgt                 >  2:902672/1‑144 (MQ=255)
                                                                                                                              ggTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCTTGAGGTGCttttttgt                 <  1:902672/147‑4 (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)
                                                                                                                                           |                                                                                                                                                    
CTACAACGTGACCTATTCCATCGGTGATCAGGTAGGCAAGGTGAAGATGGATCGCGAGCCGGGATCGACCATTCCGCTGGACAAGAATGGCAAGCTGATTCTGAGCGAAGCGGGTCAGTGATTTGCGGTGTCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCTACAGGGGCAATGCGATCCATGTAGGAGCGGGTTTACCCGCGAAGAGGTCGGTACAGGCCACGCAAAACCACCAGGCACAAAAAAGCACCTCAAGAGGTGCTTTTTGTTTGCCGGCAATCTGCG  >  NC_002947/975034‑975321

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