Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,110 T→G 10.9% intergenic (+80/‑208) PP_0367 → / → PP_0368 hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947445,1100TG10.9% 266.2 / 20.5 97intergenic (+80/‑208)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (48/38);  new base G (8/3);  total (56/41)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 107 total aligned reads displayed.

CGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGT  >  NC_002947/444961‑445258
                                                                                                                                                       |                                                                                                                                                    
cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                                      <  2:152969/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  1:779342/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  1:4272677/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  1:4272587/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  2:2591090/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                    <  1:2871040/150‑1 (MQ=255)
    cTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTTAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGt                                                                                                                                                  <  1:4239790/150‑1 (MQ=255)
     tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCGCGCATGCGGGCCCTCAGCCAACACAtttttg                                                                                                                                                  >  2:3966465/1‑148 (MQ=255)
     tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                                 <  2:858133/150‑1 (MQ=255)
     tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                                 <  1:3665901/150‑1 (MQ=255)
        gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCGCAAGCGTCCCATCACACAACCAATTTGtgttg                                                                                                                                              >  2:3149771/1‑146 (MQ=255)
        gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                              <  2:1216446/150‑1 (MQ=255)
        gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                              <  2:1216374/150‑1 (MQ=255)
           cAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAATCGTCCGATCAGACAACAAATTTTTGTGGtgg                                                                                                                                           >  2:1848665/1‑148 (MQ=255)
            aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                          <  1:1119070/150‑1 (MQ=255)
            aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                          <  2:3269220/150‑1 (MQ=255)
            aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGACTTTCCCGCGAATGCGTCCGATAAGACAACAATTTTGTCTGGTCgg                                                                                                                                          >  2:1334246/1‑150 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGAAAAAAATTTTGTCTGGGAGGGc                                                                                                                                        >  2:2167607/1‑150 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                        <  1:1285967/150‑1 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                        <  1:1286212/150‑1 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                        <  1:1286283/150‑1 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCCAAAGCGTCCGATCAGAAAAAAAATTTGTGTGGTGGGGc                                                                                                                                        >  2:2167661/1‑150 (MQ=255)
                  cgtcAAGGCATCCGACGAAATGCTTGGCCTGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCggacgcat                                                                                                                                    <  1:1362749/150‑7 (MQ=21)
                      aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                                <  1:2671737/150‑1 (MQ=255)
                      aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                                <  1:415889/150‑1 (MQ=255)
                       aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATGTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                               <  2:3833341/150‑1 (MQ=255)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACGTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTGTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                              <  2:4068803/150‑1 (MQ=255)
                         gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                             <  2:406792/150‑1 (MQ=18)
                         gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACGTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCggg                                                                                                                             <  2:1119492/150‑1 (MQ=255)
                          cATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                                            <  1:235520/150‑1 (MQ=255)
                            tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTATGGTCGGGCCACTTCGCCGGCaa                                                                                                                          >  2:2967771/1‑150 (MQ=255)
                              cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                        <  1:1532747/150‑1 (MQ=255)
                              cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAAGCGTCCGATCAGACAACCATTTTGTATGGTGGGGCCTCGTCGCGGGCAAGc                                                                                                                        >  2:2628377/1‑150 (MQ=255)
                               gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                                       <  1:1438010/150‑1 (MQ=255)
                               gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                                       <  1:1607538/150‑1 (MQ=255)
                               gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                                       <  2:837072/150‑1 (MQ=255)
                                          gggCCGCATCATCGACACCTACGCCTTAATACTTGCGACATCCCCTTTTTCAGCGCCTGCGTTGGCGGAGGGGGCTTGCCGGCTAATGCGTCCTATCTTAAAACATTTTTGTGTGGGGGGGCCTCTCCGGGGGCAAGCCCGCTCCTACaa                                                                                                            <  2:1463008/149‑1 (MQ=255)
                                             ccGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGATGCCCCCCAATTCCTACGATAAGACAACAATTTTGTCAGGTCGGGCCTCTTCCAGGGCCATCCCCCGCCGACAAcgg                                                                                                         >  2:1019611/1‑150 (MQ=255)
                                                     tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGGAAGCCGCCGGGCAGGCAACCATTTTGTGGGGGCGGGGCGCATCGCGGGCCAGCCCCCGCCGACAACGGGTTGtggg                                                                                                 >  2:2531054/1‑149 (MQ=255)
                                                     tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCCCATGCGGCCCCTCAGACAACCATTTTGTCGGGGGGGGCCGCATGGCGGGCCACCCCGCGCCGACAAGGGGTTGtggg                                                                                                 >  2:2531046/1‑149 (MQ=255)
                                                        acacCTCCGCCTCAATACTGGCGACCCCCCCTTTTGCCGCGCCTCCCTTGTGGGCGGGGGCTGGCCCGGTATGGGGCCCTATCACACACCATTTTTGTGGGGGGGGGCCGCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGtggcatg                                                                                              <  2:3895915/150‑1 (MQ=255)
                                                          cccTACGCCTGACTCCGGGCGTCACCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGACTGGGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                            <  2:2568861/149‑1 (MQ=25)
                                                          acCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCTAATGCGCCCGATCAGACAACATTTTTGTCTGGTCGGGCCTCTTGGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                            <  2:2568753/150‑1 (MQ=255)
                                                           cctcccccTCAATCCTGGCGACACCCCCTATTCCAGCCCCGGCCTTGTAGGAGGGGGCTGGCCCGCCAAGGCGCCCGAACAGACACCAATTTTGGGGGGCGGGGCCTCCCCGGGGCAAACCCCGCTCCTACAAGGGTTTGtggcatggca                                                                                           <  2:3440447/144‑1 (MQ=255)
                                                                  cTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCCAAGGCGCCGCTCAGAAAAAAAATTTTTATGGTGGGGGCTCGTAGCCGGCCACCCCCCGACCAAAACGGGGTGTGGGATGGGATGaaata                                                                                    >  2:1747790/1‑146 (MQ=255)
                                                                               tCACCCCCTATTTCCGCGCCGGCCTGGTCGGCGGGGGCTGGCCCGCGAATGCGTCCGATCATCAAAAAATTTGGCCGGGTCGGGCCTCTTCGGGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTc                                                                       <  2:3273040/150‑1 (MQ=255)
                                                                                cATCCCCTATTTCAGCGCCTGCCTTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGCCAACAATTTGGTCTGGGCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                                      <  2:4107782/150‑1 (MQ=39)
                                                                                 ccccccTCTTCCCGCCCCGGCGTTGTCGGCGCGGGCTTGCCCGCGAATGGGCCCGAACACACAACAATTTTGGGGGGGCGGGCCCCTGCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                                     <  2:3974916/148‑1 (MQ=255)
                                                                                 aTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCATACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                                     <  2:4013616/150‑1 (MQ=255)
                                                                                  tCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTGCCGCCGCAGCGTACCCTACTACAACAAATTTTTGGGGTGGGGCGGCTTAGCGGGCAACCCCGCGCCTACCAGGGTTTGTTGGATGGGCTGACAGGTTTTACAattttattt                                                                    >  2:4401789/1‑146 (MQ=255)
                                                                                      cTATTCCCGCGCCTGCCTTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGt                                                                <  2:710547/150‑1 (MQ=255)
                                                                                                ccctgcctTGTTGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTGTGTCTGGTCGGCCCTATTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCAcc                                                      <  2:4428267/149‑1 (MQ=38)
                                                                                                           attcgCGGGCAAGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                           >  2:4095695/5‑150 (MQ=255)
                                                                                                                      ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                 >  1:4325168/1‑150 (MQ=255)
                                                                                                                      ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                 >  2:3317202/1‑150 (MQ=255)
                                                                                                                      ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                 >  1:893715/1‑150 (MQ=17)
                                                                                                                       tGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAGTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAg                                >  1:3261976/1‑150 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGGCTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                             >  1:1762061/1‑150 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                             >  2:1427949/1‑150 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                             >  1:1996289/1‑150 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                             >  1:4502842/1‑150 (MQ=255)
                                                                                                                          acgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                             >  1:3543912/2‑150 (MQ=255)
                                                                                                                           cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                                      >  2:3888550/1‑108 (MQ=255)
                                                                                                                           cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                                      <  1:3888550/108‑1 (MQ=255)
                                                                                                                            gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                           >  2:3781449/1‑150 (MQ=255)
                                                                                                                             ccacggcgcCCGACCACCCAACATGTTGGTCGGGCCGGCCCCCGCCGCGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGc                          <  1:2842418/145‑1 (MQ=255)
                                                                                                                                aTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                       >  2:3854862/1‑150 (MQ=255)
                                                                                                                                    gcccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                   >  1:1581323/3‑150 (MQ=38)
                                                                                                                                    gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                   >  1:1747552/1‑150 (MQ=255)
                                                                                                                                     tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                  >  1:2979977/1‑150 (MQ=255)
                                                                                                                                     tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                  >  2:3311984/1‑150 (MQ=255)
                                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                 >  1:2766510/1‑150 (MQ=38)
                                                                                                                                      ccGACCAGACAAAATTGTTGACTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                 >  2:3440705/1‑150 (MQ=38)
                                                                                                                                        gaccAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt               >  1:4199226/4‑150 (MQ=38)
                                                                                                                                        gaccAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt               >  1:4199568/4‑150 (MQ=38)
                                                                                                                                        gATCAGACAACAATGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAAAGTCCATCCGTGACAGAAAGAAGAGATAAACGTAAGt               >  1:1694232/1‑150 (MQ=255)
                                                                                                                                         aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt              >  2:2139525/1‑150 (MQ=255)
                                                                                                                                         aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt              >  2:2139359/1‑150 (MQ=255)
                                                                                                                                         aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt              >  2:2230677/1‑150 (MQ=255)
                                                                                                                                         aTCAGACAACAATTTTGTCTGCACGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTTTGGCATGGCATGACAGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAgat              >  1:2794525/1‑148 (MQ=255)
                                                                                                                                           cAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc            >  1:1179806/1‑150 (MQ=255)
                                                                                                                                            ccacaacaATTTTGTCTGGTGGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa           <  2:3652350/148‑1 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                        >  1:2844663/1‑120 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                        <  2:2844677/120‑1 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                        <  2:2844663/120‑1 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                        >  1:2844677/1‑120 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAATTTTCAaaca       >  2:2335754/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:1602754/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  1:1319069/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:455824/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  1:3284978/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  1:703413/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:1330467/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:2457980/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:2841810/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:3007582/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTGTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca       >  2:3141458/1‑150 (MQ=255)
                                                                                                                                                  caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac     >  1:664904/1‑150 (MQ=255)
                                                                                                                                                   aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt    >  2:2658562/1‑150 (MQ=255)
                                                                                                                                                      ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt  <  2:4076612/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGT  >  NC_002947/444961‑445258

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: