Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,304 Δ1 bp 7.3% intergenic (+13/+31) 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,3040G. 7.3% 298.9 / 17.1 96intergenic (+13/+31)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base G (40/49);  new base . (3/4);  total (43/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 = 5.10e-01

Only 100 of 105 total aligned reads displayed.

CGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAA  >  NC_000913/1625155‑1625453
                                                                                                                                                     |                                                                                                                                                       
cGGGTACCGGCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGctggcctg                                                                                                                                                         <  1:356729/150‑1 (MQ=255)
   gTACCGTCTGGTTTGCGTGGCATATTCAAGGCCACACCTCCACTGCGGTTATGCAGTTTGCTGGCGACTGCGAACCGGTAGTAGCTTACGCGGTGAGGCTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCg                                                                                                                                                      <  2:349109/150‑1 (MQ=255)
   gTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGctggcctgaga                                                                                                                                                      <  2:193505/150‑4 (MQ=255)
       cGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGatat                                                                                                                                                  >  1:145747/1‑150 (MQ=255)
        gTCTGGTTTGCGTGGCATATTAAAGGCCAGACCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCATTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATAtc                                                                                                                                                 >  1:1448789/1‑150 (MQ=255)
        gTCTGGTTTGCGGGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATAtc                                                                                                                                                 <  2:1207007/150‑1 (MQ=255)
           tGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTCGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCa                                                                                                                                            <  1:1311822/150‑1 (MQ=255)
            ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                           >  1:193864/1‑150 (MQ=255)
            ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                           >  2:1019796/1‑150 (MQ=255)
            ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                           >  2:1165453/1‑150 (MQ=255)
              tttGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgc                                                                                                                                         <  2:810613/150‑1 (MQ=255)
                 gCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGACCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGccaggccag                                                                                                                                      <  1:122446/150‑1 (MQ=255)
                   gtttcATATATAAGGCCAGATCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGATACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcc                                                                                                                                    <  1:146488/146‑1 (MQ=255)
                   gTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGATACTGTATTATCTTTAGCGGTGAGGTTTGCACTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgccagcc                                                                                                                                    <  2:857120/150‑1 (MQ=255)
                   gTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcc                                                                                                                                    <  1:316056/150‑1 (MQ=255)
                    tGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGgccagcca                                                                                                                                   >  1:1337006/1‑150 (MQ=255)
                    tGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGAGATATCGCCAGgccagcct                                                                                                                                   <  1:1634317/150‑2 (MQ=255)
                     ggCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCAt                                                                                                                                  >  1:201307/1‑150 (MQ=255)
                        atatTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtgg                                                                                                                               >  1:1525464/1‑150 (MQ=255)
                        atatTAAAGGCCAGAACTACACTGCGGGTAGGCATTTTGCTGCCGACTGCGAAACGGTCTTAGCTTTAGCGGCGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtgg                                                                                                                               <  2:1638979/150‑1 (MQ=255)
                          atTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtg                                                                                                                             <  2:1182006/150‑1 (MQ=255)
                           ttAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATtggtgg                                                                                                                            <  2:1064100/150‑1 (MQ=255)
                             aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGTCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCGGCTGCAATTACTGCTATATCACGGATGGCCTGGCGATATCTCAGGCCAGCCATTGgtggtg                                                                                                                          >  2:73358/1‑150 (MQ=255)
                             aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGTCTGGCCTGGCGATATCTCAGGCCAGCCATTggtggtt                                                                                                                          >  2:1579414/1‑149 (MQ=255)
                             aaaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGgtggtg                                                                                                                          >  2:1411348/1‑150 (MQ=255)
                              aaGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGTGTTTTTTTTTTATCTCGCCGAGCTGATGGATTTATTGATATAAAATGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGgtgt                                                                                                                         <  2:445778/150‑1 (MQ=255)
                               aGGCCAGAACTACACTGCGGTTATGCATGTTGCTGGCGACTGCGAAACGGTATTAGCGTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCGGCAATTACTGCTATAACCAGGCGGGCCTGGCGATATCTCCGGCCAGCCATTGGTGGTGtt                                                                                                                        <  2:1624408/150‑1 (MQ=255)
                                ggCCAGATCTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGttt                                                                                                                       <  2:1162863/150‑1 (MQ=255)
                                 gCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTta                                                                                                                      >  2:1546687/1‑150 (MQ=255)
                                  ccAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGGGTTtat                                                                                                                     >  2:847168/1‑150 (MQ=255)
                                     gatcTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATg                                                                                                                  <  1:958219/147‑1 (MQ=255)
                                        cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGtta                                                                                                               >  1:996121/1‑149 (MQ=255)
                                        cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                               <  2:1514993/150‑1 (MQ=255)
                                        cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                               >  2:504577/1‑150 (MQ=255)
                                         tACACTTCGGTTATGCATTTTGCGGGCGACTGCGCAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCACAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCa                                                                                                              <  2:917933/150‑1 (MQ=255)
                                          acacTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGGGTTTATATGTTCaa                                                                                                             >  1:830776/1‑150 (MQ=255)
                                             cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                          >  1:170404/1‑150 (MQ=255)
                                             cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                          <  1:945039/150‑1 (MQ=255)
                                                cGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGGGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACg                                                                                                       <  1:64208/150‑1 (MQ=255)
                                                       gCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgc                                                                                                >  1:668348/1‑150 (MQ=255)
                                                       gCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgc                                                                                                >  2:1612769/1‑150 (MQ=255)
                                                         aTTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgcag                                                                                              <  1:1451628/150‑1 (MQ=255)
                                                              gCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                         >  2:862442/1‑150 (MQ=255)
                                                              gCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                         >  1:1304183/1‑150 (MQ=255)
                                                                tGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGc                                                                                       <  1:1655719/150‑1 (MQ=255)
                                                                   cGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATa                                                                                    >  2:963569/1‑150 (MQ=255)
                                                                   cGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGCTGCAACATCgtggct                                                                                    <  1:339349/150‑6 (MQ=255)
                                                                        gCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCtt                                                                               >  2:1308662/1‑150 (MQ=255)
                                                                         cGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTa                                                                              >  1:512552/1‑150 (MQ=255)
                                                                         cGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTa                                                                              <  2:1607677/150‑1 (MQ=255)
                                                                             aCGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTg                                                                          <  2:1337006/150‑1 (MQ=255)
                                                                              cGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGc                                                                         >  1:1075031/1‑150 (MQ=255)
                                                                               ggTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGcc                                                                        <  1:238507/150‑1 (MQ=255)
                                                                               ggTATTAGCTTTAGCGGTGAGGTTTGCCCGCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGcc                                                                        <  2:1066760/150‑1 (MQ=255)
                                                                               ggTATTAGCTTTAGCGGTGAGGCTTGCCCTCGCCCAGCGGCTGCAATTACTGCAATACCCAGGCCGGCCGGGCGATAGCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTAAGGTGcc                                                                        <  2:229211/150‑1 (MQ=255)
                                                                                gTATTAGCTTTAGCGGTGAGGGTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCt                                                                       <  2:831070/150‑1 (MQ=255)
                                                                                     aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCTTTGGTGGTGTTTATATGTTCTAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                  >  2:684378/1‑150 (MQ=255)
                                                                                     aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                  >  2:1155928/1‑150 (MQ=255)
                                                                                     aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                  <  2:326827/150‑1 (MQ=255)
                                                                                         ttAGCGGGGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGcc                                                              <  2:387855/150‑1 (MQ=255)
                                                                                           aGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCAt                                                            <  1:1633692/150‑1 (MQ=255)
                                                                                              ggTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCTGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATATTCTTAGGTGCCTTACCGCGCCATTGt                                                         >  1:689734/1‑150 (MQ=255)
                                                                                              ggTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACAGCGCCATTGt                                                         >  1:618006/1‑150 (MQ=255)
                                                                                               gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTTCTATAACCAGGCTGGCCTGGAGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                        <  2:467617/150‑1 (MQ=255)
                                                                                               gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCGATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                        <  1:208433/150‑1 (MQ=255)
                                                                                               gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                        >  1:1651094/1‑150 (MQ=255)
                                                                                               gTGAGGTTAGCCCTCGCACAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGAAATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTAACCGCGCCATTGTc                                                        <  1:684348/150‑1 (MQ=255)
                                                                                                tGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCg                                                       <  1:567743/150‑1 (MQ=255)
                                                                                                  aGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGAt                                                     <  1:640772/150‑1 (MQ=255)
                                                                                                   ggTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATa                                                    <  1:924670/150‑1 (MQ=255)
                                                                                                        gCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGc                                               <  1:862442/150‑1 (MQ=255)
                                                                                                           cTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGtt                                            >  1:1522774/1‑150 (MQ=255)
                                                                                                              gCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCa                                         >  1:1143700/1‑150 (MQ=255)
                                                                                                               cccAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTTGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAg                                        <  1:837399/150‑1 (MQ=255)
                                                                                                                   gctgctGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCt                                    <  1:1292773/150‑1 (MQ=255)
                                                                                                                      gctgcAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCg                                 <  1:1165453/150‑1 (MQ=255)
                                                                                                                       ctgcAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGATCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGc                                <  2:1075031/150‑1 (MQ=255)
                                                                                                                          cAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGt                             <  2:1562023/150‑1 (MQ=255)
                                                                                                                             ttACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTAc                          <  2:145747/150‑1 (MQ=255)
                                                                                                                             ttACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTAc                          >  1:1428738/1‑150 (MQ=255)
                                                                                                                               aCTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACct                        <  1:552314/150‑1 (MQ=255)
                                                                                                                                    tataACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                   >  1:964008/1‑150 (MQ=255)
                                                                                                                                    tataACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                   >  1:1673352/1‑150 (MQ=255)
                                                                                                                                    tataACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                   >  1:1666893/1‑150 (MQ=255)
                                                                                                                                    tataACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATCTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCTCTGTTACCTCTGGa                   >  2:741243/1‑150 (MQ=255)
                                                                                                                                      taACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTTGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGaaa                 <  1:1455404/150‑1 (MQ=255)
                                                                                                                                      taACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGaaa                 <  1:1474485/150‑1 (MQ=255)
                                                                                                                                        aCCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAATCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAgg               <  2:907637/150‑1 (MQ=255)
                                                                                                                                         aatggCTGGCCTGAGATATCGCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGa              >  2:242814/4‑150 (MQ=255)
                                                                                                                                         aatggCTGGCCTGAGATATCGCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGa              >  2:756524/4‑150 (MQ=255)
                                                                                                                                          cAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAAATCGTCGCTGTTACCTCTGGAAAGGAt             >  2:668967/1‑150 (MQ=255)
                                                                                                                                           aGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc            <  2:668348/150‑1 (MQ=255)
                                                                                                                                           aGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc            <  1:1546687/150‑1 (MQ=255)
                                                                                                                                            ggCTGGCCTGAGATATCGCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCg           >  2:145453/1‑150 (MQ=255)
                                                                                                                                              ctggcctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGcc         >  1:1414415/1‑150 (MQ=255)
                                                                                                                                              ctggcctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGcc         <  1:937437/150‑1 (MQ=255)
                                                                                                                                              ctggcctgagATATCGCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCGTCGGCGTAATCTTAGGTGCCTTACCGCTCCATTGTCGATACAGGCGTTCCTTATCGGCGCTGTTTGCTCTGGAAAGGATAGcc         >  2:297071/1‑150 (MQ=255)
                                                                                                                                               tggcctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCt        >  2:442071/1‑150 (MQ=255)
                                                                                                                                                  cctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcgc     >  1:998738/1‑150 (MQ=255)
                                                                                                                                                     ggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGaa  >  1:1504184/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CGGGTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAA  >  NC_000913/1625155‑1625453

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: