Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,015,961 A→T 9.8% intergenic (‑176/+58) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,015,9610AT 9.8% 234.7 / 13.2 82intergenic (‑176/+58)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base T (4/4);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

ATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGAGGA  >  NC_002947/5015816‑5016109
                                                                                                                                                   |                                                                                                                                                      
atgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTGTTACTTTCCGGGCCTTGGACAGTAGCTGATCCGAAAGTACGTCCTTGTGAACTGATGCCATCGTCTTGGCGTTGCTTGTGATAGCGCTACGAAAGAAAAAAGCCCCCCTCAAgtc                                                                                                                                                    >  1:1088872/1‑148 (MQ=255)
     tAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTAGATTGAGATAGCACTACAAAAAAAAAAAAACCGCCCTTAGGGGGGCt                                                                                                                                               >  1:773236/1‑150 (MQ=255)
      aaGCCCGAACCGTTCTCCGGCGCGGGCTTTTATTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCTTTGTGAACTAATGACAGCGTTTTGGTTTTGATTGTGATAGAATTAAAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtt                                                                                                                                              <  2:773236/150‑1 (MQ=255)
       aGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCttt                                                                                                                                             <  2:1455198/150‑2 (MQ=255)
        gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGACCGCCTAAAGGCGGGCTTaa                                                                                                                                            >  2:930568/1‑150 (MQ=255)
             cACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGTACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTAATTGTGATATCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGa                                                                                                                                       <  1:738910/150‑1 (MQ=255)
                     tcCGGCGCGCGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTACTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAACCCGCACAAACGCGGGCTTTACACGAAGCGGGTa                                                                                                                               >  1:2276533/1‑150 (MQ=255)
                       cGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAAt                                                                                                                             >  1:908218/1‑150 (MQ=255)
                             gggCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGGTAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGGAATTCAGCTg                                                                                                                      >  1:1969986/1‑150 (MQ=255)
                                    tACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGc                                                                                                                                                    <  1:1103184/114‑1 (MQ=255)
                                    tACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCCAAAGGc                                                                                                                                                    >  2:1103184/1‑114 (MQ=255)
                                      gtttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt                                                                                                              >  2:1655995/2‑150 (MQ=255)
                                       tttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                             >  1:1814152/1‑150 (MQ=255)
                                           cGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGcc                                                                                                           <  1:155610/148‑1 (MQ=255)
                                           cGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGcc                                                                                                           >  2:155610/1‑148 (MQ=255)
                                                 ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGa                                                                                                   >  2:2107651/1‑150 (MQ=255)
                                                  tGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                  <  1:930568/150‑1 (MQ=255)
                                                  tGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTCGATGCGGGCTTAAAACGCAGCGGGTAATTGAGCTGGCATCCTTGCCAACGTTGATg                                                                                                 >  1:636001/1‑150 (MQ=255)
                                                   ggACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                  <  1:1434058/150‑1 (MQ=255)
                                                      cAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCACCTACAGGCGTGCTTAAAAAGAAAAGGGGAATTTAGCTGGCCTCATTACAAAAGTTGagcgcc                                                                                              >  2:1170883/1‑145 (MQ=255)
                                                         tATCTGACACGAAAGTACGTCATTGTAAAATGATAACAGCGTATTAGATTTAATTATAAAAAAAAAAAAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCa                                                                                            <  2:1359207/150‑1 (MQ=255)
                                                          aTCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGa                                                                                          <  1:491223/150‑1 (MQ=255)
                                                           tCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa                                                                                                                                      <  2:1409405/105‑1 (MQ=255)
                                                           tCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa                                                                                                                                      >  1:1409405/1‑105 (MQ=255)
                                                              gACCCGAAAGTACGTCCTTGTGAACTGATGACATCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCtt                                                                                      <  1:198609/150‑1 (MQ=255)
                                                              gACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCtt                                                                                      >  2:1163119/1‑150 (MQ=255)
                                                              gACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCtt                                                                                      >  2:1033107/1‑150 (MQ=255)
                                                              gACCCGAAAGTACGTCCTTGTAAACTGATGACAGCGTCTTGACTTTGATTGTGATAAAAATAAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCtt                                                                                      <  2:2276533/150‑1 (MQ=255)
                                                               aCCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                             <  1:1192380/126‑1 (MQ=255)
                                                               aCCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                             >  2:1192380/1‑126 (MQ=255)
                                                                 ccGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTTAAACGAAGAGGGTAATTCAGCTGGTATCATTGCCAACGTTGATGGCTCAGACCTTCAg                                                                                   >  1:558250/1‑150 (MQ=255)
                                                                 ccGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAg                                                                                   <  2:908218/150‑1 (MQ=255)
                                                                     aaGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtca                                                                              <  1:13538/150‑1 (MQ=255)
                                                                       gTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtcatc                                                                            <  1:1554786/150‑1 (MQ=255)
                                                                           gTCCTTGTGAACTGATGACAGCCTCTTGGCTTTGATTGTGAAAAAAATAAAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAg                                                                         <  2:33734/150‑1 (MQ=255)
                                                                            tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCATAAAGGGGGGCTTAAAACGAAAAGGGTAATTCAGCTGGCATCATTTCCAACTTTGATGGGTCAGACATTCAGCTCATCAAgg                                                                        >  2:1180216/1‑150 (MQ=255)
                                                                            tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGa                                                                       >  1:2357796/1‑150 (MQ=255)
                                                                            tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGa                                                                       <  1:2107651/150‑1 (MQ=255)
                                                                            tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGACCGCCTCAAGGAGGGCTTAAAAAGAAGAGGGTAATTCAGCTTGCGTCATTTACAACCTTTATGTCTCTCACCGTTTGCTCTTCAAAGa                                                                       >  2:475164/1‑150 (MQ=255)
                                                                               ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGc                                                                    >  1:1428992/1‑150 (MQ=255)
                                                                                 cgatcttgatgaCAGCGTCTTGGCTTTGATTGCAGCTTCATGACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTc                                                                 <  2:1998864/144‑1 (MQ=255)
                                                                                    aaCTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTTCCAACGTTGATTGCTCAGACCTTCAGCTCATCAAGGACGCTCAGc                                                              >  1:2002236/1‑150 (MQ=255)
                                                                                       tgatgaCAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAAAGAAAAGGGTAATTCAGCTGACATCATTGCCAACGTTGATGGCTCAGACATTCAGCTCATCAAGGACGCTCAGaacc                                                           >  2:362094/1‑149 (MQ=255)
                                                                                        gatgaCAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCaagatcg                                                          >  1:1998864/1‑144 (MQ=255)
                                                                                               gCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgt                                                    <  1:342312/150‑1 (MQ=255)
                                                                                               gCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtg                                                   >  1:1847155/1‑150 (MQ=255)
                                                                                               gCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtg                                                   <  1:1033107/150‑1 (MQ=255)
                                                                                                   cTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAg                                                                         >  2:2039410/1‑125 (MQ=255)
                                                                                                   cTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAg                                                                         <  1:2039410/125‑1 (MQ=255)
                                                                                                       gCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGc                                           <  2:2357796/150‑1 (MQ=255)
                                                                                                        cTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTTATGGCTCAGACCTTCAGCTCCTCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCt                                          >  1:621630/1‑150 (MQ=255)
                                                                                                         tttGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAACTGGCATCAGTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtg                                         >  2:846292/1‑150 (MQ=255)
                                                                                                          ttGATTGTGATAGCACGACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtga                                        <  1:99373/150‑1 (MQ=255)
                                                                                                            gATTGTGATAGCACTACAAAAAAAAAAACCCCCCTAAAAGGGGGGTTTAAAACAAAGAGGGAAATTCAGTTGCCATATTTGCCAACTTTATTGGCTAAGCCCTTAGCTCCATAAAGGCGCTTCAGCCCACCTTCTGGTGCGATgtggata                                      >  2:2033197/1‑144 (MQ=255)
                                                                                                             aTTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtgatga                                     >  1:1655809/1‑150 (MQ=255)
                                                                                                                gtgATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGc                                  <  1:1655995/150‑1 (MQ=255)
                                                                                                                  aagcgaaCTACAAAAAAAAAAAGCCCGCCTAAAGGCAGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGa                                <  1:475164/144‑1 (MQ=255)
                                                                                                                     aGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATcc                             >  2:795961/1‑150 (MQ=255)
                                                                                                                     aGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATcc                             >  2:990764/1‑150 (MQ=255)
                                                                                                                      gCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCa                            >  1:1511326/1‑150 (MQ=255)
                                                                                                                          tACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTc                        <  1:990764/150‑1 (MQ=255)
                                                                                                                            aaaaaaaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTc                        <  2:558250/147‑1 (MQ=255)
                                                                                                                            aCAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtca                                                                              >  1:928965/1‑96 (MQ=255)
                                                                                                                            aCAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtca                                                                              <  2:928965/96‑1 (MQ=255)
                                                                                                                               aaaaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGa                                <  1:1081820/139‑1 (MQ=255)
                                                                                                                               aaaaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGa                                >  2:1081820/1‑139 (MQ=255)
                                                                                                                                 aaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGtt                   <  1:934587/150‑1 (MQ=255)
                                                                                                                                 aaaaaaaaaGCCCGCCTAAAGGCGGGCATAAAAAGAAGAGGTTAAGTCAGCGGCCATCATTGCCAACGTGGAGGGCTCCGACCTTCAGCGCATCAACGCCGCTCCGCACAGCTTCGTGTGCGATTCTGATGATGCGATGCTCTTCCTGtg                   <  2:1088872/150‑2 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCgg                >  1:261669/1‑150 (MQ=255)
                                                                                                                                     tttaaGCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCa                            <  1:732065/134‑1 (MQ=255)
                                                                                                                                     tttaaGCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCa                            >  2:732065/4‑137 (MQ=255)
                                                                                                                                      aaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTc              <  1:795961/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGt            <  2:1847155/150‑1 (MQ=255)
                                                                                                                                         aGCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTg           >  1:378175/1‑150 (MQ=255)
                                                                                                                                          gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATgg                                                                                                >  2:646169/1‑66 (MQ=39)
                                                                                                                                          gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATgg                                                                                                <  1:646169/66‑1 (MQ=39)
                                                                                                                                            ccGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACa        >  1:525604/1‑150 (MQ=255)
                                                                                                                                             cGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGCTGATGCGAGCCACTTCCTGTTCGGTCGTGACaa       <  1:1270090/150‑1 (MQ=255)
                                                                                                                                              gCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAag      <  2:1814152/150‑1 (MQ=255)
                                                                                                                                               ccTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAaga     <  1:242707/150‑1 (MQ=255)
                                                                                                                                                cTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAagag    <  2:1655809/150‑1 (MQ=255)
                                                                                                                                                 tAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAgagg   >  2:1472483/1‑150 (MQ=255)
                                                                                                                                                  aaaGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAgagga  <  2:1511326/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGAGGA  >  NC_002947/5015816‑5016109

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