Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 788,586 A→T 6.5% intergenic (+33/‑180) ettA → / → gdhA ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947788,5860AT 6.5% 897.7 / 36.3 277intergenic (+33/‑180)ettA/gdhAADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (126/131);  new base T (6/12);  total (133/144)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 282 total aligned reads displayed.

TCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAA  >  NC_002947/788437‑788735
                                                                                                                                                     |                                                                                                                                                     
tCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCAGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCa                                                                                                                                                       >  2:4266139/1‑150 (MQ=255)
  gAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAg                                                                                                                                                     >  1:845478/1‑150 (MQ=255)
  gAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAg                                                                                                                                                     >  1:2306149/1‑150 (MQ=255)
  gAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAg                                                                                                                                                     <  2:4721360/150‑1 (MQ=255)
       tggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGcc                                                                                                                                                >  1:1040039/1‑150 (MQ=255)
       tggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGcc                                                                                                                                                >  1:3336281/1‑150 (MQ=255)
        ggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTCCGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGccc                                                                                                                                               >  2:7445726/1‑150 (MQ=255)
            ttcttcGAAGGTAACTACACCGAGTACGAAGCCGAACGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAtt                                                                                                                                           >  2:2613482/1‑150 (MQ=255)
              cttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCCAGTGGCCCCAtttt                                                                                                                                         >  2:1918912/1‑150 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        <  1:197811/150‑1 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        <  2:5222590/150‑1 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        <  2:6371373/150‑1 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                        >  2:6803351/1‑150 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACTTGGGGCCCCAttctt                                                                                                                                        <  1:5151967/150‑1 (MQ=255)
               ttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCAttctt                                                                                                                                        <  1:306932/150‑1 (MQ=255)
                tcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACTTGGGGCCCCATTCTtg                                                                                                                                       <  1:1642287/150‑1 (MQ=255)
                 cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtgt                                                                                                                                      <  2:809444/150‑1 (MQ=39)
                  gAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTtgtg                                                                                                                                     >  1:2976914/1‑150 (MQ=255)
                      gTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTg                                                                                                                                 <  2:467060/150‑1 (MQ=255)
                      gTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTg                                                                                                                                 >  1:4976333/1‑150 (MQ=255)
                       tAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACTTGGGGCCCCATTCTTGTGcatca                                                                                                                                <  2:623037/150‑5 (MQ=39)
                        aaCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGc                                                                                                                               >  1:4177997/1‑150 (MQ=255)
                           tACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGa                                                                                                                            >  1:6483308/1‑150 (MQ=255)
                           tACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGa                                                                                                                            <  2:549530/150‑1 (MQ=255)
                              acCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCg                                                                                                                         >  2:3680271/1‑150 (MQ=255)
                              acCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCg                                                                                                                         >  2:1040616/1‑150 (MQ=255)
                                cGAGGACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCt                                                                                                                       >  1:285003/1‑150 (MQ=255)
                                  aggaCGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGACACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtat                                                                                                                     <  1:2176462/147‑1 (MQ=255)
                                  aGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtat                                                                                                                     <  1:4630284/150‑1 (MQ=255)
                                   gTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtata                                                                                                                    <  1:1318088/150‑1 (MQ=255)
                                       gAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                                >  2:4840174/1‑150 (MQ=255)
                                       gAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                                >  1:6809823/1‑150 (MQ=255)
                                       gAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATg                                                                                                                >  1:5004441/1‑150 (MQ=255)
                                            cGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgag                                                                                                           <  2:3623853/150‑1 (MQ=255)
                                             gACCGCAAGAAGTGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgagg                                                                                                          >  1:2637096/1‑150 (MQ=255)
                                             gACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgagg                                                                                                          >  1:585511/1‑150 (MQ=255)
                                                cGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggcgaggcga                                                                                                       <  2:5115075/150‑1 (MQ=255)
                                                cGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      <  1:3961502/150‑1 (MQ=255)
                                                 gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      <  2:3923667/150‑1 (MQ=255)
                                                 gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      <  2:2962429/150‑1 (MQ=255)
                                                 gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      <  2:5809668/150‑1 (MQ=255)
                                                 gCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGaa                                                                                                      <  2:500416/150‑1 (MQ=255)
                                                  cAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAg                                                                                                     <  1:7309275/150‑1 (MQ=255)
                                                   aagaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAgg                                                                                                    >  1:4664904/1‑150 (MQ=255)
                                                           tCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                            <  1:6261477/150‑1 (MQ=255)
                                                           tCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTAtttt                                                                                            <  1:6542655/150‑1 (MQ=255)
                                                            cTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTg                                                                                           >  2:6926125/1‑150 (MQ=255)
                                                              gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                         <  2:2324512/150‑1 (MQ=255)
                                                              gggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGta                                                                                         <  2:2888865/150‑1 (MQ=255)
                                                                gTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGtata                                                                                       >  2:3486308/1‑150 (MQ=255)
                                                                     gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACtt                                                                                  >  2:6156574/1‑150 (MQ=255)
                                                                      ctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTa                                                                                 <  1:5705181/150‑1 (MQ=255)
                                                                       tgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTAc                                                                                >  2:7084605/1‑150 (MQ=255)
                                                                        gctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACa                                                                               >  1:4495473/1‑150 (MQ=255)
                                                                         tggcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACAt                                                                              >  2:6746808/3‑150 (MQ=255)
                                                                          tgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATa                                                                             >  1:6957096/1‑150 (MQ=255)
                                                                             ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                          >  1:1701649/1‑150 (MQ=255)
                                                                             ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                          <  1:4074259/150‑1 (MQ=255)
                                                                             ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                          <  1:3214093/150‑1 (MQ=255)
                                                                             ccAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATaaaa                                                                          >  2:660021/1‑150 (MQ=255)
                                                                                      ccatgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCaaa                                                                 <  1:4527320/147‑1 (MQ=255)
                                                                                        gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAaata                                                               >  1:7080476/1‑150 (MQ=255)
                                                                                            aCGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAAttt                                                           >  1:4133649/1‑150 (MQ=255)
                                                                                              ggCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTATTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCa                                                         >  2:6982605/1‑150 (MQ=255)
                                                                                                   aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaaa                                                    >  1:3769720/1‑150 (MQ=255)
                                                                                                    aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaaaa                                                   >  1:1822016/1‑150 (MQ=255)
                                                                                                    aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaaaa                                                   <  1:2366867/150‑1 (MQ=255)
                                                                                                        aCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGAc                                               <  2:1527354/150‑1 (MQ=255)
                                                                                                            gCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACAttt                                           >  1:1102256/1‑150 (MQ=255)
                                                                                                                aGTAAGCCGGTTTTCTGATGCACACGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGcc                                       <  1:5291197/150‑1 (MQ=255)
                                                                                                                     gCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTAtt                                  >  1:1673956/1‑150 (MQ=255)
                                                                                                                      ccGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTc                                 >  2:375317/1‑150 (MQ=255)
                                                                                                                      ccGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTc                                 <  2:4235675/150‑1 (MQ=255)
                                                                                                                        ggTTTTGTGAGGCACAAGATTGGGGCCCCCCGTGGCCCCATTTTTGTGCGTGGCGGACCGCTATAAATGGCGAGGCGATGGTGTATTTTGTATACACTTTCATAAAAGGCACCAAATAATTCCAAAAAAACGACATTTCGCCCTATTCCg                               <  1:5884347/150‑1 (MQ=255)
                                                                                                                          ttttCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGt                             <  2:12643/150‑1 (MQ=255)
                                                                                                                           tttCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTg                            <  2:5369121/150‑1 (MQ=255)
                                                                                                                                gATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAtt                       <  1:3680271/150‑1 (MQ=255)
                                                                                                                                gATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGAtt                       <  1:6065230/150‑1 (MQ=255)
                                                                                                                                  ggCACAAAAATGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGc                     >  2:3958803/2‑150 (MQ=39)
                                                                                                                                     acaAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                  >  2:4033770/1‑150 (MQ=255)
                                                                                                                                     acaAGAATGGGGCCCCAAGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc                  <  1:1554238/150‑1 (MQ=255)
                                                                                                                                      caAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACTACATTTCGCCCTATTCAGGTGCGATTGcttct                 >  2:6307956/1‑150 (MQ=255)
                                                                                                                                       aaaaaTGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt                >  1:3642107/4‑150 (MQ=255)
                                                                                                                                       aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttctt                >  1:948523/1‑150 (MQ=255)
                                                                                                                                        aGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTTTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg               <  2:3212021/150‑1 (MQ=255)
                                                                                                                                        aGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg               <  2:6483308/150‑1 (MQ=255)
                                                                                                                                         gAATGGGGCCCCACGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTgg              <  2:3111258/150‑1 (MQ=255)
                                                                                                                                         gAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTgg              >  2:2101835/1‑150 (MQ=255)
                                                                                                                                          actGGGGCCCCACGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGACGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt             <  2:2439249/148‑1 (MQ=255)
                                                                                                                                          aaTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt             <  1:7377734/150‑1 (MQ=255)
                                                                                                                                          aaTGGGGCCACTTGGGGCCCCATTCTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGt             >  1:3102118/1‑150 (MQ=255)
                                                                                                                                            tGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTAc           >  1:1908677/1‑150 (MQ=255)
                                                                                                                                            tGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTAc           >  2:7314414/1‑150 (MQ=255)
                                                                                                                                              gggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAt         <  2:4525762/150‑1 (MQ=255)
                                                                                                                                              gggCCCCAAGTGGACCCTTTTTTGTGCCTGGCGGACCGATATAAATGGCGAGTCGTAGTTGTATTTTGTATAAACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAt         <  1:4266139/150‑1 (MQ=255)
                                                                                                                                               ggCCCCAGGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACAtt        >  1:2593207/1‑150 (MQ=255)
                                                                                                                                                gCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTc       >  1:6015136/1‑150 (MQ=255)
                                                                                                                                                  cccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCt     >  2:5298109/1‑150 (MQ=255)
                                                                                                                                                   ccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTg    >  2:3153588/1‑150 (MQ=255)
                                                                                                                                                    cAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGa   <  1:3486308/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGAACGTGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGAAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCTGAA  >  NC_002947/788437‑788735

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