Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,255,892 C→A 6.3% intergenic (‑112/+24) PP_2852 ← / ← PP_2853 sulfatase domain‑containing protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,255,8920CA 6.3% 284.4 / 5.7 94intergenic (‑112/+24)PP_2852/PP_2853sulfatase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/42);  new base A (4/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAA  >  NC_002947/3255744‑3256038
                                                                                                                                                    |                                                                                                                                                  
cAATCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGca                                                                                                                                                   <  1:872625/150‑1 (MQ=255)
  aGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGcaca                                                                                                                                                 <  2:1803961/150‑1 (MQ=255)
     cGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGc                                                                                                                                              <  2:1587351/150‑1 (MQ=255)
          gCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGc                                                                                                                                         <  1:778793/150‑1 (MQ=255)
              gggTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCTGTGCGAACGGGCCCGGctttg                                                                                                                                     <  1:1683075/150‑5 (MQ=255)
              gggTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGAt                                                                                                                                     <  2:2008881/150‑1 (MQ=255)
                   gtCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCat                                                                                                                                <  2:914697/150‑2 (MQ=255)
                       gACAGGTTCTTCACAAAAGCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGTTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                            <  2:942368/150‑1 (MQ=255)
                        accggTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATc                                                                                                                           >  1:2401814/4‑150 (MQ=255)
                         cAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCa                                                                                                                          <  2:2145037/150‑1 (MQ=255)
                         cAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCa                                                                                                                          <  1:2037197/150‑1 (MQ=255)
                            gTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGaa                                                                                                                       >  1:2135946/1‑150 (MQ=255)
                            gTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGaa                                                                                                                       >  2:534645/1‑150 (MQ=255)
                            gTTCTTCACAAAAACTTACCTTTTGCCATTAAATCGTCATGCATCTTTGCTTGATAATATCAGGCCAGGGCACGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGCGCCCGTTCGCACAGGCCAGGCCGATGGAAACTACCAgca                                                                                                                       >  2:467459/1‑148 (MQ=255)
                                   acaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                <  1:1971376/150‑1 (MQ=255)
                                   acaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                <  1:93775/150‑1 (MQ=255)
                                    caAAAACTTCCCTTTTGCCATTAAATCGTCGTGCATCTGTGCTTGATGATATCAGGACAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               >  1:489548/1‑150 (MQ=255)
                                    caAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTCCCGCGCGATTGCCGCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               >  2:1269317/1‑150 (MQ=255)
                                     aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAgg                                                                                                              <  1:534645/150‑1 (MQ=255)
                                               cTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                                    >  1:529379/1‑150 (MQ=255)
                                                     ccATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              <  2:1369906/150‑1 (MQ=255)
                                                     ccATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCAGTTCGCACAGTCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              <  2:1181875/150‑1 (MQ=255)
                                                       aTTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                                    >  2:1747443/1‑142 (MQ=255)
                                                       aTTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                                    <  1:1747443/142‑1 (MQ=255)
                                                       aTTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaaga                                                                                            >  1:316332/1‑150 (MQ=255)
                                                           aaTCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                         <  1:2042106/149‑1 (MQ=255)
                                                           aaTCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                         >  2:2042106/1‑149 (MQ=255)
                                                            aTCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGt                                                                                       >  2:958189/1‑150 (MQ=255)
                                                              cGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCa                                                                                     >  2:2228674/1‑150 (MQ=255)
                                                               gTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAg                                                                                    >  2:2207859/1‑150 (MQ=255)
                                                                tCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAgg                                                                                   <  1:958189/150‑1 (MQ=255)
                                                                 cATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                                  <  2:2135946/150‑1 (MQ=255)
                                                                  aTGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGcaca                                                                                                                                                 <  2:227306/86‑1 (MQ=255)
                                                                  aTGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGcaca                                                                                                                                                 >  1:227306/1‑86 (MQ=255)
                                                                   tGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATgg                                                                                                                                   >  1:1893602/1‑99 (MQ=255)
                                                                   tGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATgg                                                                                                                                   <  2:1893602/99‑1 (MQ=255)
                                                                      aTCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCTTc                                                                             <  2:2401814/150‑1 (MQ=255)
                                                                      aTCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTc                                                                             >  2:1028891/1‑150 (MQ=255)
                                                                         tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTCACGAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATc                                                                          >  1:297397/1‑150 (MQ=255)
                                                                         tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAg                                                                                                     <  2:1427055/123‑1 (MQ=255)
                                                                         tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAg                                                                                                     >  1:1427055/1‑123 (MQ=255)
                                                                         tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATAAGCTGGTGAAAGAAAAGTCAGGTGCGTCATc                                                                          <  1:1269317/150‑1 (MQ=255)
                                                                             cTTGATGATATAAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTc                                                                      <  2:609741/150‑1 (MQ=255)
                                                                              ttGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                >  1:58145/1‑107 (MQ=255)
                                                                              ttGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                <  2:58145/107‑1 (MQ=255)
                                                                              ttGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCa                                                                     >  2:653628/1‑150 (MQ=255)
                                                                                gatgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGa                                                                   >  2:916752/1‑150 (MQ=255)
                                                                                gatgatATCAGGCCAGGGCAGGGAAGGCGGCCAGAATTCATACTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGa                                                                   <  1:1331707/150‑1 (MQ=255)
                                                                                 atgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                         >  2:2355919/1‑127 (MQ=255)
                                                                                 atgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                         <  1:2355919/127‑1 (MQ=255)
                                                                                  tgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGAcc                                                                 >  2:892840/1‑150 (MQ=255)
                                                                                        cAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGa                                                                                                                                      <  2:408729/75‑1 (MQ=255)
                                                                                        cAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGa                                                                                                                                      >  1:408729/1‑75 (MQ=255)
                                                                                        cAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAAttgt                                                           <  1:467459/150‑4 (MQ=255)
                                                                                         aGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGTCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGAc                                                          <  1:653628/150‑1 (MQ=255)
                                                                                          ggCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACa                                                         <  1:1028891/150‑1 (MQ=255)
                                                                                             cagggcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGt                                                      >  1:1222492/1‑150 (MQ=255)
                                                                                               gggcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGTGTCATCGGTCAGACCAATGGACAGGTTg                                                    <  1:1603496/150‑1 (MQ=255)
                                                                                                   agggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCa                                                >  1:2275626/1‑150 (MQ=255)
                                                                                                      gAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGa                                             <  1:2228674/150‑1 (MQ=255)
                                                                                                        aGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCa                                                                                     >  2:554045/1‑108 (MQ=255)
                                                                                                        aGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCa                                                                                     <  1:554045/108‑1 (MQ=255)
                                                                                                            ctttggTTACAAGCCTTGTAACTAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCa                                       >  2:1663396/6‑150 (MQ=255)
                                                                                                             gcCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCATAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAg                                      <  2:529379/150‑1 (MQ=255)
                                                                                                               cAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGc                                                                                                                                  <  1:1934524/56‑1 (MQ=255)
                                                                                                               cAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGc                                                                                                                                  >  2:1934524/1‑56 (MQ=255)
                                                                                                                ggTTACAAGCCTTGTAACCAAAGCCGGGCCCGTTCGCGCAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAACGTCTGGGGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGtt                                   >  1:1954949/2‑150 (MQ=255)
                                                                                                                  ttACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTCAAAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGCCAGGCTGGGCACGATGGTCAGGTTGt                                 >  1:525172/1‑150 (MQ=255)
                                                                                                                   tACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTc                                <  2:1222492/150‑1 (MQ=255)
                                                                                                                     cAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCgg                              >  1:1642708/1‑150 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCa                                                                                                                          <  1:1306359/57‑1 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCa                                                                                                                          >  2:1306359/1‑57 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                             >  1:1691865/1‑150 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                             <  2:1642708/150‑1 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                             >  1:690865/1‑150 (MQ=255)
                                                                                                                       aGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTg                            >  1:1342102/1‑150 (MQ=255)
                                                                                                                           ttGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTggcgg                        >  1:707260/1‑150 (MQ=255)
                                                                                                                            tGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTggcggc                       >  1:1281503/1‑150 (MQ=255)
                                                                                                                              tAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGa                     >  2:2259708/1‑150 (MQ=255)
                                                                                                                                 ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAg                                                        <  1:1440651/112‑1 (MQ=255)
                                                                                                                                 ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAg                                                        >  2:1440651/1‑112 (MQ=255)
                                                                                                                                 ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAgg                  >  2:171199/1‑150 (MQ=255)
                                                                                                                                   tcggCCGGGCCTGTGCGAACGGGCCCGGCTGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATg                                                             >  2:123649/4‑105 (MQ=37)
                                                                                                                                   tcggCCGGGCCTGTGCGAACGGGCCCGGCTGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATg                                                             <  1:123649/102‑1 (MQ=37)
                                                                                                                                   tcggCCGGGCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTc                >  2:1962267/4‑150 (MQ=255)
                                                                                                                                   aaaGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTc                >  1:2130828/1‑150 (MQ=255)
                                                                                                                                     aGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAg              <  2:1342102/150‑1 (MQ=255)
                                                                                                                                        cGGGCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGcc           >  2:555521/1‑150 (MQ=255)
                                                                                                                                        cGGGCCTGTGCGAACGGGCCCGGCCGATGGCAAATATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGcc           >  2:177098/1‑150 (MQ=255)
                                                                                                                                        cGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGcc           >  1:966347/1‑150 (MQ=255)
                                                                                                                                          ggCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCaa         >  1:1920030/1‑150 (MQ=255)
                                                                                                                                           gCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATTGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAat        >  2:881666/1‑150 (MQ=255)
                                                                                                                                            cccGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAata       <  1:892840/150‑1 (MQ=255)
                                                                                                                                             ccGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATg                                           >  2:2185222/1‑113 (MQ=255)
                                                                                                                                             ccGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATg                                           <  1:2185222/113‑1 (MQ=255)
                                                                                                                                               gTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTg    >  1:2394458/1‑150 (MQ=255)
                                                                                                                                                 tCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa  <  2:2130828/150‑1 (MQ=255)
                                                                                                                                                 tCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa  <  1:916752/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAA  >  NC_002947/3255744‑3256038

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