Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,306 C→A 8.1% intergenic (+15/+29) 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,3060CA 8.1% 270.9 / 10.5 96intergenic (+15/+29)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base C (40/48);  new base A (3/5);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 104 total aligned reads displayed.

GTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAA  >  NC_000913/1625158‑1625455
                                                                                                                                                    |                                                                                                                                                       
gTACCGTCTGGTTTGCGTGGCATATTCAAGGCCACACCTCCACTGCGGTTATGCAGTTTGCTGGCGACTGCGAACCGGTAGTAGCTTACGCGGTGAGGCTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCg                                                                                                                                                        <  2:349109/150‑1 (MQ=255)
    cGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGatat                                                                                                                                                    >  1:145747/1‑150 (MQ=255)
     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                                                                                                                                             >  1:193864/1‑150 (MQ=255)
         ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                             >  2:1165453/1‑150 (MQ=255)
         ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                             >  2:1019796/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)
                       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)
                                     cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGtta                                                                                                                 >  1:996121/1‑149 (MQ=255)
                                     cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                 >  2:504577/1‑150 (MQ=255)
                                     cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                 <  2:1514993/150‑1 (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)
                                                                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                                                                                >  1:512552/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                                                          >  1:1651094/1‑150 (MQ=255)
                                                                                            gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                          <  2:337181/150‑1 (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                            <  2:145747/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:1673352/1‑150 (MQ=255)
                                                                                                                                 tataACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGa                     >  1:1666893/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              <  1:1546687/150‑1 (MQ=255)
                                                                                                                                        aGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc              <  2:668348/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)
                                                                                                                                                    |                                                                                                                                                       
GTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAA  >  NC_000913/1625158‑1625455

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: ‑