Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,815,600 A→G 7.4% intergenic (+28/+63) PP_5090 → / ← PP_5091 cell division protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,815,6000AG 7.4% 232.2 / 5.7 81intergenic (+28/+63)PP_5090/PP_5091cell division protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (44/31);  new base G (3/3);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGAACAA  >  NC_002947/5815455‑5815740
                                                                                                                                                 |                                                                                                                                            
gTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                          <  1:855645/150‑1 (MQ=255)
      cgTGCTTGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTTCCGCAGCGATGCGACAACGCcttc                                                                                                                                    <  2:253116/150‑1 (MQ=255)
       gTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCCGCAACCTGTTGCTGCAACAGCGGCAGATCCAGCAGTAACGTATAAGGCCTTGCCGCAGCGATGCGACAAGGCcttct                                                                                                                                   >  1:435691/1‑150 (MQ=255)
       gTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCcttct                                                                                                                                   >  2:1285326/1‑150 (MQ=255)
        tGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCcttctt                                                                                                                                  <  1:1157025/150‑1 (MQ=255)
          cTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTACTGCAACAGCGGCAGACCCGTCAGTAACGAATAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTta                                                                                                                                >  2:405617/1‑150 (MQ=255)
              ttGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCATCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCGTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATg                                                                                                                            >  2:587087/1‑150 (MQ=255)
                 ggCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACaa                                                                                                                                           >  1:1889466/1‑132 (MQ=255)
                 ggCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACaa                                                                                                                                           <  2:1889466/132‑1 (MQ=255)
                  gCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCCGCGATGCGACAAGGCCTTCTTTATATGGACa                                                                                                                        >  1:2225/1‑150 (MQ=255)
                  gCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACa                                                                                                                        >  2:1444854/1‑150 (MQ=255)
                       tCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCAGGGTTCATCAACCTGTTGGTTCGACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCAtt                                                                                                                   <  1:228436/150‑1 (MQ=255)
                        cAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTg                                                                                                                  <  2:703199/150‑1 (MQ=255)
                           cAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtat                                                                                                                               <  1:952437/134‑1 (MQ=255)
                           cAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtat                                                                                                                               >  2:952437/1‑134 (MQ=255)
                           cAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAAACTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATAGACAGCATTGcac                                                                                                               >  2:101938/1‑148 (MQ=255)
                                 cgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCAGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACACGGCCTTCCCTATAGGGACAGCATTGCCCGTGAAg                                                                                                         <  1:1500177/150‑1 (MQ=255)
                                  gAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAgg                                                                                                        <  2:1731409/150‑1 (MQ=255)
                                      aGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCt                                                                                                    >  2:625935/1‑150 (MQ=255)
                                           aCAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTaaa                                                                                               >  1:78336/1‑150 (MQ=255)
                                           aCAGGGGAACAGAAGCAGCTGGCTGGGGCCGGGTAAGGGAATCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGTGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTaaa                                                                                               <  1:101938/150‑1 (MQ=255)
                                              ggggCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagca                                                                                            <  1:1444854/150‑1 (MQ=255)
                                              ggggCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagca                                                                                            >  2:500378/1‑150 (MQ=255)
                                              ggggCACAGAAGCAGCTGGCGGGCGCCGGGTTCAGCAACCTGGTGCTGCAACAGCCGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGATAAGGCCTGCTTTATATGGACAGCTTTGCGCGTGAAGGGGCTGTAAagca                                                                                            <  2:795419/150‑1 (MQ=255)
                                                ggCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagcagc                                                                                          >  1:791531/1‑150 (MQ=255)
                                                 gCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGcc                                                                                         >  1:1229289/1‑150 (MQ=255)
                                                     aGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGGCCGTGAAGGGGCTGTAAATCAGCCCCGc                                                                                     >  2:294123/1‑150 (MQ=255)
                                                     aGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGc                                                                                     >  1:1770601/1‑150 (MQ=255)
                                                     aGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGc                                                                                     >  2:1344536/1‑150 (MQ=255)
                                                     aGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACAAAGAAGGCCTTGCCGCAGCGAGGCGACAAGGCCTTCTTTATATGGACCGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGc                                                                                     >  1:606394/1‑150 (MQ=255)
                                                         gcagcTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCAcc                                                                                 <  1:606875/150‑1 (MQ=255)
                                                              tggctggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGa                                                                                                                                              >  1:1141107/1‑84 (MQ=255)
                                                              tggctggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGa                                                                                                                                              <  2:1141107/84‑1 (MQ=255)
                                                              tggctggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAg                                                                            >  1:1287087/1‑150 (MQ=255)
                                                                  tggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgcccc                                                                        >  1:1486299/1‑150 (MQ=255)
                                                                   ggggCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCt                                                                                                    >  1:849097/1‑121 (MQ=255)
                                                                   ggggCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCt                                                                                                    <  2:849097/121‑1 (MQ=255)
                                                                   ggggCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgccccg                                                                       >  1:1764996/1‑150 (MQ=255)
                                                                       ccGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCTTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgccccgcccc                                                                   <  2:499850/150‑1 (MQ=255)
                                                                            ttCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTgcttta                                                                                              <  2:721304/118‑6 (MQ=255)
                                                                            ttCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTgcttta                                                                                              >  1:721304/1‑113 (MQ=255)
                                                                                gCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCAtt                                                              >  1:1817996/1‑146 (MQ=255)
                                                                                gCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCAtt                                                              <  2:1817996/146‑1 (MQ=255)
                                                                                gCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCtt                                                          >  2:186889/1‑150 (MQ=255)
                                                                                  aaCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTc                                                        <  2:435691/150‑1 (MQ=255)
                                                                                  aaCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTc                                                        <  1:625935/150‑1 (MQ=255)
                                                                                   aCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCg                                                       <  1:500378/150‑1 (MQ=255)
                                                                                    ccTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTCCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCgg                                                      >  2:661492/1‑150 (MQ=255)
                                                                                       gTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGaca                                                   >  1:1220163/1‑150 (MQ=255)
                                                                                           ctgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTg                                               >  1:1035322/1‑150 (MQ=255)
                                                                                            tgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGGCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                              >  1:1266966/1‑150 (MQ=255)
                                                                                            tgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCATCTATTCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGTGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                              >  2:274104/1‑150 (MQ=255)
                                                                                             gcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGc                                             >  1:1101236/1‑150 (MQ=255)
                                                                                              cAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCt                                            >  2:1930340/1‑150 (MQ=255)
                                                                                                aCAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTTAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGt                                          >  2:1856330/1‑150 (MQ=255)
                                                                                                aCAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGt                                          <  2:1486299/150‑1 (MQ=255)
                                                                                                  aGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAATGCCTTCTTTATATGGACGGCATTGCCCGTCAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCAGACACCTGGCTGTTc                                        >  1:1106489/1‑150 (MQ=255)
                                                                                                     ggCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAAc                                     <  1:1344536/150‑1 (MQ=255)
                                                                                                        aGATCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATTCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGTGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                                  <  2:2225/150‑1 (MQ=255)
                                                                                                            ccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtata                                                                                                                              >  1:1547395/1‑54 (MQ=255)
                                                                                                            ccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtata                                                                                                                              <  2:1547395/54‑1 (MQ=255)
                                                                                                             cGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGaa                             <  1:729096/150‑1 (MQ=255)
                                                                                                              gTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAg                            <  2:1764996/150‑1 (MQ=255)
                                                                                                              gTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGCAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACCCCTGGGTGTTCAACGTCCAGAAg                            >  1:121767/1‑150 (MQ=255)
                                                                                                                 aGTAACGAAGAAGACCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCa                         >  2:1914844/1‑150 (MQ=255)
                                                                                                                   tAACGAATAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATa                       <  1:780173/150‑1 (MQ=255)
                                                                                                                    tataaagaagGCCTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGac                                                    <  2:994136/116‑1 (MQ=37)
                                                                                                                    tataaagaagGCCTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACAGTCAGCCCCGCCCCGCATTCCTTTCGGac                                                    >  1:994136/5‑120 (MQ=37)
                                                                                                                     aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACa                     >  2:1052755/1‑150 (MQ=255)
                                                                                                                        aagaagGCCTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCa                  >  2:1503595/1‑150 (MQ=39)
                                                                                                                          gaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCAcc                >  2:996363/1‑150 (MQ=255)
                                                                                                                               ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                         >  1:756013/1‑104 (MQ=38)
                                                                                                                               ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                         <  2:756013/104‑1 (MQ=38)
                                                                                                                               ccTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccac           >  2:819942/1‑150 (MQ=255)
                                                                                                                                cTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccacc          >  2:261129/1‑150 (MQ=255)
                                                                                                                                cTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccacc          >  2:1480901/1‑150 (MQ=255)
                                                                                                                                cTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccacc          >  2:1488881/1‑150 (MQ=255)
                                                                                                                                  tgtcGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAg        <  1:1503595/147‑1 (MQ=255)
                                                                                                                                    ccGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGa      <  2:1229289/150‑1 (MQ=255)
                                                                                                                                     cGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGaa     <  1:1930340/150‑1 (MQ=255)
                                                                                                                                       cAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccacca         <  1:785025/144‑1 (MQ=255)
                                                                                                                                       cAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccacca         >  2:785025/1‑144 (MQ=255)
                                                                                                                                        aGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGaacaa  <  1:996363/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGAACAA  >  NC_002947/5815455‑5815740

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: