Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,065 C→G 5.3% intergenic (+29/‑111) PP_3198 → / → PP_3199 2Fe‑2S ferredoxin/monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,0650CG 5.3% 764.0 / 23.0 227intergenic (+29/‑111)PP_3198/PP_31992Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (117/98);  new base G (4/8);  total (121/106)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 235 total aligned reads displayed.

TTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTGGGCGGCGG  >  NC_002947/3629917‑3630213
                                                                                                                                                    |                                                                                                                                                     
ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                                                      >  1:3670115/1‑150 (MQ=255)
ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                                                      >  1:3265094/1‑150 (MQ=255)
ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                                                      >  2:3356469/1‑150 (MQ=255)
ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                                                      >  2:3698277/1‑150 (MQ=255)
 tCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGATGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcg                                                                                                                                                     >  2:5397509/1‑150 (MQ=255)
     gCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgac                                                                                                                                                 >  1:2250244/1‑150 (MQ=255)
      cccGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCAATCGCGCCGCGCGATAGGTCGCATCCcgcgacg                                                                                                                                                <  1:1856404/150‑1 (MQ=255)
        cGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgacgcg                                                                                                                                              >  1:4747777/1‑150 (MQ=255)
          caGACGGCACGTGCGTTCCCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAt                                                                                                                                            <  2:3922643/149‑1 (MQ=255)
           aGACGGCACTTGCGTTCATGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                                                           <  2:1084560/150‑1 (MQ=255)
           aGACGGCACTTGCGTTCATGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                                                           <  2:2865100/150‑1 (MQ=255)
           aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                                                           <  1:1381404/150‑1 (MQ=255)
           aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCAATCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGAta                                                                                                                                           <  1:3721003/150‑2 (MQ=37)
           aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCAATCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGAta                                                                                                                                           <  2:450634/150‑2 (MQ=37)
            gACGGCACTTGCGTTCCCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                                          <  2:2471489/150‑1 (MQ=255)
            gACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                                          <  2:3158368/150‑1 (MQ=255)
            gACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCAATCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGAtag                                                                                                                                          <  1:2129440/150‑3 (MQ=37)
                      gCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGAcc                                                                                                                                >  2:4373706/1‑150 (MQ=255)
                       cGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCg                                                                                                                               <  1:4785828/150‑1 (MQ=255)
                        gTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCgg                                                                                                                              >  1:4607746/1‑150 (MQ=255)
                         ttCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGgc                                                                                                                             <  1:2481329/150‑1 (MQ=255)
                          tCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGgcg                                                                                                                            >  1:4248796/1‑150 (MQ=255)
                          tCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCGATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGgcg                                                                                                                            >  1:4903965/1‑150 (MQ=255)
                            aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGgcgcc                                                                                                                          >  2:4123493/1‑149 (MQ=255)
                            aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGACGCTAGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCa                                                                                                                          >  1:345835/1‑150 (MQ=255)
                               gCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGc                                                                                                                       <  1:3184328/150‑1 (MQ=255)
                                ccccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTTCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                                                      <  2:1000035/150‑1 (MQ=255)
                                ccccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                                                      <  1:4919029/150‑1 (MQ=255)
                                    tGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                                                  >  1:3180351/1‑150 (MQ=255)
                                      ccGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                                          >  2:1525328/1‑124 (MQ=255)
                                      ccGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGc                                                                                                                >  2:1881053/1‑150 (MQ=255)
                                       cGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCt                                                                                                               <  1:663908/150‑1 (MQ=255)
                                        gAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtg                                                                                                              >  2:4813610/1‑150 (MQ=255)
                                         aGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtgt                                                                                                             <  1:4760269/150‑1 (MQ=255)
                                             gAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtgtg                                                                                                            <  2:1763016/147‑1 (MQ=255)
                                             gAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAt                                                                                                         >  1:5231274/1‑150 (MQ=255)
                                              aaGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATa                                                                                                        <  1:811878/150‑1 (MQ=255)
                                              aaGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATa                                                                                                        <  2:4697983/150‑1 (MQ=255)
                                               aGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTAGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCCTat                                                                                                       <  2:345835/150‑2 (MQ=255)
                                                gCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAgg                                                                                                      <  1:997637/150‑1 (MQ=255)
                                                 ccTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGc                                                                                                     >  2:1258937/1‑150 (MQ=255)
                                                   ttGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGcac                                                                                                   >  1:3116135/1‑150 (MQ=255)
                                                           aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCAt                                                                                                                                     >  1:1190007/1‑108 (MQ=255)
                                                           aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACg                                                                                           >  1:1770697/1‑150 (MQ=255)
                                                           aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACg                                                                                           >  1:5107440/1‑150 (MQ=255)
                                                           aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACg                                                                                           >  2:3227015/1‑150 (MQ=255)
                                                              cGATCACAAGGACGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgcc                                                                                        <  1:4721824/150‑1 (MQ=255)
                                                                aTCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgccgc                                                                                      >  2:5171195/1‑150 (MQ=255)
                                                                 tCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgccgcc                                                                                     >  2:1744198/1‑150 (MQ=255)
                                                                     aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGc                                                                                >  1:578077/1‑150 (MQ=255)
                                                                     aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGc                                                                                <  2:4720275/150‑1 (MQ=255)
                                                                       ggAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCAt                                                                              >  1:1706377/1‑150 (MQ=255)
                                                                              ggcggcTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACgggg                                                                       <  2:859957/150‑1 (MQ=255)
                                                                                    tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGcaca                                                                 >  1:2547756/1‑150 (MQ=255)
                                                                                      gctgctATTGCCGCTTGAGCAGCTCGATCGCGGACCTAGCGCGCCGCGCGATAGGTCGCACCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGTCATGCGGGgcggccgtt                                                               <  2:218464/150‑7 (MQ=255)
                                                                                        tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCAcaacaa                                                             <  1:2211972/150‑1 (MQ=255)
                                                                                        tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCAcaacaa                                                             >  2:2322819/1‑150 (MQ=255)
                                                                                          ctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCCTCGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                                           >  2:572681/1‑150 (MQ=255)
                                                                                             ttGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACa                                                        >  1:2749985/1‑150 (MQ=255)
                                                                                             ttGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACa                                                        >  1:2086896/1‑150 (MQ=255)
                                                                                               gCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACACCAAGAACaaa                                                      >  2:3644923/1‑150 (MQ=255)
                                                                                                ccGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTc                                                                                                                                       <  2:1479202/69‑1 (MQ=255)
                                                                                                ccGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTc                                                                                                                                       >  1:1479202/1‑69 (MQ=255)
                                                                                                     tGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGTTTGGTCATGGACGGGCGCTGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                                <  1:1053348/150‑1 (MQ=255)
                                                                                                     tGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGAATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                                <  1:1915594/150‑1 (MQ=255)
                                                                                                        gcagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGAc                                             >  2:535654/1‑150 (MQ=255)
                                                                                                         cagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACa                                            <  1:3034126/150‑1 (MQ=255)
                                                                                                          agcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACAt                                           >  1:4010398/1‑150 (MQ=255)
                                                                                                            cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGc                                         <  1:3756361/150‑1 (MQ=255)
                                                                                                            cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGc                                         >  2:3441479/1‑150 (MQ=255)
                                                                                                             tCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGcc                                        >  2:4613108/1‑150 (MQ=255)
                                                                                                             tCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGcc                                        >  1:93423/1‑150 (MQ=255)
                                                                                                             tCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGcc                                        >  2:5413104/1‑150 (MQ=255)
                                                                                                               gATCGCTGACCTATCGCGCCGCGCGAGAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAGAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGAcgcctctg                                      >  2:4980356/1‑143 (MQ=255)
                                                                                                                ggatgcgACCTATCGCGCGGCGCGAGAGGTCGCAACCCGCGCCGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATg                                     >  1:2391411/7‑150 (MQ=255)
                                                                                                                aTCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATg                                     >  1:2029366/1‑150 (MQ=255)
                                                                                                                     tGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATAc                                >  2:3294460/1‑150 (MQ=255)
                                                                                                                     tGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGGGAGACATGCCATGAATAc                                >  1:1010289/1‑150 (MQ=255)
                                                                                                                        ccTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt                             >  2:1158585/1‑150 (MQ=255)
                                                                                                                        ccTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt                             <  2:2250244/150‑1 (MQ=255)
                                                                                                                         cTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg                            >  2:2463260/1‑150 (MQ=255)
                                                                                                                         cTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACACCAAGAACAAAACGTGAGACATGCCATGAATACtgtg                            >  2:1184860/1‑150 (MQ=255)
                                                                                                                          tATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGa                           <  1:4373706/150‑1 (MQ=255)
                                                                                                                          tATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGa                           >  2:5299262/1‑150 (MQ=255)
                                                                                                                          tATCGCGACGCGCGATAGGTCGAATCCCGCGACGCGATTGGTCAGGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACCTGCCATGAATACTGTGa                           >  2:4909748/1‑150 (MQ=255)
                                                                                                                          aaTCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                          >  1:3777937/2‑119 (MQ=25)
                                                                                                                          aaTCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                          <  2:3777937/118‑1 (MQ=25)
                                                                                                                           aTCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaa                          >  1:2916601/1‑150 (MQ=255)
                                                                                                                           aTCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaa                          >  1:5074834/1‑150 (MQ=255)
                                                                                                                           aTCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaa                          >  1:4981864/1‑150 (MQ=255)
                                                                                                                            tCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaac                         <  2:3997398/150‑1 (MQ=255)
                                                                                                                            tCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaac                         >  2:769839/1‑150 (MQ=255)
                                                                                                                               cgccgcgcGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaacaac                      <  2:2797654/150‑1 (MQ=255)
                                                                                                                                 ccgcgcGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACgt                    >  1:3065838/1‑150 (MQ=255)
                                                                                                                                   gcgcgATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACgtgt                  <  1:1258937/150‑1 (MQ=255)
                                                                                                                                        aTAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATc             <  2:3670115/150‑1 (MQ=255)
                                                                                                                                            gTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTgg         >  2:5076395/1‑150 (MQ=255)
                                                                                                                                                 aTCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGACTACTGTGAACAACGTGTTGATCGTGggcggc    >  1:3543120/1‑150 (MQ=255)
                                                                                                                                                 aTCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg                            <  2:4909890/126‑1 (MQ=255)
                                                                                                                                                   cccGCGACGCGATTGGTCATGGACCGGCGCTGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTGGGcggcgg  >  2:5152699/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTGGGCGGCGG  >  NC_002947/3629917‑3630213

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: