Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,255,886 C→T 6.4% intergenic (‑106/+30) 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,8860CT 6.4% 269.2 / 6.8 94intergenic (‑106/+30)PP_2852/PP_2853sulfatase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/41);  new base T (4/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATA  >  NC_002947/3255744‑3256033
                                                                                                                                              |                                                                                                                                                   
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: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)
                                      aaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTc                                                                                                                                                 >  2:732228/1‑109 (MQ=255)
                                      aaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTc                                                                                                                                                 <  1:732228/109‑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                                                                                                           <  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                                                                                    >  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                                                                                                                     >  2:1306359/1‑57 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCa                                                                                                                     <  1:1306359/57‑1 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                        >  1:690865/1‑150 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                        >  1:1691865/1‑150 (MQ=255)
                                                                                                                      aaGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                        <  2:1642708/150‑1 (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                                                        <  1:123649/102‑1 (MQ=37)
                                                                                                                                   tcggCCGGGCCTGTGCGAACGGGCCCGGCTGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATg                                                        >  2:123649/4‑105 (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)
                                                                                                                                              |                                                                                                                                                   
CAAGCCGGGTGCTTGGGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATA  >  NC_002947/3255744‑3256033

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: