Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,118 G→A 16.1% intergenic (+88/‑200) 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,1180GA16.1% 497.7 / 76.5 202intergenic (+88/‑200)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (109/58);  new base A (12/20);  total (122/80)
Fisher's exact test for biased strand distribution p-value = 5.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 210 total aligned reads displayed.

CGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGC  >  NC_002947/444970‑445267
                                                                                                                                                    |                                                                                                                                                     
cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGt                                                                                                                                                      <  1:3032355/150‑1 (MQ=255)
 gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                                     <  2:4172097/150‑1 (MQ=255)
 gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                                                     <  2:7421650/150‑1 (MQ=255)
 gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                                     <  2:2511836/150‑1 (MQ=38)
  cAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCg                                                                                                                                                    <  2:4799449/150‑1 (MQ=38)
     aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                 <  1:1835506/150‑1 (MQ=255)
        tcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                                              <  1:3284679/150‑1 (MQ=255)
              aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                                        <  2:5462658/150‑1 (MQ=255)
               ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  2:6329655/150‑1 (MQ=255)
               ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  1:3210419/150‑1 (MQ=255)
               ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  1:1091651/150‑1 (MQ=255)
               ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  2:6929687/150‑1 (MQ=18)
               ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  1:1151683/150‑1 (MQ=18)
               ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTGGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  2:3829566/150‑1 (MQ=18)
                gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                                      <  1:6549812/150‑1 (MQ=18)
                gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                                      <  2:3375278/150‑1 (MQ=18)
                gCATCCGACGAAATGCTGGGCCGGATCATCGACACCTACGCCTGAATACGTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACACTTTTGTCTGGTCGGGCCTCTTCGCggg                                                                                                                                      <  2:4459994/150‑1 (MQ=255)
                   tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATGTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCaa                                                                                                                                   <  2:4590953/150‑1 (MQ=17)
                   tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACGTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCaa                                                                                                                                   <  1:1450621/150‑1 (MQ=17)
                     cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGCTCCCACAACAATGTTGTGGGGTCGGGCCGCGTCGCGGGCAAGc                                                                                                                                 >  2:6108025/1‑150 (MQ=255)
                     cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTAGTGCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                                                 <  2:5508375/150‑1 (MQ=21)
                           aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATGCCCGCGAAAGCCTCCCAACAGACAACAAATTTTTGTGGTCGGGCCGCGTCGCGGGCAAGCCCGcac                                                                                                                           >  2:3712836/1‑148 (MQ=255)
                                ttGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                 >  1:635124/1‑123 (MQ=255)
                                ttGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                 <  2:635124/123‑1 (MQ=255)
                                   gCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCCCAGGCGGCCGGACAGACAAAAAAATTTTCTGGGCGGGCCGCGTAGCCGGCCAACCCCCACcaacaaag                                                                                                                   >  2:2243333/1‑148 (MQ=255)
                                     cGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       >  1:4764120/1‑128 (MQ=14)
                                                    tACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTGCCCGCGAAAGCGTCCGGTCAGAGAACAATTTTTTGTGGTCGGGGCTCGTGGCGGGCAAGCCCCCTCCTAAAAGGGTTTGTGggatggcatg                                                                                                  >  2:1082000/1‑150 (MQ=255)
                                                     ccGCCTGAATACTTGCGTCATCCCCATTTGCCGCGCCTGCCTTGTAGGTGGGGGCTGGCCGGCTAATGCGTCGCATCACAAAACATTTTTGTCTGGGGGGGCCTCTCCGCGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGa                                                                                                 <  2:6474100/149‑1 (MQ=255)
                                                         ctcaaTACTTGCTTCAACCCCAATTTCACCCCCTGCCTTGTAGCACCGGGCTTGCCCGCGAATGCGACCTATAACAAAACAATTTTGTCTGGTCGGCCATCTCCGCGGCCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGa                                                                                             <  2:3050527/147‑1 (MQ=255)
                                                                      tCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGCATGCGTCCGGTCAGACAACAAATTTGTCTGGTCGGGCCACTTCGCGGGCAAGCCCCCGCCTACAAGGGTTTGTGGCATGGCAGGGCGGATTTTCCaaatttg                                                                                >  2:6146892/1‑144 (MQ=255)
                                                                           cccTATGCCCGCCCCTGCCTGGTCGGTGCGGGCTCCCCCGCGACGGCGTCCCACCCCACAACAATTTGGGGGGGGGGGCCCTCTCCGGGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                                           <  2:5839519/150‑1 (MQ=255)
                                                                                 ttCCGCCCCTGCATTTTAGGAGCGGGCTTCCCCGCGAATCCCTCCCAACACAAAAAAATTTTGTCTGCTCGGCCCACTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt                                                                     <  2:4011687/150‑1 (MQ=255)
                                                                                     gcgcCTGCCTTGTAGGAGGGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                                 <  2:3589263/150‑1 (MQ=255)
                                                                                          tgcctTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                            <  2:6045311/150‑1 (MQ=255)
                                                                                               tGTAGGAGCGGGCTCGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTGTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCa                                                       >  2:6809694/1‑150 (MQ=255)
                                                                                                    gtgggggCTTGCCCGCGACGGCGGCCGATCTGACAACATTTTTGTCGGGGCGGGCCTGTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                  <  1:5560259/146‑1 (MQ=255)
                                                                                                    gcgggggCTGGCCCGCTACTGCGTCCTAACACACAACATTTTGGTGTGGTCGGGCCGCCCCGCGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                  <  1:7385648/146‑1 (MQ=255)
                                                                                                        ggggTTGCCCGCGAATGCGTCCTACCACAAAACAATTTTGTCTGGTCGGGCCTCTTCGCGGTCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa                                              <  2:5631718/150‑1 (MQ=255)
                                                                                                        gggCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa                                              >  2:3319172/1‑150 (MQ=255)
                                                                                                        ggacTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa                                              >  2:7117287/4‑150 (MQ=255)
                                                                                                            ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                          >  2:1288618/1‑150 (MQ=255)
                                                                                                            ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                          >  1:893790/1‑150 (MQ=255)
                                                                                                            ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                          >  2:1316156/1‑150 (MQ=255)
                                                                                                               cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt                                                                     >  2:1409642/1‑120 (MQ=255)
                                                                                                                cggctaaTGCGTCCGATCAGACAACAATTTGGTCTGGCCGGCCCTCTCCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      <  2:4838312/145‑1 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:7060904/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:2393592/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  1:1055666/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:2583834/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  1:258968/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  1:3331389/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  1:292077/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  1:5627048/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:2765453/1‑150 (MQ=255)
                                                                                                                ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCATTGACAGAAAGCTg                                      >  2:5312306/1‑150 (MQ=255)
                                                                                                                 cgcgAATGCGTCCGATCAGACAACCATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     >  2:6181411/1‑150 (MQ=255)
                                                                                                                 cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     >  1:353836/1‑150 (MQ=255)
                                                                                                                 cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCACGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     <  1:36385/150‑1 (MQ=255)
                                                                                                                 cgcgAATGCGTCCGAACAGACAACAATTTGGTCTGGTCGGGCCTCTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     <  2:4964016/150‑1 (MQ=255)
                                                                                                                     aagaggcccGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTa                                 >  2:6702489/8‑150 (MQ=21)
                                                                                                                     aaTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTa                                 >  2:905094/1‑150 (MQ=255)
                                                                                                                       gaggcccGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGGCAGAAAGCTGAGCTaaa                               >  1:7361797/6‑150 (MQ=21)
                                                                                                                           tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACAGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                           >  1:707572/1‑150 (MQ=255)
                                                                                                                            ccGATCAGACAACAATTTTGTCTGGTCGGGACTCTTCGCGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTTATCTAAACGTaa                          >  1:1587928/1‑150 (MQ=255)
                                                                                                                            ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTCGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                          >  1:3746877/1‑150 (MQ=38)
                                                                                                                               aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                       >  2:6561130/1‑150 (MQ=255)
                                                                                                                               aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                       >  1:7173802/1‑150 (MQ=255)
                                                                                                                                  aGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa                    >  1:6977462/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                   >  1:1404579/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                   >  1:3145237/1‑150 (MQ=255)
                                                                                                                                   gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                   >  2:1541640/1‑150 (MQ=255)
                                                                                                                                    acaacaATTTTGTGGGGTCGGGCCTGTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                  <  2:5199545/150‑1 (MQ=255)
                                                                                                                                    acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                  >  1:1720789/1‑150 (MQ=255)
                                                                                                                                    acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                  >  1:4608573/1‑150 (MQ=255)
                                                                                                                                     caacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaac                 >  1:397116/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                  <  2:4870054/148‑1 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                  >  1:4870054/1‑148 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:282965/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:1167857/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:4474913/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:761136/1‑150 (MQ=255)
                                                                                                                                      aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAACGGTTTGTGGCAAGGCATGACTGATTTGCCAATTTTATTGGCGTAAAAGGCACCATAGTCCATCCGTGACAGAAAGCAGAGCTAAACGTAAGTTTCAaaca                >  1:6237276/1‑150 (MQ=255)
                                                                                                                                      aaaattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  2:1876962/8‑150 (MQ=37)
                                                                                                                                       acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa               >  1:6297545/1‑150 (MQ=255)
                                                                                                                                       aaatttttGTCTGTTGGGGCCTCTTCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa               <  2:1843935/146‑1 (MQ=255)
                                                                                                                                        caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              >  2:6976119/1‑150 (MQ=255)
                                                                                                                                        aatttgtgtgtgGTGGGGCCTGGTGGGGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              <  1:4114113/144‑1 (MQ=255)
                                                                                                                                        aattgttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              >  2:2616389/6‑150 (MQ=255)
                                                                                                                                        aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              >  2:1968833/6‑150 (MQ=38)
                                                                                                                                        aattgttGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                          <  2:4751095/133‑1 (MQ=255)
                                                                                                                                         aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt             >  1:6898962/1‑150 (MQ=255)
                                                                                                                                         aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt             >  1:226099/1‑150 (MQ=255)
                                                                                                                                           ttttGTGGGGTGGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt           <  2:78387/150‑1 (MQ=255)
                                                                                                                                           ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCGATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt           >  2:5088642/1‑150 (MQ=255)
                                                                                                                                            tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGCAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt          >  1:5230420/1‑150 (MQ=255)
                                                                                                                                            tttGTCGGGGGGGGCCTCTGCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt          <  1:2434898/150‑1 (MQ=255)
                                                                                                                                                 cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc     <  2:7173802/150‑1 (MQ=255)
                                                                                                                                                   ggTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTc                                                         >  1:944789/1‑96 (MQ=255)
                                                                                                                                                    gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc  >  1:2346531/1‑150 (MQ=255)
                                                                                                                                                    gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc  >  1:4677201/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGC  >  NC_002947/444970‑445267

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: