Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,625,313 T→C 6.0% intergenic (+22/+22) 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,3130TC 6.0% 301.0 / 5.9 100intergenic (+22/+22)ydeJ/dcpinactive PncC family protein/dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (44/49);  new base C (3/3);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 104 total aligned reads displayed.

TGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCC  >  NC_000913/1625166‑1625462
                                                                                                                                                    |                                                                                                                                                      
tGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTCGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑Ca                                                                                                                                                     <  1:1311822/150‑1 (MQ=255)
 ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                                    >  2:1165453/1‑150 (MQ=255)
 ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAg                                                                                                                                                    >  1:193864/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: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:170404/1‑150 (MQ=255)
                                  cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                                   <  1:945039/150‑1 (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:326827/150‑1 (MQ=255)
                                                                          aGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACcg                                                                           >  2:1155928/1‑150 (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                                   >  1:1428738/1‑150 (MQ=255)
                                                                                                                  ttACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTAc                                   <  2:145747/150‑1 (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)
                                                                                                                                             gATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAAc        >  2:21567/1‑150 (MQ=255)
                                                                                                                                                    t‑cAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGcc  <  1:1522951/150‑1 (MQ=255)
                                                                                                                                                    t‑cAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGcc  >  1:1205329/1‑150 (MQ=255)
                                                                                                                                                    t‑cAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGcc  <  1:963569/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATC‑T‑CAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCC  >  NC_000913/1625166‑1625462

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