Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,962,078 A→T 5.0% intergenic (+41/‑155) PP_2588 → / → PP_2589 class III aminotransferase/aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,962,0780AT 5.0% 975.3 / 26.8 304intergenic (+41/‑155)PP_2588/PP_2589class III aminotransferase/aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base A (135/142);  new base T (2/13);  total (143/161)
Fisher's exact test for biased strand distribution p-value = 7.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 305 total aligned reads displayed.

ATGATCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAG  >  NC_002947/2961930‑2962227
                                                                                                                                                    |                                                                                                                                                     
atgatCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCAc                                                                                                                                                      <  2:6214814/150‑1 (MQ=255)
    tCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGccc                                                                                                                                                  >  1:1281017/1‑150 (MQ=255)
     cATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCcg                                                                                                                                                 >  1:7445733/1‑150 (MQ=255)
     cATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCcg                                                                                                                                                 >  2:4968599/1‑150 (MQ=255)
       tGGCCCCGGCCTTGATCGACAACCACAGCGAGATGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGTGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCcgcg                                                                                                                                               >  2:6474794/1‑150 (MQ=255)
         gCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGaa                                                                                                                                             >  1:4260331/1‑150 (MQ=255)
         gCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGaa                                                                                                                                             <  2:5116340/150‑1 (MQ=255)
         caacgagaccgacgtcctgaccgacgctcagcgcaactatgccgaacagcagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGCACgcccgctcc                                                                                                                                             <  1:4871971‑M2/20‑4 (MQ=255)
              tccgaggcgcgttgttgccagaattcgagagtgtacatttgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCg                                                                                                                                        >  1:680427‑M1/121‑150 (MQ=255)
               gCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCgg                                                                                                                                       <  1:1616995/150‑1 (MQ=255)
                ccTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGc                                                                                                                                      <  2:7413303/150‑1 (MQ=255)
                  ttGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGa                                                                                                                                    <  2:3545422/150‑1 (MQ=255)
                   tGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGaa                                                                                                                                   <  1:3551346/150‑1 (MQ=255)
                    gATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAg                                                                                                                                  >  2:1444843/1‑150 (MQ=255)
                       cGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCg                                                                                                                               <  1:991054/150‑1 (MQ=255)
                        gCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGTGGGCACGCCCGCGAACACCGGCGAAGCCgg                                                                                                                              >  1:2168008/1‑150 (MQ=255)
                        gCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCgg                                                                                                                              >  1:7339596/1‑150 (MQ=255)
                        gCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCgg                                                                                                                              >  2:448898/1‑150 (MQ=255)
                           aaCCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGc                                                                                                                           <  2:1364345/150‑1 (MQ=255)
                               acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       >  1:101079/1‑150 (MQ=255)
                               acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       >  2:7087942/1‑150 (MQ=255)
                               acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       <  1:3613032/150‑1 (MQ=255)
                                 aGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCa                                                                                                                     <  1:322262/150‑1 (MQ=255)
                                 aGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCa                                                                                                                     >  1:6813522/1‑150 (MQ=255)
                                 aGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCa                                                                                                                     <  1:577894/150‑1 (MQ=255)
                                  gCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCAc                                                                                                                    <  1:3751053/150‑1 (MQ=255)
                                          gATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCACGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGc                                                                                                            <  2:1146553/150‑1 (MQ=255)
                                           aTGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  1:6908062/150‑1 (MQ=255)
                                           aTGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  2:6940553/150‑1 (MQ=255)
                                           aTGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  1:4800820/150‑1 (MQ=255)
                                              aGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGa                                                                                                        <  2:2104510/150‑1 (MQ=255)
                                                cTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGAtt                                                                                                      >  2:1178139/1‑150 (MQ=255)
                                                acatttgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGAtt                                                                                                      <  1:7086631‑M1/64‑1 (MQ=255)
                                                acatgtgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGAtt                                                                                                      <  1:3347070‑M1/64‑1 (MQ=255)
                                                 tggtggAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTc                                                                                                     <  2:4353771/150‑1 (MQ=255)
                                                 gccgaacagcagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctccccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTc                                                                                                     <  2:2860649‑M2/60‑1 (MQ=255)
                                                  ggtggAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTcg                                                                                                    <  1:2737667/150‑1 (MQ=255)
                                                   gtggAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTcgc                                                                                                   <  2:4567229/150‑1 (MQ=255)
                                                    gaacagcagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTcgcg                                                                                                  <  1:1763133‑M2/63‑1 (MQ=255)
                                                     tgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCgg                                                                                                 >  2:3347070‑M1/82‑150 (MQ=255)
                                                      gtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCATGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                >  2:2451704‑M1/81‑150 (MQ=255)
                                                      gagaAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                <  1:7268372/150‑1 (MQ=255)
                                                       tactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGc                                                                                               >  1:4558698‑M1/80‑150 (MQ=255)
                                                       agaagaCCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGc                                                                                               >  2:5261998/1‑150 (MQ=255)
                                                        actccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCa                                                                                              >  2:5139871‑M1/79‑150 (MQ=255)
                                                         ctccaaacgaccgtcgtggcggtccttcttgcgcaagtaacagggaccacacagtccgcaagcacgccgcactccccCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCAt                                                                                             >  2:6553310‑M1/78‑150 (MQ=255)
                                                          cagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGCGTACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGgcacg                                                                                            <  2:4006629‑M2/69‑3 (MQ=255)
                                                          cagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGCGTACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGgcacg                                                                                            <  2:5794537‑M2/69‑3 (MQ=255)
                                                           agtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGc                                                                                           <  2:3320842‑M2/70‑1 (MQ=255)
                                                            caaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCgg                                                                                                 >  2:6170851‑M1/75‑143 (MQ=255)
                                                            aCCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGTTTCGCGGGCATGcc                                                                                          >  1:1136263/1‑150 (MQ=255)
                                                            aCCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGcc                                                                                          <  1:6134489/150‑1 (MQ=255)
                                                             cccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGccc                                                                                         >  2:2167385/1‑150 (MQ=255)
                                                              ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCCCGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCg                                                                                        >  1:2025599/1‑150 (MQ=255)
                                                              ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCa                                                                                              >  2:3454815/1‑144 (MQ=255)
                                                              ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCg                                                                                        >  2:1094422/1‑150 (MQ=255)
                                                              ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCg                                                                                        >  2:2587471/1‑150 (MQ=255)
                                                                 gATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTc                                                                                     <  2:4077182/150‑1 (MQ=255)
                                                                  accgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTcc                                                                                    >  2:5979918‑M1/69‑150 (MQ=255)
                                                                      ccGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCcaca                                                                                >  2:5388449/1‑150 (MQ=255)
                                                                       cGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACag                                                                               >  1:7217237/1‑150 (MQ=255)
                                                                        gTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACaga                                                                              <  1:3691836/150‑1 (MQ=255)
                                                                        gTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACaga                                                                              <  2:2679674/150‑1 (MQ=255)
                                                                        gTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGCGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACaga                                                                              <  2:5040664/150‑1 (MQ=255)
                                                                            agctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTcccacagggaga                                                                          <  2:1014409‑M2/87‑6 (MQ=255)
                                                                               tggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGggagagca                                                                       <  1:5918016‑M2/90‑3 (MQ=255)
                                                                               gCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTg                                                                       <  2:2262879/150‑1 (MQ=255)
                                                                               gCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTg                                                                       >  1:541338/1‑150 (MQ=255)
                                                                                cACGGCGCGGCTGGTCGGCAAGCTGTAAGCAGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAACCCGGTGCCATAAAACCTGTCGCCTGATTCGCGGGCATGACCGCTCCCACAGAGGGAGTgc                                                                      >  1:6540294/1‑150 (MQ=38)
                                                                                   ggCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTgcggc                                                                   <  1:4603587/150‑1 (MQ=255)
                                                                                   ggCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTgcggc                                                                   <  1:4076554/150‑1 (MQ=255)
                                                                                    tcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcgcttccCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTgcggcg                                                                  >  1:5084315‑M1/51‑150 (MQ=255)
                                                                                      gcgGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       >  1:3462305/1‑95 (MQ=255)
                                                                                            ggTCGGCCAGCTGTAAGCCGCGGACCCGGGCAACGCGCTTCCCTGGGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCg                                                          <  1:6453033/150‑1 (MQ=37)
                                                                                            ggTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCg                                                          <  1:2317074/150‑1 (MQ=255)
                                                                                                ggCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACt                                                      >  2:2946232/1‑150 (MQ=255)
                                                                                                   aaGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtg                                                   <  1:2088087/150‑1 (MQ=255)
                                                                                                   aaGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtg                                                   >  2:1655757/1‑150 (MQ=255)
                                                                                                    aGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtgt                                                  <  1:2800979/150‑1 (MQ=255)
                                                                                                      cTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTgg                                                >  1:4324525/1‑150 (MQ=255)
                                                                                                      cTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTgg                                                <  2:1280759/150‑1 (MQ=255)
                                                                                                            gCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGACCACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTg                                          <  2:3109875/150‑1 (MQ=255)
                                                                                                             ccGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCCAACACCCGCGAAGCCCGGGCCCTCCCCCGTGTCCCCAGATTCCCGGGCCATCCCCCGACCACAGAGGGCGTGCGGCGTGCTTTCGGACAGTTTTGTACCAGt                                         >  1:2640617/1‑150 (MQ=255)
                                                                                                               gCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTa                                       <  1:1959982/150‑1 (MQ=255)
                                                                                                                  gACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGCACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCAGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCAGTTACtt                                    >  1:2610884/1‑150 (MQ=255)
                                                                                                                     aaGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtgtg                                                 <  1:38517/134‑1 (MQ=255)
                                                                                                                     aaGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtgtg                                                 >  2:38517/1‑134 (MQ=255)
                                                                                                                          cAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGGGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCaaga                            >  2:881042/1‑150 (MQ=255)
                                                                                                                            aCGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCaagaag                          >  1:6026455/1‑150 (MQ=255)
                                                                                                                                gcTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGAcc                      <  2:1813280/150‑1 (MQ=255)
                                                                                                                                gcTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGAcc                      <  2:2243389/150‑1 (MQ=255)
                                                                                                                                  ttCCCTGTGGGAGCGGGCCCGCCCGCGGACACCCGCGAAGCCCGTGCCATCCACCCTGTCGCCTGATTCGCCGGCATGCCCCCACCCACAGAGGGAGTGCGGCGTGCGTTCGGCCGGTGTGGTAACAGTTACATGCGCAAGAAAGACCGc                    >  1:6290869/1‑150 (MQ=21)
                                                                                                                                  ttCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCACCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGc                    >  1:2989999/1‑150 (MQ=255)
                                                                                                                                         TGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                <  1:5605528‑M1/67‑1 (MQ=255)
                                                                                                                                             AGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCg         <  2:4773731‑M1/150‑1 (MQ=255)
                                                                                                                                              GCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGt        <  2:2749230‑M1/150‑1 (MQ=255)
                                                                                                                                               CGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGtt       <  1:2648124‑M1/150‑1 (MQ=255)
                                                                                                                                                 GGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGtcgtgtg     <  1:1826038‑M1/150‑4 (MQ=255)
                                                                                                                                                 GGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTg     >  1:14573‑M1/1‑150 (MQ=255)
                                                                                                                                                    ACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAg  >  2:1471634‑M1/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGATCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAG  >  NC_002947/2961930‑2962227

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: