Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,398,502 G→T 5.3% G489V (GGA→GTA)  hsdS → type I restriction modification system specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,398,5020GT 5.3% 236.1 / 2.5 76G489V (GGA→GTA) hsdStype I restriction modification system specificity protein
Reads supporting (aligned to +/- strand):  ref base G (31/40);  new base T (2/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

CAGAGACTAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACC  >  NC_002947/5398355‑5398651
                                                                                                                                                   |                                                                                                                                                     
cAGAGACTAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCTGTTTTCATCCTGTTGAGCTGGGTGTTTATGGaa                                                                                                                                                     <  1:461002/150‑1 (MQ=255)
 agagACTAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAg                                                                                                                                                    <  2:692506/150‑1 (MQ=255)
   agaCTAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGt                                                                                                                                                  >  1:348003/1‑150 (MQ=255)
    gaCTAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGtt                                                                                                                                                 >  1:912768/1‑150 (MQ=255)
       tAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGt                                                                                                                                              <  2:1953130/150‑1 (MQ=255)
           gTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGct                                                                                                                                          >  2:1650365/1‑150 (MQ=255)
            tCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGctc                                                                                                                                         <  2:1489376/150‑1 (MQ=255)
                   tataCTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCtt                                                                                                                                  <  2:352909/150‑1 (MQ=255)
                   tataCTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCtt                                                                                                                                  <  1:1974386/150‑1 (MQ=255)
                    ataCTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGctct                                                                                                                                        <  1:2243228/143‑1 (MQ=255)
                    ataCTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGctct                                                                                                                                        >  2:2243228/1‑143 (MQ=255)
                     taCTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGa                                                                                                                                <  2:2316970/150‑1 (MQ=255)
                         aTTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGc                                                                                                                            >  2:2263691/1‑150 (MQ=255)
                           tCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTg                                                                                                                          <  2:1436841/150‑1 (MQ=255)
                            cACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTgg                                                                                                                         <  1:1965935/150‑1 (MQ=255)
                             aCTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGga                                                                                                                        >  2:1911019/1‑150 (MQ=255)
                             aCTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGga                                                                                                                        <  2:1974823/150‑1 (MQ=255)
                               tCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGgaga                                                                                                                      >  1:650846/1‑150 (MQ=255)
                                  ggCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAAc                                                                                                                   <  2:348003/150‑1 (MQ=255)
                                         ccTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCaaa                                                                                                            >  1:1172468/1‑150 (MQ=255)
                                             tGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATAtc                                                                                                        <  2:650846/150‑1 (MQ=255)
                                              gAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGAGCATATGCAAGGGTTTTCGCTCTCGGCTTGAACGCTGGAGAAACTGCCAAATAtct                                                                                                       >  2:62686/1‑150 (MQ=255)
                                              gAACAGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGCGAAACTGCCAAATAtct                                                                                                       >  2:286521/1‑150 (MQ=255)
                                               aaCCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATAtctc                                                                                                      <  2:1036103/150‑1 (MQ=255)
                                               aaCCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATAtctc                                                                                                      <  1:1655940/150‑1 (MQ=255)
                                               aaCCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATAtctc                                                                                                      <  2:860814/150‑1 (MQ=255)
                                                aCCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCg                                                                                                     >  2:2392119/1‑150 (MQ=255)
                                                     gTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAg                                                                                                <  1:2224880/150‑1 (MQ=255)
                                                     gTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAg                                                                                                <  2:1765817/150‑1 (MQ=255)
                                                      tGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAgg                                                                                               <  1:1911019/150‑1 (MQ=255)
                                                       gggAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGc                                                                                              >  2:125346/1‑150 (MQ=255)
                                                        ggAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGcc                                                                                             <  1:2263691/150‑1 (MQ=255)
                                                               cTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTACGGAAGGTATGTCGCTCTTTGCTCGAACGCTGGAGAAACTTCGAAATGTCTCGAACAGGCCAAGTCCg                                                                                      <  2:868372/150‑1 (MQ=255)
                                                                  gCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCa                                                                                   >  1:752309/1‑150 (MQ=255)
                                                                     tCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGaaa                                                                                                                    >  2:234018/1‑114 (MQ=255)
                                                                     tCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGaaa                                                                                                                    <  1:234018/114‑1 (MQ=255)
                                                                     tCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATgg                                                                                <  1:2392119/150‑1 (MQ=255)
                                                                       gTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCttt                                                                                                                                      >  1:1610845/1‑94 (MQ=255)
                                                                       gTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCttt                                                                                                                                      <  2:1610845/94‑1 (MQ=255)
                                                                       gTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATggtt                                                                              <  2:1172468/150‑3 (MQ=255)
                                                                           aaaGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCaa                                                                          <  2:1164417/150‑1 (MQ=255)
                                                                             aGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGt                                                                        <  1:125346/150‑1 (MQ=255)
                                                                                      ccGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGATCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAat                                                               >  1:416332/1‑150 (MQ=255)
                                                                                      ccGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCCAGACCCGCATGGCTGCAAGTCAGGTCAat                                                               >  2:1135590/1‑150 (MQ=255)
                                                                                      ccGAAATTAATGACTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAat                                                               <  2:551215/150‑1 (MQ=255)
                                                                                       cGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAata                                                              >  1:39817/1‑150 (MQ=255)
                                                                                             tAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCttt                                                                                                                                      <  1:1510231/72‑1 (MQ=255)
                                                                                             tAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCttt                                                                                                                                      >  2:1510231/1‑72 (MQ=255)
                                                                                               aTGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCaa                                                      >  1:1640482/1‑150 (MQ=255)
                                                                                                     tCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAg                                                >  1:2175205/1‑150 (MQ=255)
                                                                                                         aCAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTga                                            <  2:2175205/150‑1 (MQ=255)
                                                                                                         aCAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTga                                            <  2:883517/150‑1 (MQ=255)
                                                                                                         aCAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCGCTTTGCTTGAACGCTGGAGAAACTGCCAAATAGCTCGACTAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTga                                            <  2:912768/150‑1 (MQ=255)
                                                                                                             gATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGAtt                                        >  2:1005551/1‑150 (MQ=255)
                                                                                                                catcaGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCaa                                                                          >  2:1044075/1‑113 (MQ=255)
                                                                                                                catcaGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCaa                                                                          <  1:1044075/113‑1 (MQ=255)
                                                                                                                catcaGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGc                                     >  1:1164305/1‑150 (MQ=255)
                                                                                                                catcaGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGc                                     >  1:2435399/1‑150 (MQ=255)
                                                                                                                   caGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAata                                                                                                          >  2:2296154/1‑78 (MQ=255)
                                                                                                                   caGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAata                                                                                                          <  1:2296154/78‑1 (MQ=255)
                                                                                                                       tttCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATc                              <  2:609325/150‑1 (MQ=255)
                                                                                                                       tttCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATc                              >  1:1258313/1‑150 (MQ=255)
                                                                                                                       tttCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATc                              >  2:878189/1‑150 (MQ=255)
                                                                                                                            tCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCt                                                                              >  2:333359/1‑97 (MQ=255)
                                                                                                                            tCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCt                                                                              <  1:333359/97‑1 (MQ=255)
                                                                                                                            tCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGtt                         <  1:2080682/150‑1 (MQ=255)
                                                                                                                               tGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGc                                                                            >  1:677278/1‑96 (MQ=255)
                                                                                                                               tGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGc                                                                            <  2:677278/96‑1 (MQ=255)
                                                                                                                                gTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCa                     >  1:2095868/1‑150 (MQ=255)
                                                                                                                                   gAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGgaga                                                                                                                      >  1:1599303/1‑50 (MQ=255)
                                                                                                                                   gAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGgaga                                                                                                                      <  2:1599303/50‑1 (MQ=255)
                                                                                                                                            gtTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAg         <  1:1717674/150‑1 (MQ=255)
                                                                                                                                               tATGTAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTc                                                                       <  2:1600087/85‑1 (MQ=255)
                                                                                                                                               tATGTAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTc                                                                       >  1:1600087/1‑85 (MQ=255)
                                                                                                                                               tATGTAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTc                                                                       >  1:1600093/1‑85 (MQ=255)
                                                                                                                                               tATGTAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTc                                                                       <  2:1600093/85‑1 (MQ=255)
                                                                                                                                                  ggAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAacac   >  1:847173/1‑150 (MQ=255)
                                                                                                                                                   gAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGGTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACAcc  >  2:438658/1‑150 (MQ=255)
                                                                                                                                                   gAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACAcc  >  1:774930/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGAGACTAAGGTCGGTGATATACTATTCACTCGGGCTGGCCCTATGAACCGCGTGGGAATTTCTTGCCTCGTTGAAAGTACGCGCCCGAAATTAATGATTTCCGACAAGATCATCAGGTTTCATCCTGTTGAGCTGGGTGTTTATGGAAGGTTTGTCGCTCTTTGCTTGAACGCTGGAGAAACTGCCAAATATCTCGAACAGGCCAAGTCCGGCATGGCTGCAAGTCAGGTCAATATCTCGCAAGAAAAGCTGAGATTGGCTCCGATCCCGTTAGCACCGCTCCGCGAGCAACACC  >  NC_002947/5398355‑5398651

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