Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,120 C→G 6.5% intergenic (+90/‑198) 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,1200CG 6.5% 551.1 / 3.3 184intergenic (+90/‑198)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (94/74);  new base G (4/11);  total (99/85)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

GCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTT  >  NC_002947/444971‑445269
                                                                                                                                                       |                                                                                                                                                     
gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                                       <  2:7992943/150‑1 (MQ=255)
gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                                       <  2:6678256/150‑1 (MQ=25)
 cAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCg                                                                                                                                                      <  2:3693461/150‑1 (MQ=255)
   aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCggg                                                                                                                                                    <  2:5695189/150‑1 (MQ=255)
   aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCggg                                                                                                                                                    <  2:1108682/150‑1 (MQ=255)
    aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                   <  1:7717578/150‑1 (MQ=255)
         gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                              >  2:8354663/1‑150 (MQ=255)
          tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCGAAAGCGTCCGGGCGGACAACAAATATTTGTGGTGGGGccgcgtc                                                                                                                                             >  1:6800816/1‑145 (MQ=255)
            aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTTCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                                           <  1:7924046/150‑1 (MQ=255)
              ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCTCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcgt                                                                                                                                         <  2:2562934/150‑2 (MQ=255)
              ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                         <  2:3904066/150‑1 (MQ=255)
              ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                         <  1:4546444/150‑1 (MQ=18)
              ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                         <  1:7031647/150‑1 (MQ=18)
              ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                         <  1:2370291/150‑1 (MQ=18)
              ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                         <  2:3488369/150‑1 (MQ=18)
               gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                                        <  1:5705929/150‑1 (MQ=18)
                 aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCa                                                                                                                                      <  2:7224138/150‑1 (MQ=255)
                    cgacgaAATGCTTGGCCGGATCTTCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                                   <  1:41504/150‑1 (MQ=255)
                           aTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                                     >  1:4954651/1‑125 (MQ=255)
                           aTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                                     <  2:4954651/125‑1 (MQ=255)
                                tGGCCGGACCATCGACACCTACGCCTGAATCCTGGCTTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGACCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                       <  2:7632886/150‑1 (MQ=37)
                                        acatcGACACCTACGCCTCAATACTTGCGACATCCCCTATCCCCGCGCCGGCCTGGTCGGAGCGGGCTTGCCCGCGACTGCGCCCGATCACACAACAATTTTGTCTGGTGGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTtgt                                                                                                               <  2:3768562/149‑1 (MQ=37)
                                         catcGACACCTCCGCCTTAATACTTGCGACATCCCCTATTTCAGCGCCTGCATTGTAGCTGCGGGCTTGCCCGCGAATGCGACCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTtgtg                                                                                                              <  2:8205146/150‑1 (MQ=38)
                                              acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                   <  2:286698/108‑1 (MQ=255)
                                               cacCTCCGCGTTATTCCTGGCGTCAGCCCCTTTTCCCGCGCCTGCCTTGTCGGCGGGGGCTGGCCCGCTACTGCGCCCGAACACAAAAAAATTTGGTCTGGGCGGGCCACTGCGGGGGAACGCCCGCTCCTACAAGGGTTTGtggcatgg                                                                                                        <  2:8073700/150‑1 (MQ=255)
                                                          gAATACTTGCGACAGCCCCTTTTTCCGCGCCGGCCTTGTAGGTGGGGGCTTGCCCGCGACTGCGTCCTAACAGACAACAATTTGGTCTGGTCGGGCCTGTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                                                             <  2:6640801/150‑1 (MQ=255)
                                                           atgcggggCGTCACCCCCTATGCCCGGCCCGGCGTTGTCGGCGGGGGGTGCCCCGCGACGGCGCCCTCTCACACAACACTTGTGTGTGGCGGGGCCGCGGGGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                            <  2:7900112/143‑1 (MQ=255)
                                                                        cccccAATTCCCGCGCGGGCCTTGTAGGCGGGGGCTTGCCCGCGATTGCTGCCGAACCCCAAACATTGTTGTGTGGTCGGGCCTCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                               <  2:8689215/149‑1 (MQ=255)
                                                                           ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAATTTGTCTGGTCGGGCCTCGTAGCGGGCAATCCCGCTCCTACAAGGGTTTGTTGCATGGCATGACTGATTTGCCAATTTTCGTTGCg                                                                            >  2:9155378/1‑150 (MQ=38)
                                                                           ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATTCCCGCGAAAGCCTACCATAAGACAACAATTATGTCTGGTCGGGCCGCTTCGCCGGCAAGACCCCACCGACAAGGGTTTTTGGCAAGGCATGACTGATTTTCCAATTTTCGTGGCg                                                                            >  2:7524415/1‑150 (MQ=255)
                                                                                ttCCGCGCCGGCCTTGTCGGCGGGGGCTGCCCCGCCACTGCGCCCGAACACACAAAATTTTGTTCGGGTCGGGCCTCTCCGGGGCCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt                                                                       <  2:3296616/150‑1 (MQ=255)
                                                                                         gcctttgttgGCGCGGGCGGCCCCGCGACTGCGTCCCAACCAACACCATTTTGGGGTGGTGGGGCCTCCCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                              <  1:7537879/146‑1 (MQ=255)
                                                                                                 tggAGCGGGGTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                      <  2:6740798/149‑1 (MQ=255)
                                                                                                 tggAGCGGGCTTGCGCTCTAAGGCTTCCTCTCCTCAAACATTTTGGTCGGGGGGGGCCTCTTGGCGGGAATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                      <  2:1968839/149‑1 (MQ=255)
                                                                                                 aggaggGGGCTTGCCCGCGACTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                                      <  2:2124509/150‑1 (MQ=255)
                                                                                                   gAGCGGGCTTGCCCGCTAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                    <  2:4564600/150‑1 (MQ=255)
                                                                                                    gcgggggTTGCCCGGGTATTGCGTCCTATCACACAAAATTTTTGTGGGGCGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                    <  1:2986595/147‑1 (MQ=255)
                                                                                                    aGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTACCTCCTTg                                                   >  2:5885032/1‑150 (MQ=255)
                                                                                                     gCGGGCTTGCCCGCGAAAGCGTCCGAACAGACAACAATTTTGTCTGGGCGGGCCTCTACGCGGGCAAGCCCCCTCCTACAAGGGGTTGTGGCATGGCATGACTGATTTGCCAATTTTCGTGGCGTAAATGGCACCAAAGTCCATACTTGa                                                  >  1:8819618/1‑150 (MQ=255)
                                                                                                         gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                                              >  2:694331/1‑150 (MQ=255)
                                                                                                           ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                            >  1:7281357/1‑150 (MQ=255)
                                                                                                              cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                         >  2:8784295/1‑150 (MQ=255)
                                                                                                              cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                         >  2:7230839/1‑150 (MQ=255)
                                                                                                              cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                         >  2:2994741/1‑150 (MQ=18)
                                                                                                              cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                         >  1:5512363/1‑150 (MQ=18)
                                                                                                               ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                       <  1:7185886/71‑1 (MQ=255)
                                                                                                               ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                                       >  2:7185886/1‑71 (MQ=255)
                                                                                                               ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGTCGTAAAt                                                                       >  2:7172161/1‑119 (MQ=255)
                                                                                                               ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                        >  2:6175340/1‑150 (MQ=255)
                                                                                                               ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                        >  1:6503236/1‑150 (MQ=255)
                                                                                                               ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                        >  2:3094776/1‑150 (MQ=255)
                                                                                                                cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                       >  2:1368129/1‑150 (MQ=255)
                                                                                                                   gAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCt                                    >  2:1437650/1‑150 (MQ=255)
                                                                                                                      agaggcccGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                                  >  1:9042438/7‑150 (MQ=21)
                                                                                                                          gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                                  >  2:1601660/1‑146 (MQ=255)
                                                                                                                          gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                              >  2:400243/1‑150 (MQ=255)
                                                                                                                          gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                              >  1:1923021/1‑150 (MQ=255)
                                                                                                                           tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                             >  1:7493975/1‑150 (MQ=255)
                                                                                                                           tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                             >  2:9122252/1‑150 (MQ=255)
                                                                                                                            ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                       <  1:1263837/139‑1 (MQ=255)
                                                                                                                            ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                            >  2:7055818/1‑86 (MQ=18)
                                                                                                                            ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                            <  1:7055818/86‑1 (MQ=18)
                                                                                                                             cGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                           >  1:7225067/1‑150 (MQ=38)
                                                                                                                               aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                         >  1:5121193/1‑150 (MQ=255)
                                                                                                                               aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                         >  1:6246004/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                     >  2:979310/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                     >  1:1677932/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                     >  1:6678504/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                     >  2:6478689/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                                          <  2:4220457/62‑1 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAtaca                  >  1:3635340/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  1:5587028/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  2:7732837/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  2:511333/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAAGTAAGTTTCAaaca                  >  1:6986146/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTAGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  2:7671976/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGAAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  2:7672240/1‑150 (MQ=255)
                                                                                                                                      aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  1:1189921/8‑150 (MQ=37)
                                                                                                                                      aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                  >  1:4642264/8‑150 (MQ=37)
                                                                                                                                       acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGTCGGAATGCTGAGCTAAAGGTAAGTTTCaatctc                 >  1:8122393/1‑146 (MQ=255)
                                                                                                                                       aaatttttGTAGGGTCGGGCCTCCCCGGGGCCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGATAGCTGAGCTAAACGTAAGTTTCAaacaa                 <  1:3649191/146‑1 (MQ=255)
                                                                                                                                        aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac                >  2:8556896/6‑150 (MQ=38)
                                                                                                                                         attgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt               >  1:2539654/5‑150 (MQ=38)
                                                                                                                                         aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt               >  1:8377065/1‑150 (MQ=255)
                                                                                                                                         aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt               >  2:8170980/1‑150 (MQ=255)
                                                                                                                                         aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGATGAGCTAAACGTAAGTTTCAAACAACt               >  1:7083362/1‑150 (MQ=255)
                                                                                                                                           ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt             >  2:20738/1‑150 (MQ=255)
                                                                                                                                           ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACGGt             >  2:7853607/1‑150 (MQ=255)
                                                                                                                                            tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt            >  1:1407600/1‑150 (MQ=255)
                                                                                                                                            tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt            <  1:979310/150‑1 (MQ=255)
                                                                                                                                            ttggtcgggtcggGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt            <  2:5074692/147‑1 (MQ=255)
                                                                                                                                                  cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc       >  2:345749/1‑150 (MQ=255)
                                                                                                                                                  cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc       >  1:8031584/1‑150 (MQ=255)
                                                                                                                                                  cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTATCAAACAGTTGTTTGaaa       <  2:6706467/150‑3 (MQ=255)
                                                                                                                                                   tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg      >  1:4381200/1‑150 (MQ=255)
                                                                                                                                                     gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc    >  1:3131552/1‑150 (MQ=255)
                                                                                                                                                     gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAATGTTTGATCGGc    >  1:8459563/1‑150 (MQ=255)
                                                                                                                                                      tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc    >  2:4840708/1‑149 (MQ=255)
                                                                                                                                                      tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc    <  1:4840708/149‑1 (MQ=255)
                                                                                                                                                      tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCt   >  1:1291141/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTT  >  NC_002947/444971‑445269

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: