Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,488,989 G→T 11.4% H398N (CAC→AAC)  argD ← N‑acetylornithine aminotransferase/N‑succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,488,9890GT11.4% 458.0 / 29.0 146H398N (CAC→AAC) argDN‑acetylornithine aminotransferase/N‑succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (60/69);  new base T (7/10);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

CGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACTGTG  >  NC_000913/3488840‑3489138
                                                                                                                                                     |                                                                                                                                                     
cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                       >  1:74055/1‑150 (MQ=255)
cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                       >  2:1953264/1‑150 (MQ=255)
 gACCGCATCACGGGATTTGGCTGAGTCGGAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgg                                                                                                                                                      <  1:1599142/150‑1 (MQ=255)
 gACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgg                                                                                                                                                      <  1:1215389/150‑1 (MQ=255)
  aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTggg                                                                                                                                                     <  2:383887/150‑1 (MQ=255)
   ccGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGc                                                                                                                                                    <  1:2217485/150‑1 (MQ=255)
   ccGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGc                                                                                                                                                    >  2:192577/1‑150 (MQ=255)
    cGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCg                                                                                                                                                   <  2:298928/150‑1 (MQ=255)
        tCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaagg                                                                                                                                               <  1:2413931/150‑3 (MQ=255)
        tCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaagg                                                                                                                                               <  1:1429322/150‑3 (MQ=255)
         cACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggt                                                                                                                                              <  1:2180375/150‑4 (MQ=255)
         cACGGGATTTGGCGGAGTCGCAAGTGGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCACCCACCTTCGCCCACGCGGTGGCGaaggt                                                                                                                                              <  2:2278390/150‑4 (MQ=255)
          aCGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtg                                                                                                                                             <  1:2063083/150‑5 (MQ=255)
          aCGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtg                                                                                                                                             <  2:708114/150‑5 (MQ=255)
            gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtggt                                                                                                                                           <  2:381678/150‑7 (MQ=255)
            gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTg                                                                                                                                           <  2:1132567/150‑1 (MQ=255)
            gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTg                                                                                                                                           <  2:2176299/150‑1 (MQ=255)
            gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCAGTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTg                                                                                                                                           <  1:1896509/150‑1 (MQ=255)
               aTTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCAt                                                                                                                                        <  1:1232067/150‑1 (MQ=255)
                  tGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATccc                                                                                                                                     <  2:1317887/150‑1 (MQ=255)
                   ggCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCt                                                                                                                                    >  1:1685700/1‑150 (MQ=255)
                       gAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCAt                                                                                                                                <  1:617927/150‑1 (MQ=255)
                        aGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATc                                                                                                                               <  1:1932486/150‑1 (MQ=255)
                         gTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCg                                                                                                                              <  2:1764428/150‑1 (MQ=255)
                         gTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCg                                                                                                                              >  1:2297274/1‑150 (MQ=255)
                         gTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCg                                                                                                                              >  2:544482/1‑150 (MQ=255)
                           cGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACAACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGat                                                                                                                            >  2:1766426/1‑150 (MQ=255)
                            gCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGata                                                                                                                           >  1:1133250/1‑150 (MQ=255)
                              aaGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATc                                                                                                                         <  1:769058/150‑1 (MQ=255)
                                  tGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCa                                                                                                                     <  1:74287/150‑1 (MQ=255)
                                    cGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATc                                                                                                                   <  1:1933937/150‑1 (MQ=255)
                                    cGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTAACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATc                                                                                                                   <  2:2349751/150‑1 (MQ=255)
                                        atgcgtcggCGTTATGGAGGCTAAACCCTCGGATAAGGCATTCACCCCGCTTGCGACATCTTTTGCCGGTTGCTTCGTATTTACGCCCCAAGCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTcc                                                                                                               <  2:1568620/143‑1 (MQ=255)
                                        aTTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTcc                                                                                                               <  1:476300/150‑1 (MQ=255)
                                           cGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTccacc                                                                                                            >  2:142386/1‑150 (MQ=255)
                                                 cGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGc                                                                                                      >  2:104807/1‑150 (MQ=255)
                                                  gTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCg                                                                                                     <  2:196475/150‑1 (MQ=255)
                                                       gCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGGTGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGt                                                                                                >  2:2018102/1‑150 (MQ=255)
                                                         aGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGc                                                                                              <  1:2468504/150‑1 (MQ=255)
                                                          ggCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCa                                                                                             >  2:1004535/1‑150 (MQ=255)
                                                          ggCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCa                                                                                             >  1:1243928/1‑150 (MQ=255)
                                                          ggCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCa                                                                                             >  1:45606/1‑150 (MQ=255)
                                                           gCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCAGCTTCCACCACCAGCGACGGTGCaa                                                                                            <  2:1103659/150‑1 (MQ=255)
                                                            cTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCaaa                                                                                           >  1:1607657/1‑150 (MQ=255)
                                                              aaaCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAAc                                                                                         <  1:909488/150‑1 (MQ=255)
                                                                 ccATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCa                                                                                      >  1:1371335/1‑150 (MQ=255)
                                                                   aTCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATc                                                                                    <  1:2113382/150‑1 (MQ=255)
                                                                   aTCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATc                                                                                    <  2:1940019/150‑1 (MQ=255)
                                                                      ggATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATcaca                                                                                 >  2:126093/1‑150 (MQ=255)
                                                                        aTAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtc                                                                               >  2:1736212/1‑150 (MQ=255)
                                                                         tAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCATCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtcc                                                                              <  2:234957/150‑1 (MQ=255)
                                                                           aGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccgg                                                                            <  2:1731616/150‑1 (MQ=255)
                                                                              cATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccggtcc                                                                         >  2:52705/1‑150 (MQ=255)
                                                                               aTTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccggtccg                                                                        <  1:31264/150‑1 (MQ=255)
                                                                                 tCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGc                                                                      <  1:1004535/150‑1 (MQ=255)
                                                                                  cACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCa                                                                     >  1:1128884/1‑150 (MQ=255)
                                                                                  cACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCa                                                                     >  2:483800/1‑150 (MQ=255)
                                                                                   aCGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCAt                                                                    <  2:45606/150‑1 (MQ=255)
                                                                                   aCGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCAt                                                                    <  1:1463802/150‑1 (MQ=255)
                                                                                    cgccgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACAAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCAtt                                                                   >  2:75208/1‑150 (MQ=255)
                                                                                     gccgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTc                                                                  <  1:1673663/150‑1 (MQ=255)
                                                                                        gcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGc                                                               <  1:579974/150‑1 (MQ=255)
                                                                                         cATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCa                                                              <  1:188410/150‑1 (MQ=255)
                                                                                         cATCCGACATCTTTTGCCTGATGATTCGTATTTACGCCCCTACCACGTTCGCCACCGCGTGGGCGATACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCa                                                              <  2:74055/150‑1 (MQ=255)
                                                                                          aTCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCAc                                                             <  2:2073716/150‑1 (MQ=255)
                                                                                                  tCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGc                                                     >  2:2273968/1‑150 (MQ=255)
                                                                                                    ttttGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCa                                                   <  2:2456414/150‑1 (MQ=255)
                                                                                                          cTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGACGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctct                                             >  2:1543417/1‑150 (MQ=255)
                                                                                                           tGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTg                                            <  2:1689405/150‑1 (MQ=255)
                                                                                                            gATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTgc                                           >  2:176420/1‑150 (MQ=255)
                                                                                                                  tCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGc                                     <  1:2148813/150‑1 (MQ=255)
                                                                                                                     tATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATa                                  >  1:2469598/1‑150 (MQ=255)
                                                                                                                      aTTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATAc                                 >  2:1390016/1‑150 (MQ=255)
                                                                                                                        ttACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAg                               <  1:1977216/150‑1 (MQ=255)
                                                                                                                        ttACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAg                               >  1:2389761/1‑150 (MQ=255)
                                                                                                                         tACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAgg                              <  1:142386/150‑1 (MQ=255)
                                                                                                                           cGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCTTACAGGaa                            >  2:1274975/1‑150 (MQ=255)
                                                                                                                           cGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaa                            >  2:1819736/1‑150 (MQ=255)
                                                                                                                            gCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaaa                           >  2:818768/1‑150 (MQ=255)
                                                                                                                            gCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaaa                           >  2:103947/1‑150 (MQ=255)
                                                                                                                             ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCCCCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                          >  2:913906/1‑150 (MQ=255)
                                                                                                                             ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                          >  2:1412969/1‑150 (MQ=255)
                                                                                                                             ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                          >  2:2407988/1‑150 (MQ=255)
                                                                                                                              cccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGCCGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCCCCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATc                         >  2:1138061/1‑150 (MQ=255)
                                                                                                                               ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCCTCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                        >  1:1976379/1‑150 (MQ=255)
                                                                                                                                cAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAc                       >  1:1354324/1‑150 (MQ=255)
                                                                                                                                 aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg                      <  1:1215424/150‑1 (MQ=255)
                                                                                                                                 aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg                      <  2:789197/150‑1 (MQ=255)
                                                                                                                                 aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGAGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg                      <  1:1005223/150‑1 (MQ=255)
                                                                                                                                   ccaccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcgcg                    >  2:728478/1‑150 (MQ=255)
                                                                                                                                     cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcgcg                    >  1:972653/4‑148 (MQ=255)
                                                                                                                                     cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc                  >  2:54856/4‑150 (MQ=255)
                                                                                                                                     cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc                  >  1:2023897/4‑150 (MQ=255)
                                                                                                                                     cgtttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGcc                  >  1:810745/4‑150 (MQ=255)
                                                                                                                                        ttCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGa               >  1:275158/1‑150 (MQ=255)
                                                                                                                                         tCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGAc              <  2:874947/150‑1 (MQ=255)
                                                                                                                                         tCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGAc              >  2:316210/1‑150 (MQ=255)
                                                                                                                                            ccACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCtt           <  1:535554/150‑1 (MQ=255)
                                                                                                                                                  cgTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACt     <  1:884726/150‑1 (MQ=255)
                                                                                                                                                     gggCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACtgtg  <  2:796636/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTGTACTGTG  >  NC_000913/3488840‑3489138

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