Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,255,895 A→G 6.4% intergenic (‑115/+21) 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,8950AG 6.4% 273.1 / 6.8 94intergenic (‑115/+21)PP_2852/PP_2853sulfatase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/42);  new base G (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 = 3.21e-01

Only 100 of 102 total aligned reads displayed.

AGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAAGTTGGC  >  NC_002947/3255746‑3256044
                                                                                                                                                     |                                                                                                                                                     
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:93775/150‑1 (MQ=255)
                                 acaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                      <  1:1971376/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                                                                                                          <  1:1747443/142‑1 (MQ=255)
                                                     aTTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                                          >  2:1747443/1‑142 (MQ=255)
                                                     aTTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaaga                                                                                                  >  1:316332/1‑150 (MQ=255)
                                                         aaTCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                               <  1:2042106/149‑1 (MQ=255)
                                                         aaTCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                               >  2:2042106/1‑149 (MQ=255)
                                                            cGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCa                                                                                           >  2:2228674/1‑150 (MQ=255)
                                                             gTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAg                                                                                          >  2:2207859/1‑150 (MQ=255)
                                                               cATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                                        <  2:2135946/150‑1 (MQ=255)
                                                                aTGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGcaca                                                                                                                                                       >  1:227306/1‑86 (MQ=255)
                                                                aTGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGcaca                                                                                                                                                       <  2:227306/86‑1 (MQ=255)
                                                                 tGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATgg                                                                                                                                         <  2:1893602/99‑1 (MQ=255)
                                                                 tGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATgg                                                                                                                                         >  1:1893602/1‑99 (MQ=255)
                                                                    aTCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCTTc                                                                                   <  2:2401814/150‑1 (MQ=255)
                                                                    aTCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTc                                                                                   >  2:1028891/1‑150 (MQ=255)
                                                                       tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTCACGAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATc                                                                                >  1:297397/1‑150 (MQ=255)
                                                                       tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAg                                                                                                           >  1:1427055/1‑123 (MQ=255)
                                                                       tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAg                                                                                                           <  2:1427055/123‑1 (MQ=255)
                                                                       tgtgCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATAAGCTGGTGAAAGAAAAGTCAGGTGCGTCATc                                                                                <  1:1269317/150‑1 (MQ=255)
                                                                           cTTGATGATATAAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTc                                                                            <  2:609741/150‑1 (MQ=255)
                                                                            ttGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                      <  2:58145/107‑1 (MQ=255)
                                                                            ttGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                      >  1:58145/1‑107 (MQ=255)
                                                                            ttGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCa                                                                           >  2:653628/1‑150 (MQ=255)
                                                                              gatgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGa                                                                         >  2:916752/1‑150 (MQ=255)
                                                                              gatgatATCAGGCCAGGGCAGGGAAGGCGGCCAGAATTCATACTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGa                                                                         <  1:1331707/150‑1 (MQ=255)
                                                                               atgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                               <  1:2355919/127‑1 (MQ=255)
                                                                               atgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGaaaa                                                                                               >  2:2355919/1‑127 (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:690865/1‑150 (MQ=255)
                                                                                                                    aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                                   <  2:1642708/150‑1 (MQ=255)
                                                                                                                    aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                                   >  1:1691865/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                                                 <  1:2185222/113‑1 (MQ=255)
                                                                                                                                           ccGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATg                                                 >  2:2185222/1‑113 (MQ=255)
                                                                                                                                             gTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTg          >  1:2394458/1‑150 (MQ=255)
                                                                                                                                               tCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa        <  2:2130828/150‑1 (MQ=255)
                                                                                                                                               tCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa        <  1:916752/150‑1 (MQ=255)
                                                                                                                                                   acaGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGGATTTCAGGCTGAAGGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAAGTTg    >  2:1727549/1‑150 (MQ=255)
                                                                                                                                                   acaGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa        >  2:2274015/1‑146 (MQ=255)
                                                                                                                                                   acaGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa        <  1:2274015/146‑1 (MQ=255)
                                                                                                                                                   acaGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGGTGTTGGTGGCGGCGAAGGTCAGGCCCATATTGACGTTg    >  1:277816/1‑150 (MQ=255)
                                                                                                                                                     aGGGCCGGCCGATTACAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGAACAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAAGTTggc  <  2:297397/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAAGTTGGC  >  NC_002947/3255746‑3256044

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: