Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,962,080 G→A 5.7% intergenic (+43/‑153) 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,0800GA 5.7% 934.7 / 30.2 298intergenic (+43/‑153)PP_2588/PP_2589class III aminotransferase/aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (139/142);  new base A (2/15);  total (141/157)
Fisher's exact test for biased strand distribution p-value = 2.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 301 total aligned reads displayed.

TGATCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAG  >  NC_002947/2961931‑2962227
                                                                                                                                                     |                                                                                                                                                   
   tCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGccc                                                                                                                                                  >  1:1281017/1‑150 (MQ=255)
    cATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCcg                                                                                                                                                 >  1:7445733/1‑150 (MQ=255)
      tGGCCCCGGCCTTGATCGACAACCACAGCGAGATGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGTGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCcgcg                                                                                                                                               >  2:6474794/1‑150 (MQ=255)
        gCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGaa                                                                                                                                             >  1:4260331/1‑150 (MQ=255)
        gCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGaa                                                                                                                                             <  2:5116458/150‑1 (MQ=255)
        caacgagaccgacgtcctgaccgacgctcagcgcaactatgccgaacagcagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGCACgcccgctcc                                                                                                                                             <  1:4871971‑M2/20‑4 (MQ=255)
            cGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGACGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACAcc                                                                                                                                         <  2:4783957/150‑1 (MQ=255)
             tccgaggcgcgttgttgccagaattcgagagtgtacatttgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCg                                                                                                                                        >  1:680427‑M1/121‑150 (MQ=255)
              gCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCgg                                                                                                                                       <  1:1616995/150‑1 (MQ=255)
               ccTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGc                                                                                                                                      <  2:7413303/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)
                          aaCCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGc                                                                                                                           <  2:1364345/150‑1 (MQ=255)
                            ccACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCa                                                                                                                         >  2:2737667/1‑150 (MQ=255)
                              acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       >  1:101079/1‑150 (MQ=255)
                              acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       <  1:3613032/150‑1 (MQ=255)
                              acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       <  2:4598023/150‑1 (MQ=255)
                              acaGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATc                                                                                                                       >  2:7087942/1‑150 (MQ=255)
                                aGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCa                                                                                                                     <  1:577894/150‑1 (MQ=255)
                                aGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCa                                                                                                                     >  1:6813522/1‑150 (MQ=255)
                                 gCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCAc                                                                                                                    <  1:3751053/150‑1 (MQ=255)
                                        ggATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCg                                                                                                             >  2:3501433/1‑150 (MQ=255)
                                         gATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCACGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGc                                                                                                            <  2:1146553/150‑1 (MQ=255)
                                          aTGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  1:4800820/150‑1 (MQ=255)
                                          aTGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  1:6908062/150‑1 (MQ=255)
                                          aTGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGcc                                                                                                           <  2:6940553/150‑1 (MQ=255)
                                              gCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCCCGCCCGCGAACACCGGCGAAGCCGGTGCCATCCCCCGTGTCGCCTGag                                                                                                       >  2:381073/1‑149 (MQ=255)
                                               cTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGAtt                                                                                                      >  2:1178139/1‑150 (MQ=255)
                                               cTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGAtt                                                                                                      <  2:1702499/150‑1 (MQ=255)
                                               acatttgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGAtt                                                                                                      <  1:7086631‑M1/64‑1 (MQ=255)
                                                   ttgtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTcgcg                                                                                                  >  2:2248646‑M1/83‑150 (MQ=255)
                                                   tggAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTcgcg                                                                                                  <  2:2936086/150‑1 (MQ=255)
                                                   gaacagcagtacgcgctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGCGTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTcgcg                                                                                                  <  1:1763133‑M2/63‑1 (MQ=255)
                                                    agtacaccaaaccaccgtcgtggcgcacatacttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCATGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCgg                                                                                                 >  2:3877752‑M1/82‑150 (MQ=255)
                                                     gtactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCATGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                >  2:6199821‑M1/81‑150 (MQ=255)
                                                     gagaAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                <  1:299719/150‑1 (MQ=255)
                                                     gagaAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                <  1:7268372/150‑1 (MQ=255)
                                                      tactccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGc                                                                                               >  1:4558698‑M1/80‑150 (MQ=255)
                                                        ctccaaacgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCAt                                                                                             >  2:1642596‑M1/78‑150 (MQ=255)
                                                        ctccaaacgaccgtcgtggcggtccttcttgcgcaagtaacagggaccacacagtccgcaagcacgccgcactccccCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCAt                                                                                             >  2:6553310‑M1/78‑150 (MQ=255)
                                                        aagaCCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCAt                                                                                             <  2:4985142/150‑1 (MQ=255)
                                                          gaCCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGc                                                                                           <  2:5153167/150‑1 (MQ=255)
                                                           aCCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGTTTCGCGGGCATGcc                                                                                          >  1:1136263/1‑150 (MQ=255)
                                                            cccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGccc                                                                                         >  2:4287668/1‑150 (MQ=255)
                                                             ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCCCGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCg                                                                                        >  1:2025599/1‑150 (MQ=255)
                                                             ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCTCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCg                                                                                        >  1:6161063/1‑150 (MQ=255)
                                                             ccGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCg                                                                                        >  2:1094422/1‑150 (MQ=255)
                                                               cgaccgtcgtggcggtccttcttgcgcaagtaacaggtaccacacagtccgcaagcacgccgcactccctCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCt                                                                                      >  2:256637‑M1/71‑150 (MQ=255)
                                                                 gctgttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTcc                                                                                    >  2:264635‑M2/74‑150 (MQ=255)
                                                                     ttccgcagctggtatgccagctggtcggccgagcaaaacctcgcctgagcagccactgctctccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCcaca                                                                                <  1:2653765‑M2/81‑1 (MQ=255)
                                                                     ttccgcagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCcaca                                                                                <  1:264635‑M2/81‑1 (MQ=255)
                                                                      cGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACag                                                                               >  1:7217237/1‑150 (MQ=255)
                                                                       gTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACaga                                                                              <  1:3691836/150‑1 (MQ=255)
                                                                        tAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACagag                                                                             <  2:2350189/150‑1 (MQ=255)
                                                                           cagctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTcccacagggaga                                                                          <  2:3287619‑M2/87‑6 (MQ=255)
                                                                           agctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTcccacagggaga                                                                          <  2:1014409‑M2/87‑6 (MQ=255)
                                                                             ctggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtaGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGgagagc                                                                        <  2:5042873‑M2/89‑2 (MQ=255)
                                                                             cGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGt                                                                        >  2:5934083/1‑150 (MQ=255)
                                                                              tggtatgccagctggtcggccgagcaaaacatcgcctgagcagccactgctctccctgtgGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGggagagca                                                                       <  1:5918016‑M2/90‑3 (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)
                                                                                    cgcgGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGt                                                                 <  2:3431754/150‑1 (MQ=255)
                                                                                     gcggcgGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTg                                                                <  1:3338421/150‑1 (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:439720/1‑150 (MQ=255)
                                                                                                 cAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgt                                                    >  2:2176713/1‑150 (MQ=255)
                                                                                                  aaGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtg                                                   <  1:2088087/150‑1 (MQ=255)
                                                                                                     cTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTgg                                                <  2:1280759/150‑1 (MQ=255)
                                                                                                     cTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTgg                                                >  1:4324525/1‑150 (MQ=255)
                                                                                                        tAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGc                                                           <  1:2831819/136‑1 (MQ=255)
                                                                                                        tAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGc                                                           >  2:2831819/1‑136 (MQ=255)
                                                                                                        tAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTAc                                             <  2:3344531/150‑1 (MQ=255)
                                                                                                        tAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTAc                                             <  2:1973516/150‑1 (MQ=255)
                                                                                                              gCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTa                                       <  1:1959982/150‑1 (MQ=255)
                                                                                                                 gACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGCACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCAGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCAGTTACtt                                    >  1:2610884/1‑150 (MQ=255)
                                                                                                                 gACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACtt                                    <  1:34020/150‑1 (MQ=255)
                                                                                                                    aaGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGGAGCCGGTGCCATCCACCCTGTCGCCTGATTCGCGGGCATGCCCCCGCCCACAGCGGGAGTGCGGCGGGCTTGCGGGCGGTGGGGTGCCGGGTACGTgcg                                 >  1:277334/1‑150 (MQ=37)
                                                                                                                    aaGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACtgtgtg                                                 <  1:38517/134‑1 (MQ=255)
                                                                                                                         cAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGGGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCaaga                            >  2:881042/1‑150 (MQ=255)
                                                                                                                                 ttCCCTGTGGGAGCGGGCCCGCCCGCGGACACCCGCGAAGCCCGTGCCATCCACCCTGTCGCCTGATTCGCCGGCATGCCCCCACCCACAGAGGGAGTGCGGCGTGCGTTCGGCCGGTGTGGTAACAGTTACATGCGCAAGAAAGACCGc                    >  1:6290869/1‑150 (MQ=21)
                                                                                                                                 ttCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGc                    <  2:3098015/150‑1 (MQ=255)
                                                                                                                                   cccTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCa                  <  2:6099011/150‑1 (MQ=255)
                                                                                                                                     CTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTATGGTACCTGTTACTTGCGCAAGAAGGACCGCCacg                <  2:4060570‑M1/150‑1 (MQ=255)
                                                                                                                                        TGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                >  2:5605528‑M1/1‑67 (MQ=255)
                                                                                                                                        TGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCggg                                                                                                <  1:5605528‑M1/67‑1 (MQ=255)
                                                                                                                                                GGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGtcgtgtg     <  1:1826038‑M1/150‑4 (MQ=255)
                                                                                                                                                GGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTg     >  1:6061384‑M1/1‑150 (MQ=255)
                                                                                                                                                GGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTg     >  1:14573‑M1/1‑150 (MQ=255)
                                                                                                                                                 GCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTgg    <  2:4494308‑M1/150‑1 (MQ=255)
                                                                                                                                                  CACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGa   >  2:5291698‑M1/1‑150 (MQ=255)
                                                                                                                                                  CACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGa   >  2:2161923‑M1/1‑150 (MQ=255)
                                                                                                                                                   ACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAg  >  2:1471634‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGATCATGGCCCCGGCCTTGATCGCCAACCACAGCGAGCTGGATGAGCTGGTGGAGAAGACCCGGATAGCCGTAGACCGCACGGCGCGGCTGGTCGGCAAGCTGTAAGCCGCTGACAAGGGCAACGCGCTTCCCTGTGGGAGCGGGCACGCCCGCGAACACCGGCGAAGCCGGTGCCATCCACCGTGTCGCCTGATTCGCGGGCATGCCCGCTCCCACAGAGGGAGTGCGGCGTGCTTGCGGACTGTGTGGTACCTGTTACTTGCGCAAGAAGGACCGCCACGACGGTCGTTTGGAG  >  NC_002947/2961931‑2962227

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: