Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,270,602 T→C 16.8% intergenic (+57/‑124) mepS → / → rtn murein DD‑endopeptidase, space‑maker hydrolase, mutational suppressor of prc thermosensitivity, outer membrane lipoprotein, weak murein LD‑carboxypeptidase/resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,270,6020TC16.8% 210.3 / 29.1 93intergenic (+57/‑124)mepS/rtnmurein DD‑endopeptidase, space‑maker hydrolase, mutational suppressor of prc thermosensitivity, outer membrane lipoprotein, weak murein LD‑carboxypeptidase/resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (38/39);  new base C (7/9);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCG  >  NC_000913/2270503‑2270701
                                                                                                   |                                                                                                   
gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCt                                                                                                     <  1:933198/100‑1 (MQ=255)
gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCt                                                                                                     >  1:37691/1‑100 (MQ=255)
    cGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAgg                                                                                                 <  2:37691/100‑1 (MQ=255)
     gTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGc                                                                                                <  1:397175/100‑1 (MQ=255)
     gTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGc                                                                                                <  2:827219/100‑1 (MQ=255)
         cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTg                                                                                            >  2:2207183/1‑100 (MQ=255)
         cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTg                                                                                            >  1:189928/1‑100 (MQ=255)
         cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTg                                                                                            <  2:894994/100‑1 (MQ=255)
           aCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCt                                                                                          <  2:2103138/100‑1 (MQ=255)
           aCGAAGCACGCCGGGTTCTCAGCCGCAGCTAAAAAACCGATTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCt                                                                                          <  2:1828326/100‑1 (MQ=255)
           aCGAAGCACGCCGGGGTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGGTAACTGAAAGCACTGCTTAGGCAGTGCt                                                                                          >  2:851019/1‑100 (MQ=255)
              aaGCAGGCCGGGTTCTCAGCCGCACGTAATAAACCGGTTGGATGCAATGTCTTGGCTAGCCTGACTAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttc                                                                                       <  2:1050395/100‑2 (MQ=255)
              aaGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttc                                                                                       <  1:490151/100‑2 (MQ=255)
                gCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttcag                                                                                     <  2:1958053/100‑4 (MQ=39)
                gCACGCCGGGTTATCAGCTGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGTTATCCTGACGAGTTAACTTAAAGCACTGCCTAAGCAGTGCtttcag                                                                                     <  1:1276195/100‑4 (MQ=37)
                gCACGCCGGGTGCTCAGCCGCAGCTAATGAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttcag                                                                                     <  2:2016627/100‑4 (MQ=38)
                 cACGCCGGGTTCTCAGCCGCAGCTAATAACCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                    <  2:1556394/100‑1 (MQ=255)
                 cACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                    <  2:255820/100‑1 (MQ=255)
                 cACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttcagt                                                                                    <  1:11705/100‑5 (MQ=39)
                  aCGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGACAGCACTGCTTAGGCAGTGCTTTTTTGt                                                                                   >  2:2458034/1‑100 (MQ=255)
                   cGCCGGGTTCTCAGCCGCAGCAAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtt                                                                                         <  2:305950/93‑1 (MQ=255)
                     ccGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGtttt                                                                                <  2:354078/100‑1 (MQ=255)
                            ctcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTtcattcat                                                                         >  1:416130/1‑100 (MQ=255)
                             tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatc                                                                        <  1:2185915/100‑1 (MQ=255)
                             tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatc                                                                        <  2:1973889/100‑1 (MQ=255)
                             tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatc                                                                        <  2:1915455/100‑1 (MQ=255)
                               aGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatc                                                                        >  2:485938/1‑98 (MQ=255)
                               aGCCGCAGCTAATAAACAGTTTGGATGCAATCCCTCGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGTCAGTGCTTTTTTGTTTTCATTCATCag                                                                      >  2:487653/1‑100 (MQ=255)
                                gCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCaga                                                                     >  2:2068836/1‑100 (MQ=255)
                                     aGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGaaaa                                                                >  1:1683319/1‑100 (MQ=255)
                                      gCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                               <  2:1437020/100‑1 (MQ=255)
                                      gCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                               <  2:1435307/100‑1 (MQ=255)
                                      gCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                               <  1:1057127/100‑1 (MQ=255)
                                      gCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                               <  1:723312/100‑1 (MQ=255)
                                       cTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAGatg                                                              >  2:1463497/1‑100 (MQ=255)
                                       cTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatg                                                              >  2:919705/1‑100 (MQ=255)
                                         aataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                            <  2:482410/100‑1 (MQ=255)
                                         aataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                            <  2:2112282/100‑1 (MQ=255)
                                           taaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGt                                                          >  1:1166793/1‑100 (MQ=255)
                                               ccGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTcc                                                      <  1:483940/100‑1 (MQ=255)
                                                cGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCcg                                                     <  1:956232/100‑1 (MQ=255)
                                                cGTTTGCATGCAATCACTTGGCTAACCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTATTTTTTTTCATCACAAGAAATTATGttttag                                                     >  2:281140/1‑97 (MQ=255)
                                                  tttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCcgcg                                                   <  2:2281808/100‑1 (MQ=255)
                                                  tttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCcgcg                                                   >  1:2020937/1‑100 (MQ=255)
                                                    tGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGt                                                  >  2:2024733/1‑99 (MQ=255)
                                                    tGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTc                                                 >  2:2500098/1‑100 (MQ=255)
                                                    tGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTc                                                 >  1:2330336/1‑100 (MQ=255)
                                                     ggATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCt                                                <  1:2437543/100‑1 (MQ=255)
                                                      gATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCtt                                               <  2:1039986/100‑1 (MQ=255)
                                                      gATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCtt                                               >  2:2014393/1‑100 (MQ=255)
                                                       aTGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTg                                              >  2:1247030/1‑100 (MQ=255)
                                                       aTGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTg                                              >  1:1095982/1‑100 (MQ=255)
                                                        tGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGa                                             >  1:721841/1‑100 (MQ=255)
                                                        tGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGa                                             >  2:117936/1‑100 (MQ=255)
                                                          cAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATc                                           >  1:698227/1‑100 (MQ=255)
                                                              cccTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAgg                                       <  2:352591/100‑1 (MQ=255)
                                                               ccTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGc                                      <  2:2000974/100‑1 (MQ=255)
                                                                   ggCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCtata                                  >  1:2049443/1‑100 (MQ=255)
                                                                     cTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTTTTTTTTTTTTTCAGAGAAAAGGATGTTTCCGCGTTTTGATCTATGCtgttg                                 >  2:1911487/1‑95 (MQ=255)
                                                                       aTCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTcc                              <  2:1512370/100‑1 (MQ=255)
                                                                       aTCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTcc                              >  1:166857/1‑100 (MQ=255)
                                                                          cTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCgg                            >  2:803655/1‑99 (MQ=255)
                                                                           tGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTc                          <  1:324564/100‑1 (MQ=255)
                                                                            gACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCa                         >  1:2244653/1‑100 (MQ=255)
                                                                             aCGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCAt                        <  1:1310276/100‑1 (MQ=255)
                                                                               gAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTg                      <  1:1091751/100‑1 (MQ=255)
                                                                               gAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTg                      >  1:634907/1‑100 (MQ=255)
                                                                                aGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGt                     <  2:1273453/100‑1 (MQ=255)
                                                                                  ttAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTa                   <  2:320121/100‑1 (MQ=255)
                                                                                  ttAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTa                   <  1:472683/100‑1 (MQ=255)
                                                                                  ttAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTa                   <  2:1095977/100‑1 (MQ=255)
                                                                                  ttAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTa                   <  1:658981/100‑1 (MQ=255)
                                                                                    acaaaaaaGCACTGCCTCAGCCGTGCTTTTTTGTTTTCATTCATCAGAGCAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATc                 >  2:2591827/6‑100 (MQ=255)
                                                                                    aaCTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATc                 <  2:2259765/100‑1 (MQ=255)
                                                                                     caaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCTAGGCTATAGTCCGGTCATTGTTATCt                >  2:946639/5‑100 (MQ=38)
                                                                                     caaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCt                >  2:1897532/5‑100 (MQ=39)
                                                                                     caaaaaaGCACTGCCTAAGCAGTGCTTTCTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCt                >  1:160938/5‑100 (MQ=38)
                                                                                     aCTGAAAGCACTGCTTAGGCAGTGCTTTTTTTTTTTTTTTTTTCAGAGAAAATGATGTTTCCGCGTCTTGTTTTAGGCTATAGTCTGGTCATTGtttttt                >  2:2511612/1‑96 (MQ=255)
                                                                                     aCTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCt                >  1:354813/1‑100 (MQ=255)
                                                                                      cTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCtt               >  1:2492788/1‑100 (MQ=255)
                                                                                      aaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCtt               >  1:2174475/4‑100 (MQ=39)
                                                                                       aaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCttt              >  1:787454/3‑100 (MQ=255)
                                                                                        gAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCtttt             <  1:2458047/100‑1 (MQ=255)
                                                                                          aaGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTaa           >  1:2371310/1‑100 (MQ=255)
                                                                                           aGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTaaa          <  1:520535/100‑1 (MQ=255)
                                                                                            gCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAt         >  1:1820614/1‑100 (MQ=255)
                                                                                             cACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtg        <  2:2140821/100‑1 (MQ=255)
                                                                                              aCTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgt       >  1:1915598/1‑100 (MQ=255)
                                                                                              aCTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgt       >  1:613077/1‑100 (MQ=255)
                                                                                                tGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgttg     <  1:2387688/100‑1 (MQ=255)
                                                                                                tGCTTAGGCAGTGCTTTTTTGTTTTCATTCATAAGAGAAAATGCTGTTTCCGAGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgttg     >  2:577413/1‑100 (MQ=255)
                                                                                                 gCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgttgt    <  2:613074/100‑1 (MQ=255)
                                                                                                 gCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgttgt    >  1:2052122/1‑100 (MQ=255)
                                                                                                 gCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgttgt    <  2:1915589/100‑1 (MQ=255)
                                                                                                  cTTAGGCAGTGCTTTTTTTTTTTTTTTTATAAGAGAAACTGATGTTTCCGCGTCTTGATTCAGTCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc   >  2:1600980/1‑100 (MQ=255)
                                                                                                   ttAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCg  >  1:2514731/1‑100 (MQ=255)
                                                                                                   ttAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCg  >  2:2077533/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCG  >  NC_000913/2270503‑2270701

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: