Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,406,098 C→T 6.6% D40D (GAC→GAT)  smrA → DNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,406,0980CT 6.6% 375.9 / 9.3 120D40D (GAC→GAT) smrADNA endonuclease
Reads supporting (aligned to +/- strand):  ref base C (53/59);  new base T (4/4);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 123 total aligned reads displayed.

ATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAA  >  NC_000913/1405949‑1406244
                                                                                                                                                     |                                                                                                                                                  
atCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                                                                    >  1:1854169/1‑150 (MQ=255)
atCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCAGCTGATACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                                                                    >  2:735277/1‑150 (MQ=255)
  cgcgCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGacac                                                                                                                                                  <  1:1063508/150‑1 (MQ=255)
    cgcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTAGACGCCAAGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTCTCGACACGc                                                                                                                                                >  2:1861248/1‑150 (MQ=255)
     gcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTACTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGGGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCt                                                                                                                                               <  2:333318/150‑1 (MQ=255)
        ctTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCa                                                                                                                                            >  1:1323373/1‑150 (MQ=255)
         ttAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                                                                           <  1:530202/150‑3 (MQ=255)
         ttAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCAGCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                                                                           <  1:1275822/150‑3 (MQ=255)
          tacttGCGGGAGGAGGTACCATGACACTTGACGACAAATCGCTGTTTCTAGACGCCATGGACGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                                                                          <  1:353812/147‑1 (MQ=255)
          tAATTGCGGGAGGATGTAACATGAACCTTGACGACAATTCGCTTTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                                                                          <  1:1175217/150‑1 (MQ=255)
          tAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                                                                          <  1:748553/150‑1 (MQ=255)
          tAATTGCGGGAGGAGGTAACATGAACGTTGACGACAAATCGCTGTTTCTTGACGCCATGGCAGATGTCCAGCCGCTGAAACGTGCTACCGATGACCACTGGCTTCCATCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                                                                          <  1:1281226/150‑1 (MQ=255)
           aaTAGCGGGAGGATGTAACATGAACCTAGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                                                                         <  1:1562078/150‑1 (MQ=255)
             ttGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTg                                                                                                                                       <  1:1829538/150‑1 (MQ=255)
              tGCGGGAGGATGTAACATGAACGTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGGGCTACCGATGTCCACTGGCATCCAACGCGTAACCGACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                                                                      <  2:190921/150‑1 (MQ=255)
                  ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                                                                                  >  2:1691237/1‑150 (MQ=255)
                  ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                                                                                  >  1:1416217/1‑150 (MQ=255)
                  ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACACCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGAAGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGCCACGCGGCAGCTTGATAAt                                                                                                                                  >  1:221826/1‑150 (MQ=255)
                      gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCCGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                              >  1:158133/1‑150 (MQ=255)
                      gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                              <  2:1699018/150‑1 (MQ=255)
                      gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                              <  1:506559/150‑1 (MQ=255)
                      gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGAGGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                              <  2:907562/150‑1 (MQ=255)
                             aaTGAACCTTGAAGTCAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCAACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGTCACGCTGCAGCTTGATAATTTCCTcacctc                                                                                                                       <  2:45078/149‑3 (MQ=255)
                                  aCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAt                                                                                                                  >  1:885958/1‑150 (MQ=255)
                                  aCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGATGCAGCTTGATAATTTCCTCACAACCGGAt                                                                                                                  >  2:1320166/1‑150 (MQ=255)
                                   ccTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAtt                                                                                                                 >  1:658277/1‑150 (MQ=255)
                                        acgacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                            >  2:1753685/1‑150 (MQ=255)
                                             aaaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatc                                                                                                       <  2:1797869/150‑1 (MQ=255)
                                              aaTCGCTGTTTCTTGACGCCTTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCCCCACCGGATTTCTCGAcatca                                                                                                      >  2:1422393/1‑150 (MQ=255)
                                               aTCGCTGTTTCTTGACGCCATGGAAGATGACCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatcat                                                                                                     <  2:1112604/150‑1 (MQ=255)
                                                  gCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                                                                  <  2:25339/150‑1 (MQ=255)
                                                    tGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                <  2:1191490/150‑1 (MQ=255)
                                                    tGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                                <  1:1651728/150‑1 (MQ=255)
                                                      tttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                              <  1:424101/150‑1 (MQ=255)
                                                         cTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGt                                                                                           <  2:3990/150‑1 (MQ=255)
                                                          ttGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTc                                                                                          >  2:1137610/1‑150 (MQ=255)
                                                           tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                                                                         <  2:1598891/150‑1 (MQ=255)
                                                           tGACGCCATGGAAGATGTACAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                                                                         <  2:1591343/150‑1 (MQ=255)
                                                            gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAg                                                                                        <  2:1058062/150‑1 (MQ=255)
                                                              cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                                                                      >  1:799380/1‑150 (MQ=255)
                                                                ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                                                                    <  1:1216678/150‑1 (MQ=255)
                                                                 cATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCt                                                                                   >  1:1414066/1‑150 (MQ=255)
                                                                  aTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTg                                                                                  <  1:1864789/150‑1 (MQ=255)
                                                                   tGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTgg                                                                                 >  2:1398316/1‑150 (MQ=255)
                                                                     gAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCAACGTATTTCTCGCCATCATCCCACTAAGTCAGCAGCTGGAg                                                                               >  1:917525/1‑150 (MQ=255)
                                                                          tGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcg                                                                          <  2:1399686/150‑1 (MQ=255)
                                                                              cAGCCGCTGAATCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGAGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCAGGAGTTTcggcgg                                                                      <  1:1361538/150‑1 (MQ=255)
                                                                              cAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTACCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCCGCCGCTGGAGTTTcggcgg                                                                      >  2:546307/1‑150 (MQ=255)
                                                                                gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGa                                                                    <  2:94432/150‑1 (MQ=255)
                                                                                   gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                                                                 >  2:95002/1‑150 (MQ=255)
                                                                                    cTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAggg                                                                >  1:1563606/1‑150 (MQ=255)
                                                                                     tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                                                               >  1:1800648/1‑150 (MQ=255)
                                                                                      gAAACGTGCTACCGATGTCCACTGGCATCCACCGCGTAACCAACGTGCGCCGCAGCGTATCCCCCCCCTGCAGCTTGCTAATTTCCTCACCCCCGGAGGTCGGGACATCTTCCCACTTCGTCAGACGCGGGAGTTTCGGCAGGAAGGGtt                                                              >  1:302101/1‑150 (MQ=255)
                                                                                      gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt                                                              <  2:1323373/150‑1 (MQ=255)
                                                                                      gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt                                                              >  2:1056908/1‑150 (MQ=255)
                                                                                        aaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGc                                                            <  1:1166981/150‑1 (MQ=255)
                                                                                         aCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGCCATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGCGTTGCa                                                           >  2:1576754/1‑150 (MQ=255)
                                                                                               tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgg                                                     <  2:1302015/150‑1 (MQ=255)
                                                                                               tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgg                                                     <  2:409384/150‑1 (MQ=255)
                                                                                                aCCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg                                                    >  1:1195288/1‑150 (MQ=255)
                                                                                                 ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg                                                   >  1:1214211/1‑150 (MQ=255)
                                                                                                  cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                  <  2:1344931/150‑1 (MQ=255)
                                                                                                   gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACCTCATCCCACTAAGTAAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACAGgggggg                                                 >  2:487248/1‑148 (MQ=255)
                                                                                                     tGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTATGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCt                                               <  1:238441/150‑1 (MQ=255)
                                                                                                       tCCACTGGCATCCAACGCGTACCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCATGATAAATTCCTCACCACCGGAGTTCTCGAAAACATCCCACTAAGTCAGCCGCTTGAGTTTCGGCGGGAAGGGTTGCAACAAGGGGTGCTgg                                             >  2:524165/1‑150 (MQ=255)
                                                                                                         cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                           >  2:1163752/1‑150 (MQ=255)
                                                                                                         cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                           <  2:618635/150‑1 (MQ=255)
                                                                                                          aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGTGTGCTGGATa                                          <  1:735277/150‑1 (MQ=255)
                                                                                                           cTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCCCCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATaa                                         >  1:1708635/1‑150 (MQ=255)
                                                                                                             ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                       >  2:390952/1‑150 (MQ=255)
                                                                                                             ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                       <  1:257605/150‑1 (MQ=255)
                                                                                                                  ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                  >  1:1719375/1‑150 (MQ=255)
                                                                                                                    aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAg                                >  2:1520995/1‑150 (MQ=255)
                                                                                                                    aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAg                                <  2:1171926/150‑1 (MQ=255)
                                                                                                                    aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAg                                <  1:1398316/150‑1 (MQ=255)
                                                                                                                     aCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGt                               >  1:587261/1‑150 (MQ=255)
                                                                                                                     aCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGGTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGt                               <  1:1247255/150‑1 (MQ=255)
                                                                                                                       gcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTgg                             <  2:1404837/150‑1 (MQ=255)
                                                                                                                        cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                            >  2:692417/1‑150 (MQ=255)
                                                                                                                        cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                            <  2:522155/150‑1 (MQ=255)
                                                                                                                         gTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTa                           >  2:1321876/1‑150 (MQ=255)
                                                                                                                          tAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTaa                          <  1:1605591/150‑1 (MQ=255)
                                                                                                                           aaCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTAGACATCATCCCAATAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCACTGGTaaa                         >  1:963109/1‑150 (MQ=255)
                                                                                                                            aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATATCAACACCACCGGATTTCTCGACATCCTCCCACTAAGACAGCCGCTGGAGTTTCGGAGGCAAGGGATCCAACATGGGGTGAAGGATAAGCTGCGCAGTggaatat                        >  1:238241/1‑145 (MQ=255)
                                                                                                                                aCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATcc                    <  1:1188620/150‑1 (MQ=255)
                                                                                                                                 cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCg                   <  1:896857/150‑1 (MQ=255)
                                                                                                                                  gTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGc                  <  1:1691237/150‑1 (MQ=255)
                                                                                                                                    gcgcCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGCCATCATCCCCCTAAGTCAGCCTCTGGCGTTTCGGCGGTAATGGTGGGACCATGGGGTGCTGGATAAGCAGCGCAGTGCGAAATAGCCGCaa                >  1:382492/1‑150 (MQ=255)
                                                                                                                                     agctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTACGGCGGGAAGGGTTGCAAGATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAAc               >  2:276527/5‑150 (MQ=255)
                                                                                                                                       ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCCCCACCGGATTTCTCGACATCCTCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAg             >  2:1738126/3‑150 (MQ=255)
                                                                                                                                        tgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAAATGGGGTGGTGGATAAGCTGCGCTGTGGTAAATATCCGCAACAgg            >  2:920567/2‑150 (MQ=255)
                                                                                                                                          cAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCg          <  1:1058903/150‑1 (MQ=255)
                                                                                                                                            gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAg        <  1:1137610/150‑1 (MQ=255)
                                                                                                                                               tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAACTCAGCCGCTGGAGGTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt     <  1:1154378/150‑1 (MQ=255)
                                                                                                                                                aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCCGATTTCTCGACATCATCCCACTATGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCTACAGGCGAGCCTg    <  2:1233449/150‑1 (MQ=255)
                                                                                                                                                 tCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCCCTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGATGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGa   >  1:58361/1‑150 (MQ=255)
                                                                                                                                                 tCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGa   >  1:9454/1‑150 (MQ=255)
                                                                                                                                                  cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCCTCCCACTAAGTCTGCCGCTGGAGTTTCGGCGGGTAGGGTTGCAACATGGGGTGCGCGATAAGCTGACCGGTGGTATATATGCGCAACAGGCGAGCATGaa  >  2:90120/1‑150 (MQ=255)
                                                                                                                                                  cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGaa  >  1:1631176/1‑150 (MQ=255)
                                                                                                                                                  cGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGaa  >  1:366737/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAA  >  NC_000913/1405949‑1406244

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