Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,761 G→A 5.0% I226I (ATC→ATT PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,7610GA 5.0% 570.1 / 12.3 179I226I (ATC→ATTPP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (95/75);  new base A (7/2);  total (102/77)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 180 total aligned reads displayed.

CCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCG  >  NC_002947/3975613‑3975910
                                                                                                                                                    |                                                                                                                                                     
   gCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATAc                                                                                                                                                   >  1:938817/1‑150 (MQ=255)
   gCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATAc                                                                                                                                                   <  1:2380274/150‑1 (MQ=255)
       ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                                                                               >  2:3908569/1‑150 (MQ=255)
          tACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGt                                                                                                                                            <  1:3381727/150‑1 (MQ=255)
           aCATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTc                                                                                                                                           >  1:1134461/1‑150 (MQ=255)
             aTCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCtt                                                                                                                                         <  1:928219/150‑1 (MQ=255)
             aTCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCtt                                                                                                                                         <  1:1094978/150‑1 (MQ=255)
             aTCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACACCCGGTCtt                                                                                                                                         <  2:4428422/150‑1 (MQ=255)
               ccACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTgc                                                                                                                                       >  1:1080147/1‑150 (MQ=255)
                cACCGGTATGGGCTGGCTGCCACCCACCCTCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTgcg                                                                                                                                      <  1:1937564/150‑1 (MQ=255)
                    ggTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTgg                                                                                                                                  >  2:4886472/1‑150 (MQ=255)
                     gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                                                                                 <  1:1330509/150‑1 (MQ=255)
                     gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                                                                                 <  1:5106500/150‑1 (MQ=255)
                     gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                                                                                 >  1:4285380/1‑150 (MQ=255)
                       aTGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGa                                                                                                                               <  2:316839/150‑1 (MQ=255)
                       aTGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGa                                                                                                                               >  2:1286386/1‑150 (MQ=255)
                             tggctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCg                                                                                                                         >  2:4612366/1‑150 (MQ=255)
                                 tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                                                                     >  2:1041606/1‑150 (MQ=255)
                                   ccacccacCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGGCTTGCGCTGGTGCGCCCCGGTCAAg                                                                                                                   >  2:847367/1‑150 (MQ=255)
                                     acccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgc                                                                                                                 <  2:2960856/150‑1 (MQ=255)
                                      cccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcg                                                                                                                >  2:1573875/1‑150 (MQ=255)
                                      cccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcg                                                                                                                >  2:1850655/1‑150 (MQ=255)
                                       ccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcgc                                                                                                               >  1:2239164/1‑150 (MQ=255)
                                          cccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCt                                                                                                            >  1:1811512/1‑150 (MQ=255)
                                           ccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCtt                                                                                                           <  2:4973011/150‑1 (MQ=255)
                                            cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                                                                          >  2:2538998/1‑150 (MQ=255)
                                             gCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtct                                                                                                         <  1:3458948/150‑1 (MQ=255)
                                             gCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtct                                                                                                         >  2:3485753/1‑150 (MQ=255)
                                               tCGATGTGCTTCTCCTGCAAGAAACGCATCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCa                                                                                                       <  2:2465288/150‑1 (MQ=255)
                                                 gATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATg                                                                                                     >  1:4566574/1‑150 (MQ=255)
                                                      gCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGa                                                                                                >  1:140543/1‑150 (MQ=255)
                                                       cTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGaa                                                                                               <  1:3562036/150‑1 (MQ=255)
                                                         tctcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCGAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACa                                                                                             <  2:4393121/150‑1 (MQ=255)
                                                            ccTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGtt                                                                                          >  2:3756960/1‑150 (MQ=255)
                                                            ccTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGtt                                                                                          >  1:4562792/1‑150 (MQ=255)
                                                             cTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTc                                                                                         <  2:2018507/150‑1 (MQ=255)
                                                             cTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCATTCTCATGGCCGAACAGTTc                                                                                         >  1:1981135/1‑150 (MQ=255)
                                                              tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                                                                        >  1:2638464/1‑150 (MQ=255)
                                                              tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                                                                        <  2:4111686/150‑1 (MQ=255)
                                                              tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                                                                        >  1:2124301/1‑150 (MQ=255)
                                                              tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                                                                        >  2:437604/1‑150 (MQ=255)
                                                                 aacacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGAc                                                                                     <  1:2823338/145‑1 (MQ=255)
                                                                 aaGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGAc                                                                                     <  1:1078259/150‑1 (MQ=255)
                                                                 aaGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGAc                                                                                     <  2:1548444/150‑1 (MQ=255)
                                                                    acaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGACtcg                                                                                  >  2:3288087/3‑148 (MQ=255)
                                                                     aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                                                                                 >  1:1198042/1‑150 (MQ=255)
                                                                      aCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGAt                                                                                >  2:37849/1‑150 (MQ=255)
                                                                        gcagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCa                                                                              >  1:2419573/1‑150 (MQ=255)
                                                                        gcagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCa                                                                              >  2:1218966/1‑150 (MQ=255)
                                                                         cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                             <  2:1968957/150‑1 (MQ=255)
                                                                         cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                                                                             >  1:3113905/1‑150 (MQ=255)
                                                                          agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                                                                            >  1:3506049/1‑150 (MQ=255)
                                                                          agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                                                                            <  1:130801/150‑1 (MQ=255)
                                                                          agcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAg                                                                            >  1:1907044/1‑150 (MQ=255)
                                                                          agcaAGATGCCCAGCAATACAAACGGCAGCTCGCCGAATTCATCAAGAAACAAAGAGCCGCCAACGGCCGCTTCGATACGCCCGGACTAGCGCTGGTCAGCACAGCTAAAGGCGCCCTTCACAAGGCCGAACAGATCGCACAGGAACAAg                                                                            >  2:2654712/1‑150 (MQ=255)
                                                                              aGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCt                                                                        >  1:669605/1‑150 (MQ=255)
                                                                               gTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTc                                                                       <  2:930780/150‑1 (MQ=255)
                                                                                ttGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCg                                                                      >  1:3114448/1‑150 (MQ=255)
                                                                                ttGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCg                                                                      <  2:1756825/150‑1 (MQ=255)
                                                                                  ggCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCggg                                                                    >  1:2990506/1‑150 (MQ=255)
                                                                                  ggCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCggg                                                                    <  1:4159167/150‑1 (MQ=255)
                                                                                      tGCAATTCCAACGGCAGGTCGCCGATTTCTTCAGGTTAGAATGTGCCGCCAGCGGCCGCTTCGTTTCGCCGGGTCGTGCGCTGGTGAGCACCGGTCGTGGCGCCCGGCGCATGGCCGTACGGTGCGGACTGGATGGGGTGCTCggggggc                                                                >  2:720779/1‑147 (MQ=255)
                                                                                      tGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATc                                                                <  2:2784241/150‑1 (MQ=255)
                                                                                      tGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATc                                                                >  1:1797962/1‑150 (MQ=255)
                                                                                         aaTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcc                                                             >  1:876291/1‑150 (MQ=255)
                                                                                                 cGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGtt                                                     <  1:2551567/150‑1 (MQ=255)
                                                                                                     aGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATc                                                 <  1:4848736/150‑1 (MQ=255)
                                                                                                      ggTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATcg                                                <  1:978778/150‑1 (MQ=255)
                                                                                                         cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                             <  2:1431586/150‑1 (MQ=255)
                                                                                                         cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                             <  1:2626022/150‑1 (MQ=255)
                                                                                                         cGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGa                                             <  1:1491825/150‑1 (MQ=255)
                                                                                                           ccGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATg                                           >  2:2217996/1‑150 (MQ=255)
                                                                                                              aTTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAAc                                        >  1:2579043/1‑150 (MQ=255)
                                                                                                                ttCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACgg                                      <  2:3263917/150‑1 (MQ=255)
                                                                                                                 tcatcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgc                                     <  1:1311083/150‑1 (MQ=255)
                                                                                                                  catcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgct                                    <  1:1286386/150‑1 (MQ=255)
                                                                                                                  catcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgct                                    <  1:4149238/150‑1 (MQ=255)
                                                                                                                   atcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTAAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctg                                   >  1:61696/1‑150 (MQ=255)
                                                                                                                    tcaAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctgc                                  >  2:1103009/1‑150 (MQ=255)
                                                                                                                     caAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGgctgct                                 >  2:4460693/1‑150 (MQ=255)
                                                                                                                      aaGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTc                                >  1:59426/1‑150 (MQ=255)
                                                                                                                        gAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGc                              >  1:824196/1‑150 (MQ=255)
                                                                                                                          aagaaTGTGGGGCCATCGGGCGCTTCGATACGCCCGGTCTTGCGGTGGTGAGCACCGGTAATGGCGCCGTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCGCCGCTGTTGATCGCGATGAACGGGTGCTCGCTg                            >  2:3388507/4‑150 (MQ=255)
                                                                                                                           aCAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGc                           <  2:3609471/150‑1 (MQ=255)
                                                                                                                            cAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCg                          <  1:3756960/150‑1 (MQ=255)
                                                                                                                              ctgtgCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCGGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGtt                        <  2:3868219/149‑1 (MQ=255)
                                                                                                                                      ccATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTgg                <  1:4612366/150‑1 (MQ=255)
                                                                                                                                        aTCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCg              <  1:915473/150‑1 (MQ=255)
                                                                                                                                        aTCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCg              <  1:2124629/150‑1 (MQ=255)
                                                                                                                                          cGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGt            <  1:1573875/150‑1 (MQ=255)
                                                                                                                                            gCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGc          >  1:2861966/1‑150 (MQ=255)
                                                                                                                                              cGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCaa        <  2:2084764/150‑1 (MQ=255)
                                                                                                                                               gCTTCGATACGCCCGGTCTTGCGCTGGTGTGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCTGTTGTTCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg       >  1:3590835/1‑150 (MQ=255)
                                                                                                                                               gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg       <  1:3908569/150‑1 (MQ=255)
                                                                                                                                               gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg       >  1:3485719/1‑150 (MQ=255)
                                                                                                                                               gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg       <  1:3908424/150‑1 (MQ=255)
                                                                                                                                               gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg       <  1:1712448/150‑1 (MQ=255)
                                                                                                                                               gCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAg       >  1:3485761/1‑150 (MQ=255)
                                                                                                                                                cTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAgg      <  2:1080307/150‑1 (MQ=255)
                                                                                                                                                 ttCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGt     >  2:2621267/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAGTTGATCGCGATGAACGGCTGCTCGCTGCGTTGCGACTGGCGGTGCAAGGTGCG  >  NC_002947/3975613‑3975910

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