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.

GTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAA  >  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)
        tGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTCGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCa                                                                                                                                              <  1:1311822/150‑1 (MQ=255)
         ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                             >  1:193864/1‑150 (MQ=255)
         ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                             >  2:1165453/1‑150 (MQ=255)
         ggTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAg                                                                                                                                             >  2:1019796/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)
                       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)
                                     cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGtta                                                                                                                 >  1:996121/1‑149 (MQ=255)
                                     cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                 >  2:504577/1‑150 (MQ=255)
                                     cTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTc                                                                                                                 <  2:1514993/150‑1 (MQ=255)
                                      tACACTTCGGTTATGCATTTTGCGGGCGACTGCGCAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCACAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCa                                                                                                                <  2:917933/150‑1 (MQ=255)
                                       acacTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGGGTTTATATGTTCaa                                                                                                               >  1:830776/1‑150 (MQ=255)
                                          cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                            <  1:945039/150‑1 (MQ=255)
                                          cTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGcc                                                                                                            >  1:170404/1‑150 (MQ=255)
                                             cGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGGGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACg                                                                                                         <  1:64208/150‑1 (MQ=255)
                                                    gCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgc                                                                                                  >  2:1612769/1‑150 (MQ=255)
                                                    gCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTgc                                                                                                  >  1:668348/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)
                                                           gCTGGCGACTGCGAAACGGTATTAGCTTTAACGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCg                                                                                           >  1:659051/1‑150 (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)
                                                                                            gTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTc                                                          <  2:337181/150‑1 (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:1579414/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              <  1:1546687/150‑1 (MQ=255)
                                                                                                                                        aGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATc              <  2:668348/150‑1 (MQ=255)
                                                                                                                                         ggCTGGCCTGAGATATCGCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCg             >  2:145453/1‑150 (MQ=255)
                                                                                                                                           ctggcctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGcc           <  1:937437/150‑1 (MQ=255)
                                                                                                                                           ctggcctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGcc           >  1:1414415/1‑150 (MQ=255)
                                                                                                                                           ctggcctgagATATCGCCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCGTCGGCGTAATCTTAGGTGCCTTACCGCTCCATTGTCGATACAGGCGTTCCTTATCGGCGCTGTTTGCTCTGGAAAGGATAGcc           >  2:297071/1‑150 (MQ=255)
                                                                                                                                            tggcctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCt          >  2:442071/1‑150 (MQ=255)
                                                                                                                                              gcctggcGATATCTCAGGACAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcg        <  2:1178093/150‑1 (MQ=255)
                                                                                                                                               cctggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTcgc       >  1:998738/1‑150 (MQ=255)
                                                                                                                                                  ggcGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGaa    >  1:1504184/1‑150 (MQ=255)
                                                                                                                                                    cGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGaaaa  >  1:1380622/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GTACCGTCTGGTTTGCGTGGCATATTAAAGGCCAGAACTACACTGCGGTTATGCATTTTGCTGGCGACTGCGAAACGGTATTAGCTTTAGCGGTGAGGTTTGCCCTCGCCCAGCTGCTGCAATTACTGCTATAACCAGGCTGGCCTGGCGATATCTCAGGCCAGCCATTGGTGGTGTTTATATGTTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAA  >  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: