Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,312:1 +G 6.9% intergenic (+21/+23) ydeJ → / ← dcp inactive PncC family protein/dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,625,3121.G 6.9% 404.6 / 24.5 101intergenic (+21/+23)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base . (45/49);  new base G (3/4);  total (48/53)
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.15e-01

Only 100 of 103 total aligned reads displayed.

GTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAAC  >  NC_000913/1625163‑1625456
                                                                                                                                                      |                                                                                                                                                 
gTCTGGTTTGCGTGGCATATTAAAGGCCAGACCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCATTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATAtc                                                                                                                                                    >  1:1448789/1‑150 (MQ=255)
gTCTGGTTTGCGGGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATAtc                                                                                                                                                    <  2:1207007/150‑1 (MQ=255)
   tGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTCGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑Ca                                                                                                                                               <  1:1311822/150‑1 (MQ=255)
    ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                              >  2:1019796/1‑150 (MQ=255)
    ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                              >  2:1165453/1‑150 (MQ=255)
    ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                              >  1:193864/1‑150 (MQ=255)
      tttGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGgc                                                                                                                                            <  2:810613/150‑1 (MQ=255)
         gCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGACCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCT‑GAGATATCGc‑caggccag                                                                                                                                         <  1:122446/150‑1 (MQ=255)
           gtttcATATATAAGGCCAGATCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGATACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCT‑GAGATATCGC‑CAGgccagcc                                                                                                                                       <  1:146488/146‑1 (MQ=255)
           gTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGATACTGTATTATCTTTAGCGGTGAGGTTTGCACTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGgccagcc                                                                                                                                       <  2:857120/150‑1 (MQ=255)
           gTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCT‑GAGATATCGC‑CAGgccagcc                                                                                                                                       <  1:316056/150‑1 (MQ=255)
            tGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGgccagcca                                                                                                                                      >  1:1337006/1‑150 (MQ=255)
            tGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCT‑GAGATATCGC‑CAGgccagcct                                                                                                                                      <  1:1634317/150‑2 (MQ=255)
             ggCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCAt                                                                                                                                     >  1:201307/1‑150 (MQ=255)
                atatTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATtgg                                                                                                                                  >  1:1525464/1‑150 (MQ=255)
                atatTAAAGGCCAGAACTACACTGCGGGTAGGCATTTTGCTGCCGACTGCGAAACGGTCTTAGCTTTAGCGGCGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATtgg                                                                                                                                  <  2:1638979/150‑1 (MQ=255)
                  atTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATtggtg                                                                                                                                <  2:1182006/150‑1 (MQ=255)
                   ttAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATtggtgg                                                                                                                               <  2:1064100/150‑1 (MQ=255)
                     aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGTCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCGGCTGCAATTACTGCTATATCACGGATGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGgtggtg                                                                                                                             >  2:73358/1‑150 (MQ=255)
                     aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGTCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTggtggtt                                                                                                                             >  2:1579414/1‑149 (MQ=255)
                     aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGgtggtg                                                                                                                             >  2:1411348/1‑150 (MQ=255)
                      aaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGTGTTTTTTTTTTATCTCGCCGAGCTGATGGATTTATTGATATAAAATGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGgtgt                                                                                                                            <  2:445778/150‑1 (MQ=255)
                       aGGCCAGAACTACACTGCGGTTATGCATGTTGCTGGCGACTGCGAAACGGTATTAGCGTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCGGCAATTACTGCTATAACCAGGCGGGCCTGGCGATATC‑T‑CCGGCCAGCCATTGGTGGTGtt                                                                                                                           <  2:1624408/150‑1 (MQ=255)
                        ggCCAGATCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGttt                                                                                                                          <  2:1162863/150‑1 (MQ=255)
                         gCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTta                                                                                                                         >  2:1546687/1‑150 (MQ=255)
                          ccAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGGGTTtat                                                                                                                        >  2:847168/1‑150 (MQ=255)
                             gatcTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATg                                                                                                                     <  1:958219/147‑1 (MQ=255)
                                cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGtta                                                                                                                  >  1:996121/1‑149 (MQ=255)
                                cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                  <  2:1514993/150‑1 (MQ=255)
                                cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                  >  2:504577/1‑150 (MQ=255)
                                 tACACTTCGGTTATGCATTTTGCGGGCGACTGCGCAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCACAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCa                                                                                                                 <  2:917933/150‑1 (MQ=255)
                                  acacTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGGGTTTATATGTTCaa                                                                                                                >  1:830776/1‑150 (MQ=255)
                                     cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                             <  1:945039/150‑1 (MQ=255)
                                     cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                             >  1:170404/1‑150 (MQ=255)
                                        cGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGGGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACg                                                                                                          <  1:64208/150‑1 (MQ=255)
                                               gCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgc                                                                                                   >  2:1612769/1‑150 (MQ=255)
                                               gCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgc                                                                                                   >  1:668348/1‑150 (MQ=255)
                                                 aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgcag                                                                                                 <  1:1451628/150‑1 (MQ=255)
                                                      gCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                            >  2:862442/1‑150 (MQ=255)
                                                      gCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                            >  1:1304183/1‑150 (MQ=255)
                                                      gCTGGCGACTGCGAAACGGTATTAGCTTTAACGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                            >  1:659051/1‑150 (MQ=255)
                                                        tGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGc                                                                                          <  1:1655719/150‑1 (MQ=255)
                                                           cGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATa                                                                                       >  2:963569/1‑150 (MQ=255)
                                                           cGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGCTGCAACATCgtggct                                                                                       <  1:339349/150‑6 (MQ=255)
                                                                gCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCtt                                                                                  >  2:1308662/1‑150 (MQ=255)
                                                                 cGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTa                                                                                 <  2:1607677/150‑1 (MQ=255)
                                                                     aCGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTg                                                                             <  2:1337006/150‑1 (MQ=255)
                                                                      cGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGc                                                                            >  1:1075031/1‑150 (MQ=255)
                                                                       ggTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGcc                                                                           <  1:238507/150‑1 (MQ=255)
                                                                       ggTATTAGCTTTAGCGGTGAGGTTTGCCCGCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGcc                                                                           <  2:1066760/150‑1 (MQ=255)
                                                                       ggTATTAGCTTTAGCGGTGAGGCTTGCCCTCGCCCAGCGGCTGCAATTACTGCAATACCCAGGCCGGCCGGGCGATAGC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTAAGGTGcc                                                                           <  2:229211/150‑1 (MQ=255)
                                                                        gTATTAGCTTTAGCGGTGAGGGTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCt                                                                          <  2:831070/150‑1 (MQ=255)
                                                                             aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCTTTGGTGGTGTTTATATGTTCTAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                     >  2:684378/1‑150 (MQ=255)
                                                                             aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                     >  2:1155928/1‑150 (MQ=255)
                                                                             aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                     <  2:326827/150‑1 (MQ=255)
                                                                                 ttAGCGGGGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGcc                                                                 <  2:387855/150‑1 (MQ=255)
                                                                                   aGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCAt                                                               <  1:1633692/150‑1 (MQ=255)
                                                                                      ggTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CTGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATATTCTTAGGTGCCTTACCGCGCCATTGt                                                            >  1:689734/1‑150 (MQ=255)
                                                                                      ggTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACAGCGCCATTGt                                                            >  1:618006/1‑150 (MQ=255)
                                                                                       gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTTCTATAACCAGGCTGGCCTGGAGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                           <  2:467617/150‑1 (MQ=255)
                                                                                       gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCGATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                           <  1:208433/150‑1 (MQ=255)
                                                                                       gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                           <  2:337181/150‑1 (MQ=255)
                                                                                       gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                           >  1:1651094/1‑150 (MQ=255)
                                                                                       gTGAGGTTAGCCCTCGCACAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGAAATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTAACCGCGCCATTGTc                                                           <  1:684348/150‑1 (MQ=255)
                                                                                        tGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCg                                                          <  1:567743/150‑1 (MQ=255)
                                                                                          aGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGAt                                                        <  1:640772/150‑1 (MQ=255)
                                                                                           ggTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATa                                                       <  1:924670/150‑1 (MQ=255)
                                                                                                gCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGc                                                  <  1:862442/150‑1 (MQ=255)
                                                                                                   cTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGtt                                               >  1:1522774/1‑150 (MQ=255)
                                                                                                      gCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCa                                            >  1:1143700/1‑150 (MQ=255)
                                                                                                       cccAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTTGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAg                                           <  1:837399/150‑1 (MQ=255)
                                                                                                           gctgctGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCt                                       <  1:1292773/150‑1 (MQ=255)
                                                                                                              gctgcAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCg                                    <  1:1579414/150‑1 (MQ=255)
                                                                                                              gctgcAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCg                                    <  1:1165453/150‑1 (MQ=255)
                                                                                                               ctgcAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGATCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGc                                   <  2:1075031/150‑1 (MQ=255)
                                                                                                                  cAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGt                                <  2:1562023/150‑1 (MQ=255)
                                                                                                                     ttACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTAc                             >  1:1428738/1‑150 (MQ=255)
                                                                                                                       aCTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACct                           <  1:552314/150‑1 (MQ=255)
                                                                                                                            tataACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                      >  1:964008/1‑150 (MQ=255)
                                                                                                                            tataACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                      >  1:1666893/1‑150 (MQ=255)
                                                                                                                            tataACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                      >  1:1673352/1‑150 (MQ=255)
                                                                                                                            tataACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATCTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCTCTGTTACCTCTGGa                      >  2:741243/1‑150 (MQ=255)
                                                                                                                              taACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTTGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGaaa                    <  1:1455404/150‑1 (MQ=255)
                                                                                                                              taACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGaaa                    <  1:1474485/150‑1 (MQ=255)
                                                                                                                                aCCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAATCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAgg                  <  2:907637/150‑1 (MQ=255)
                                                                                                                                 aatggCTGGCCT‑GAGATATCGC‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGa                 >  2:242814/4‑150 (MQ=255)
                                                                                                                                 aatggCTGGCCT‑GAGATATCGC‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGa                 >  2:756524/4‑150 (MQ=255)
                                                                                                                                  cAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAAATCGTCGCTGTTACCTCTGGAAAGGAt                >  2:668967/1‑150 (MQ=255)
                                                                                                                                   aGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc               <  2:668348/150‑1 (MQ=255)
                                                                                                                                   aGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc               <  1:1546687/150‑1 (MQ=255)
                                                                                                                                    ggCTGGCCT‑GAGATATCGC‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCg              >  2:145453/1‑150 (MQ=255)
                                                                                                                                      ctggcctggcGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGcc            <  1:937437/150‑1 (MQ=255)
                                                                                                                                      ctggcctggcGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGcc            >  1:1414415/1‑150 (MQ=255)
                                                                                                                                      ctggcct‑gagATATC‑GCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCGTCGGCGTAATCTTAGGTGCCTTACCGCTCCATTGTCGATACAGGCGTTCCTTATCGGCGCTGTTTGCTCTGGAAAGGATAGcc            >  2:297071/1‑150 (MQ=255)
                                                                                                                                       tggcctggcGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCt           >  2:442071/1‑150 (MQ=255)
                                                                                                                                         gcctggcGATATC‑T‑CAGGACAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcg         <  2:1178093/150‑1 (MQ=255)
                                                                                                                                          cctggcGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcgc        >  1:998738/1‑150 (MQ=255)
                                                                                                                                             ggcGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGaa     >  1:1504184/1‑150 (MQ=255)
                                                                                                                                               cGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGaaaa   >  1:1380622/1‑150 (MQ=255)
                                                                                                                                                gATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAAc  >  2:21567/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAAC  >  NC_000913/1625163‑1625456

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