Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,062,423 T→G 13.4% intergenic (+81/+100) PP_5311 → / ← PP_5312 signal transduction protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,062,4230TG13.4% 438.8 / 56.2 173intergenic (+81/+100)PP_5311/PP_5312signal transduction protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (79/69);  new base G (9/14);  total (89/84)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

Only 100 of 180 total aligned reads displayed.

GTGAACAAGGTGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACA  >  NC_002947/6062274‑6062572
                                                                                                                                                     |                                                                                                                                                     
gtgAACAAGGTGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCt                                                                                                                                                       >  2:2241172/1‑150 (MQ=255)
    aCAAGGTGCTCGAGGCCTAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCTCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCC                                                                                                                                                   <  1:3430202‑M1/150‑1 (MQ=255)
     cAAGGTGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCC                                                                                                                                                  <  2:197638‑M1/150‑1 (MQ=255)
       aGGTGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGC                                                                                                                                                <  1:1943948‑M1/150‑1 (MQ=255)
        ggTGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCG                                                                                                                                               >  2:2733298‑M1/1‑150 (MQ=255)
          tGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAA                                                                                                                                             <  1:1002972‑M1/150‑1 (MQ=255)
           gCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAac                                                                                                                                            >  2:1018568/1‑150 (MQ=255)
             ttcaggacaagcacatccagatcaagggcaacccggcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGcacgcccgctcccac                                                                                                                                          >  2:3993222‑M1/123‑144 (MQ=255)
              tcaggacaagcacatccagatcaagggcaacccggcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGgcatgcccgctcccaca                                                                                                                                         >  2:1228133‑M1/122‑143 (MQ=255)
              cGAGGCCGAGAACCTGGTGCGTGCGCTGGCCTCGCAGTTCCAGGCACCGCATGAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCAAGCCCGCtcccaca                                                                                                                                         <  2:4871454‑M1/150‑8 (MQ=255)
                aggacaagcacatccagatcaagggcaacccggcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCGTGCCCGCGAACACCgg                                                                                                                                       >  2:2615177‑M1/120‑150 (MQ=255)
                aGGCCGAGAACCTGTTGCGAGAGCTGGCTTCGCAGTTCCATGCACCTCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTACTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGATTGCCCGCGAACACCgg                                                                                                                                       <  2:5388662/150‑1 (MQ=255)
                 ggacaagcacatccagatcaagggcaacccggcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCTTGCCCGCGAACACCGGc                                                                                                                                      >  2:3934228‑M1/119‑150 (MQ=255)
                  gCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCg                                                                                                                                     >  2:806160/1‑150 (MQ=255)
                   ccGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGTAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGa                                                                                                                                    <  2:236316/150‑1 (MQ=255)
                    cGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGaa                                                                                                                                   >  1:2470929/1‑150 (MQ=255)
                    cGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGaa                                                                                                                                   >  1:5306159/1‑150 (MQ=255)
                      agaACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCTCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGc                                                                                                                                 <  2:745422/150‑1 (MQ=255)
                      agaACCTGGTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGc                                                                                                                                 <  2:36301/150‑1 (MQ=255)
                            tGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTgt                                                                                                                           <  1:3910652/150‑2 (MQ=255)
                            tGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATGTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGc                                                                                                                           >  2:3644425/1‑150 (MQ=255)
                              ttGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCa                                                                                                                         <  1:3629160/150‑1 (MQ=255)
                              ttGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCa                                                                                                                         <  1:3629848/150‑1 (MQ=255)
                                 cGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGc                                                                                                                      >  1:3579885/1‑150 (MQ=255)
                                      gggcaacccggcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACcgc                                                                                                                 >  2:5292776‑M1/98‑150 (MQ=255)
                                       cTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATGTCTAACAGCTGGCGCGAGGTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACcgcg                                                                                                                >  2:2859115/1‑150 (MQ=255)
                                         ggCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTc                                                                                                              <  2:2406114/150‑1 (MQ=255)
                                            ttCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCTTCGCCGGTg                                                                                                                            <  2:3829230/133‑1 (MQ=255)
                                            ttCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCTTCGCCGGTg                                                                                                                            >  1:3829230/1‑133 (MQ=255)
                                             tCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCtt                                                                                                          >  1:5367388/1‑150 (MQ=255)
                                             ccggcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCtt                                                                                                          >  1:2869594‑M1/91‑150 (MQ=255)
                                               gCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAAGTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGt                                                                                                        >  2:3697951/1‑150 (MQ=255)
                                                gcgctggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGtt                                                                                                       >  2:1885041‑M1/88‑150 (MQ=255)
                                                 aGTTCCATGCACCGCATTAAGCATTAACACATCGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACCGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAAacaagcgacgcggtgcgtggcaccggcttcgccggtgtt                                                                                                      <  1:2886818‑M1/150‑40 (MQ=255)
                                                  gTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCGTGCCCGCGAAacaagcgacgtggtgcatggcaccggcttcgccggtgttc                                                                                                     <  1:4739876‑M1/150‑41 (MQ=255)
                                                    tggtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTcgc                                                                                                   >  1:4134359‑M1/84‑150 (MQ=255)
                                                     ccATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcg                                                                                                  <  1:4107661‑M1/150‑44 (MQ=255)
                                                     ccATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcg                                                                                                  >  1:181625‑M2/1‑144 (MQ=255)
                                                      gtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCgg                                                                                                 >  1:3944145‑M1/82‑150 (MQ=255)
                                                      gtgatgtggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCgg                                                                                                 >  2:1301137‑M1/82‑150 (MQ=255)
                                                      cATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGACCACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgg                                                                                                 >  2:2090448‑M2/1‑143 (MQ=255)
                                                      cATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCGTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgg                                                                                                 <  1:4956763‑M1/150‑45 (MQ=255)
                                                      cATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCGTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgg                                                                                                 <  1:2416981‑M1/150‑45 (MQ=255)
                                                       aTGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcggg                                                                                                <  1:2259754‑M2/150‑9 (MQ=255)
                                                       aTGAACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCGTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcggg                                                                                                <  2:5100315‑M1/150‑46 (MQ=255)
                                                           aCCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACg                                                                                            >  1:3983898/1‑150 (MQ=255)
                                                           aCCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACg                                                                                            >  2:4963520/1‑150 (MQ=255)
                                                           aCCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACg                                                                                            >  2:4871015/1‑150 (MQ=255)
                                                            tggttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttaaCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGc                                                                                           >  1:5243119‑M1/76‑150 (MQ=255)
                                                              gttccaggggttgatgaaatacctgaaaccaaagaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGccc                                                                                         >  1:1415198‑M1/74‑150 (MQ=255)
                                                              gcattaagcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTACCAGCTGGCGCGATTTTGGACGGTGGGAGCGGGCGTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcgggcaagccc                                                                                         <  1:3653398‑M1/150‑53 (MQ=255)
                                                              gcattaagcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGTGTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGccc                                                                                         >  2:2684141/1‑150 (MQ=255)
                                                               cattaagcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcg                                                                                                  <  2:4146798‑M2/140‑7 (MQ=255)
                                                               cattaagcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcg                                                                                                  <  2:4146860‑M2/140‑7 (MQ=255)
                                                               cattaagcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcg                                                                                                  >  1:4146860‑M2/1‑134 (MQ=255)
                                                                     gcattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTAGCGGGCACGCCCGCTCCca                                                                                  >  1:791876/1‑150 (MQ=255)
                                                                      cattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgc                                                                                                   <  2:1323658‑M1/132‑43 (MQ=255)
                                                                      cattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAAacaagcgacgcggtgcatggcaccggcttcgccggtgttcgc                                                                                                   >  1:1323658‑M1/1‑90 (MQ=255)
                                                                        ttaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAg                                                                               <  1:2260218/150‑1 (MQ=255)
                                                                        ttaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAg                                                                               <  1:2132054/150‑1 (MQ=255)
                                                                         taaCACATAGCCAGCGGGAGCAACTGTCTTGCTCAATATCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCATGCCCGCGAACACCGGCGAACCCGGAGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCACCCACAgg                                                                              >  2:3007691/1‑150 (MQ=37)
                                                                               aTAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGc                                                                        >  1:3504855/1‑150 (MQ=255)
                                                                                   ccTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtca                                                                    >  1:4875064/1‑150 (MQ=255)
                                                                                     tGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtcatc                                                                  >  2:3586023/1‑150 (MQ=255)
                                                                                              agaagaaaggctgaagcccgacggggatgatgacgctttacCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTcgcg                                                                                                  >  2:5229550‑M1/42‑109 (MQ=255)
                                                                                              aaCTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTc                                                         >  2:2377293/1‑150 (MQ=255)
                                                                                                 tGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCggg                                                      <  1:1743584/150‑1 (MQ=255)
                                                                                                 tGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCggg                                                      >  2:433505/1‑150 (MQ=255)
                                                                                                 tGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCggg                                                      >  2:3208198/1‑150 (MQ=255)
                                                                                                     ttGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTc                                                  <  2:2470929/150‑1 (MQ=255)
                                                                                                     ttGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTc                                                  >  2:1424714/1‑150 (MQ=255)
                                                                                                     ttGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCCTCGGGCTTc                                                  >  1:4440202/1‑150 (MQ=255)
                                                                                                        cTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGc                                               >  1:1233085/1‑150 (MQ=255)
                                                                                                        cTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGc                                               <  2:1554757/150‑1 (MQ=255)
                                                                                                          cAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCt                                                                                      >  2:1337479/1‑109 (MQ=255)
                                                                                                                cTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAAcacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcg                                                                                                 >  2:3086591‑M1/1‑48 (MQ=255)
                                                                                                                cTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAAcacaagcgacgcggtgcatggcaccggcttcgccggtgttcgcg                                                                                                 <  1:3086591‑M1/92‑45 (MQ=255)
                                                                                                                  aaCCGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTcttctt                                     <  1:4331276/150‑1 (MQ=255)
                                                                                                                  aaCAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTc                                                         <  2:384423/130‑1 (MQ=255)
                                                                                                                  aaCAGCTGGCGCGAGTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTcttctt                                     <  1:5034893/150‑1 (MQ=255)
                                                                                                                     aGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTgg                                  <  1:420032/150‑1 (MQ=255)
                                                                                                                      gCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAATCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGt                                 >  1:5261686/1‑150 (MQ=255)
                                                                                                                          gcgcGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCa                             <  2:4072242/150‑1 (MQ=255)
                                                                                                                             cgATTTTGGACTGTGGGAGCGGGCTTTCCCGCGAACACCGGCGAAGCCGGTGCCAAGCACCGCGTCGCTTTTTTAGCGGGCACCCCCGCTCCCACAGGTAAAGCGTCATCAAACCCCACGGGCTTCAGCCTTTCTTCTTTGGTTTCACGt                          >  1:2803656/1‑150 (MQ=25)
                                                                                                                               aTTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCGCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGcc                                              >  2:17529/1‑128 (MQ=255)
                                                                                                                               aTTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCCCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGcc                                              <  1:17529/128‑1 (MQ=255)
                                                                                                                                  ttGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtc                     <  1:2241172/150‑1 (MQ=255)
                                                                                                                                             GAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTg          >  2:283412‑M1/1‑150 (MQ=255)
                                                                                                                                               GCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCGTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGa        <  1:3200174‑M1/150‑1 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccaccgtccaaaatcgcgccagctggtagaaattgagcaagacagttgctcccgcaggctatgtgttaatgcttaatgc      <  1:3653398‑M2/150‑99 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagtttctcccgcaggctatgtgttaatgcttaatgc      <  2:422894‑M2/150‑99 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagttgctcccgccggctatgtgttaatgcttaatgc      <  2:3786199‑M2/150‑99 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagttgctcccgcaggctatgtgttaattcttaatgc      <  1:3653474‑M2/150‑99 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagttgctcccgcaggctatgtgttaatgcttaatgc      <  2:3786302‑M2/150‑99 (MQ=255)
                                                                                                                                                 gggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagttgctcccgcaggctatgtgttaatgcttaatgc      <  2:255592‑M2/150‑99 (MQ=255)
                                                                                                                                                 GGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAAc      <  2:4443046‑M1/150‑1 (MQ=255)
                                                                                                                                                  ggCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTttcgcgggcacgcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagttgctcccgcaggctatgtgttaatgcttaatgcg     <  2:4873365‑M2/150‑100 (MQ=255)
                                                                                                                                                   gCTTGCCCGCGAAcaccggcgaagccggtgccatgcaccgcgtcgccggtgttcgcgggcaagcccgctcccacagtccaaaatcgcgccagctgttagaaattgagcaagacagttgctcccgcaggctatgtgttaatgcttaatgcg    >  2:2723556‑M1/1‑13 (MQ=255)
                                                                                                                                                    CTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACcac   >  2:3564352‑M1/1‑150 (MQ=255)
                                                                                                                                                     ttGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCACCCACAGGTAAAGCGTCATCATACCCGTCGGGCTTCAGCCTTTCTTATTTGGTTTCAGGTATTTAATAAACCCCTGGAACcaca  >  1:2879488/1‑150 (MQ=37)
                                                                                                                                                     |                                                                                                                                                     
GTGAACAAGGTGCTCGAGGCCGAGAACCTGTTGCGTGAGCTGGCTTCGCAGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACA  >  NC_002947/6062274‑6062572

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