Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,270,602 T→C 10.5% 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,6020TC10.5% 549.8 / 44.4 190intergenic (+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 (84/86);  new base C (10/10);  total (94/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-02

Only 100 of 194 total aligned reads displayed.

CCACCAGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTG  >  NC_000913/2270454‑2270750
                                                                                                                                                    |                                                                                                                                                    
  accaGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTa                                                                                                                                                    <  1:177468/149‑1 (MQ=255)
    cagcagTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAgg                                                                                                                                                  <  3:206535/149‑1 (MQ=255)
         gtggtgTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTg                                                                                                                                             <  3:1287806/149‑1 (MQ=255)
            gtgtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtt                                                                                                                                          >  2:404379/1‑149 (MQ=255)
            gtgtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtt                                                                                                                                          >  3:1460086/1‑149 (MQ=255)
             tgtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCttt                                                                                                                                         <  3:1055964/149‑1 (MQ=255)
             tgtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCttt                                                                                                                                         >  4:1542129/1‑149 (MQ=255)
             tgtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCCCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCttt                                                                                                                                         <  4:887226/149‑1 (MQ=255)
              gtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttc                                                                                                                                        <  4:1366785/149‑2 (MQ=255)
               ttattaTTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCa                                                                                                                                                <  2:1541458/140‑1 (MQ=255)
               ttattaTTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCa                                                                                                                                                >  1:1541458/1‑140 (MQ=255)
               ttattaTTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCttttt                                                                                                                                       >  4:126564/1‑149 (MQ=255)
                tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                      <  3:121274/149‑1 (MQ=255)
                tattatTGCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                      <  1:325538/149‑1 (MQ=255)
                tatcatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCtttttt                                                                                                                                      <  3:1184446/149‑1 (MQ=255)
                  ttattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGt                                                                                                                                    >  3:1398079/1‑149 (MQ=255)
                  ttattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttca                                                                                                                                       >  4:401751/1‑144 (MQ=255)
                  ttattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttca                                                                                                                                       <  3:401751/146‑3 (MQ=255)
                         cAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTttcatt                                                                                                                             <  4:144841/149‑1 (MQ=255)
                              tgaatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAg                                                                                                                                                   <  4:1424383/122‑1 (MQ=255)
                              tgaatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                        >  1:1367701/1‑149 (MQ=255)
                              tgaatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                        <  3:773990/149‑1 (MQ=255)
                                aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCaga                                                                                                                      <  1:507895/149‑1 (MQ=255)
                                aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATACAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCaga                                                                                                                      >  2:801702/1‑149 (MQ=255)
                                 atgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGt                                                                                                                                    <  4:1202028/134‑1 (MQ=255)
                                 atgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCagag                                                                                                                     >  3:230125/1‑149 (MQ=255)
                                 atgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCagag                                                                                                                     >  4:1133030/1‑149 (MQ=255)
                                  tgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAgaga                                                                                                                    >  4:74744/1‑149 (MQ=255)
                                       cGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatg                                                                                                               >  1:1035402/1‑149 (MQ=255)
                                        gTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatga                                                                                                              <  2:682443/149‑1 (MQ=255)
                                            tGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgatggt                                                                                                          >  1:562960/1‑147 (MQ=255)
                                            tGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGc                                                                                                                                            >  2:1130827/1‑115 (MQ=255)
                                            tGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGc                                                                                                                                            <  1:1130827/115‑1 (MQ=255)
                                               aagaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGACGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTcc                                                                                                       >  2:1060066/1‑149 (MQ=255)
                                                 gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                        <  3:487419/130‑1 (MQ=255)
                                                 gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATtcatca                                                                                                                        >  4:487419/1‑130 (MQ=255)
                                                 gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCcgc                                                                                                     <  3:1498999/149‑1 (MQ=255)
                                                    gCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTc                                                                                                  >  1:734845/1‑149 (MQ=255)
                                                       ttACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTg                                                                                               >  2:372480/1‑149 (MQ=255)
                                                          cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATc                                                                                            <  1:285878/149‑1 (MQ=255)
                                                          cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATc                                                                                            <  3:819719/149‑1 (MQ=255)
                                                           aaCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATcc                                                                                           >  4:1489408/1‑149 (MQ=255)
                                                           aaCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATcc                                                                                           >  3:666887/1‑149 (MQ=255)
                                                           aaCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttcagtt                                                                                                                                    >  3:845883/1‑103 (MQ=39)
                                                           aaCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttcagtt                                                                                                                                    <  4:845883/108‑6 (MQ=39)
                                                            aCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCa                                                                                          >  4:1393213/1‑149 (MQ=255)
                                                             cGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAg                                                                                         >  2:128284/1‑149 (MQ=255)
                                                              gAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAgg                                                                                        <  3:982320/149‑1 (MQ=255)
                                                               aaGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGAt                                                                                             >  2:1244581/1‑143 (MQ=255)
                                                               aaGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGc                                                                                       <  1:1034889/149‑1 (MQ=255)
                                                                aGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCt                                                                                      >  1:338202/1‑149 (MQ=255)
                                                                     gCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGt                                                                                 >  2:380110/1‑149 (MQ=255)
                                                                        gggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCg                                                                              <  2:420573/149‑1 (MQ=255)
                                                                         ggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCgg                                                                             <  1:376574/149‑1 (MQ=255)
                                                                              tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCAtt                                                                        >  2:21877/1‑149 (MQ=255)
                                                                              tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCAtt                                                                        >  4:580395/1‑149 (MQ=255)
                                                                              tcAGCCGAAGCTAATAATCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATACATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCAtt                                                                        <  2:739149/149‑1 (MQ=255)
                                                                                 gCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGtt                                                                     <  1:877946/149‑1 (MQ=255)
                                                                                    gCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATc                                                                  <  1:1267529/149‑1 (MQ=255)
                                                                                         taataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTa                                                             <  3:746914/149‑1 (MQ=255)
                                                                                         taataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTa                                                             >  3:114427/1‑149 (MQ=255)
                                                                                            taaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAt                                                          <  4:666887/149‑1 (MQ=255)
                                                                                               aCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgtt                                                       <  4:1526947/149‑1 (MQ=255)
                                                                                                 cGTTTGGATGCAATCCCTTGGCTATACTGACGAGTTCACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGATAAAATGATGGTTCCGCGTCTTGATCTATGCTATAGTTCGGGCATTGTTATCTTTTAATtgtggt                                                     >  2:1319330/1‑149 (MQ=255)
                                                                                                  gTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                    <  3:1248490/149‑1 (MQ=255)
                                                                                                  gTTTGGATGCAAGCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTGCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                    >  4:1328402/1‑149 (MQ=255)
                                                                                                   tttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                    <  3:183944/148‑1 (MQ=255)
                                                                                                   tttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                    >  4:183944/1‑148 (MQ=255)
                                                                                                    ttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGt                                                  >  1:1330881/1‑149 (MQ=255)
                                                                                                       gATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgtt                                                       >  1:1081431/1‑141 (MQ=255)
                                                                                                       gATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgtt                                                       <  2:1081431/141‑1 (MQ=255)
                                                                                                       gATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAAt                                               >  4:1015451/1‑149 (MQ=255)
                                                                                                            aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                                                                             <  2:538547/82‑1 (MQ=255)
                                                                                                            aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTc                                                                                >  4:481870/1‑111 (MQ=255)
                                                                                                            aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTc                                                                                <  3:481870/111‑1 (MQ=255)
                                                                                                            aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCg                                                                              <  4:624641/113‑1 (MQ=255)
                                                                                                            aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCg                                                                              >  3:624641/1‑113 (MQ=255)
                                                                                                                  ttGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAAtt                                    >  2:1030701/1‑149 (MQ=255)
                                                                                                                     gCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAAc                                 <  1:372480/149‑1 (MQ=255)
                                                                                                                           cTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGa                              >  4:1237200/1‑146 (MQ=255)
                                                                                                                            tGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCa                          >  3:539847/1‑149 (MQ=255)
                                                                                                                                gAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGtt                      >  2:298749/1‑149 (MQ=255)
                                                                                                                                   ttAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCAt                   >  1:1405587/1‑149 (MQ=255)
                                                                                                                                   ttAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCAt                   >  1:121802/1‑149 (MQ=255)
                                                                                                                                   aaacaaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCAt                   >  1:822510/8‑149 (MQ=255)
                                                                                                                                    tAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATa                  <  1:404379/149‑1 (MQ=255)
                                                                                                                                         gAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGtt                      <  2:922532/140‑1 (MQ=255)
                                                                                                                                         gAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGtt                      >  1:922532/1‑140 (MQ=255)
                                                                                                                                         gAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACgcgc             >  4:813353/1‑149 (MQ=255)
                                                                                                                                          aaaGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCt                                                                                                 <  1:450358/64‑1 (MQ=255)
                                                                                                                                          aaaGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATc                                                                  <  1:197015/95‑1 (MQ=255)
                                                                                                                                          aaaGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCt            >  3:1159824/1‑149 (MQ=255)
                                                                                                                                          aaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCt            >  4:289944/1‑149 (MQ=255)
                                                                                                                                           aaGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTc           <  4:1460086/149‑1 (MQ=255)
                                                                                                                                               aCTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCaa       <  2:1465888/149‑1 (MQ=255)
                                                                                                                                               aCTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCaa       >  3:855944/1‑149 (MQ=255)
                                                                                                                                                 tGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAtt     >  1:644150/1‑149 (MQ=255)
                                                                                                                                                  gCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAttt    >  1:1491879/1‑149 (MQ=255)
                                                                                                                                                    ttAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTg  >  2:652925/1‑149 (MQ=255)
                                                                                                                                                    ttAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTg  <  3:1542129/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCACCAGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTG  >  NC_000913/2270454‑2270750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: