Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,963 G→C 5.3% T65R (ACG→AGG) ‡ ygdB ← DUF2509 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,9630GC 5.3% 303.6 / 7.5 107T65R (ACG→AGG) ‡ygdBDUF2509 family protein
Reads supporting (aligned to +/- strand):  ref base G (56/45);  new base C (4/2);  total (60/47)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 108 total aligned reads displayed.

CAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAG  >  NC_000913/2962817‑2963108
                                                                                                                                                  |                                                                                                                                                 
cAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTc                                                                                                                                                 >  2:719837/1‑149 (MQ=255)
       tCCATCAATGACTTCGCCTGTTCGCCACAACGCAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                          >  2:1496001/1‑149 (MQ=255)
       tCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGGTTCAGCGTAATGCGAGCACTGAACTGCTGGCTGCGTCTGCcagc                                                                                                                                          >  1:1447342/1‑149 (MQ=255)
        ccATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAATCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagca                                                                                                                                         >  1:673366/1‑149 (MQ=255)
        ccATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCcagca                                                                                                                                         <  1:147565/149‑1 (MQ=255)
         cATCAATGACTTCGCCTGTTCGCCACAACGCAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGAAAACAAACCTGGGCATAGGTTTCAGCGTACTGCGAGCAATGAACTGCTGGCTGCGTCTGAcaggat                                                                                                                                        >  2:1483906/1‑146 (MQ=255)
           tCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGtt                                                                                                                                      <  1:440494/149‑2 (MQ=255)
             aaTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCa                                                                                                                                    >  1:374254/1‑149 (MQ=255)
             aaTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGACGCAGCCAGCAGTTCa                                                                                                                                    <  1:1472895/149‑1 (MQ=255)
              aTGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTACACGCAAACACACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGACCTGCTGGCTGCGTCTGCCAGCAGAGCAt                                                                                                                                   >  1:799360/1‑149 (MQ=255)
                gACTTCGCCTGTTCGCCACAACGAAACGCCGTCATAACCGGCACTCGAGAAGGCGTCATTCTCTGCCAGTAAACGCAAACAAACCGGGGCATCGGGTTCAGCGTACTGCGAGCACTGCACTGCTGGCTGCGTCTGCCAGCAGTGCAttt                                                                                                                                 <  1:1125534/149‑1 (MQ=255)
                  cTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAttttt                                                                                                                               >  1:63027/1‑149 (MQ=255)
                  cTTCGCCTGTTAGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCTGCAGTGCAatttt                                                                                                                               >  2:514357/1‑149 (MQ=255)
                   ttCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCGTCATTATCTGCCAGTAAACGCCAACACCCCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                              >  1:1550671/1‑149 (MQ=255)
                   ttCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGCCAATAAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                              <  2:565753/149‑1 (MQ=255)
                     cGCCTGTGCGCCACTACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTccc                                                                                                                            <  1:685700/149‑1 (MQ=255)
                      gCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTc                                                                                                                              <  1:1552832/146‑1 (MQ=255)
                          gTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTTTCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGc                                                                                                                       >  2:1438030/1‑149 (MQ=255)
                          gTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAAAAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGc                                                                                                                       <  2:1393898/149‑1 (MQ=255)
                              gCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTACCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAgc                                                                                                                   >  1:603207/1‑149 (MQ=255)
                                 acaACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                <  2:481496/149‑1 (MQ=255)
                                   aaCTAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACACACCTGGGCATGGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGa                                                                                                              <  1:277219/149‑1 (MQ=255)
                                   aaCGAAACGCCTTCATACCCGGCACTCAATAACGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGTTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGa                                                                                                              <  2:884051/149‑1 (MQ=255)
                                     cGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACt                                                                                                            <  2:913830/149‑1 (MQ=255)
                                      gAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTg                                                                                                           >  1:1466071/1‑149 (MQ=255)
                                            ccTTCATAACCGGCAATAAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGAt                                                                                                     <  2:433773/149‑1 (MQ=255)
                                             cTTCATCACCGGCACTCCATAAGGCTTCATTATCTGCCAGTCACCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGCCCGGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATg                                                                                                    <  2:625243/149‑1 (MQ=255)
                                                cATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGcc                                                                                                 >  1:1552798/1‑149 (MQ=255)
                                                 aTAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCCGCGCCGACTGAACGATGGCCt                                                                                                >  2:1034969/1‑149 (MQ=255)
                                                    aCCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCCGCGCCGACTGAACGATGGCCTggc                                                                                             >  1:831528/1‑149 (MQ=255)
                                                       ggcatacaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCTAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTggcggc                                                                                          <  1:1163990/143‑1 (MQ=255)
                                                       ggCAATCAATAAGGCTTCATTATCTGCCTGTAAACGCAAACAACCCTGGGCTTCGCTTTCAGCGTCCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCCGTGCAGTTTTCCCCACGCCAGCGCGGACTGACCGATGGCCTggcggc                                                                                          >  1:128765/1‑149 (MQ=255)
                                                       ggCAATCAATAAGGCTTCATTATCTGCCACTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGAGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGACTcgcggc                                                                                          >  2:473715/1‑149 (MQ=255)
                                                        gCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATGGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGGCTGCCCGCAGTGCATTTTTCCCCACGCCACCGCCGCCTGAACGCTGGCCAggcggcc                                                                                         >  2:388522/1‑148 (MQ=255)
                                                        gCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGTGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGgcggcg                                                                                         >  1:566826/1‑149 (MQ=255)
                                                         caatcaatAAGGCTTCATTATCTGCCAGTAAACGCAACCAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                <  1:785809/125‑1 (MQ=255)
                                                         caatcaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGcc                                                                                                                >  2:785809/1‑125 (MQ=255)
                                                         caatcaatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGgcgc                                                                                        <  1:627543/149‑1 (MQ=255)
                                                          aatcaatATGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACAGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGTCTGAACGATGGCCTGGCGGCGCa                                                                                       >  2:411583/1‑149 (MQ=255)
                                                              aatAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCTTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtga                                                                                   >  1:1263271/1‑149 (MQ=255)
                                                               atAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgac                                                                                  >  1:427029/1‑149 (MQ=255)
                                                                tAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgact                                                                                 >  1:916180/1‑149 (MQ=255)
                                                                tAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAtgact                                                                                 >  1:1214606/1‑149 (MQ=255)
                                                                    gCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATgactgact                                                                             >  1:285354/1‑149 (MQ=255)
                                                                        cATTATCTGCAAGTAAACGCAAACAACCCTGGGCAGCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTcc                                                                         <  2:1359771/149‑1 (MQ=255)
                                                                             tcagCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCt                                                                    <  1:195209/146‑1 (MQ=255)
                                                                                gCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGAACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTACCCCACGCCACCGCCGACTGCCCGAAGGCCTGGCGGCCCAATGACAGACTTTACATGCacac                                                                 >  1:1025227/1‑149 (MQ=255)
                                                                                   aGTAAACGCAAACAAACCTGGGCATCGGTTTCCGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGCCTGACTTTCCATGCTCAcgcg                                                              >  1:1459854/1‑149 (MQ=255)
                                                                                    gTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCga                                                             <  1:522292/149‑1 (MQ=255)
                                                                                    gTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTACCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAAAGACTGACTTTCCATGCTCAcgcca                                                             >  1:1307762/1‑147 (MQ=255)
                                                                                     tAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCgag                                                            <  1:901573/149‑1 (MQ=255)
                                                                                       aaCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGagag                                                          <  1:1319240/149‑1 (MQ=255)
                                                                                        aCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCTATGACTGACTTTCCATGCTCACGCGAGAgg                                                         >  1:1339934/1‑149 (MQ=255)
                                                                                          gCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCCCGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCGCACGCGAGAGGCa                                                       >  1:1568771/1‑149 (MQ=255)
                                                                                            aaacaagcCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaa                                                     >  2:700031/1‑149 (MQ=255)
                                                                                            aaacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaa                                                     >  2:449598/1‑149 (MQ=255)
                                                                                            aaacaaacCTGGGCATCGGATTCAGCGTACTGCGAGCACTGAACAGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACAGAACGATGGCCTGGCGGCGCAATGCCTGCATTTCCATGCTCACGCGAGAGGCaaa                                                     >  2:1534416/1‑149 (MQ=255)
                                                                                             aacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaaa                                                    >  1:1063076/1‑149 (MQ=255)
                                                                                             aacaaaaGTGGGCATCGGTTTCAGCGTACTGCGAGCTCTGAACTGCTGGCTGCGTCTTCCAGCAGTGAATTTTTCACCACGCCAGCGCCTTCTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCaaaa                                                    <  1:569825/149‑1 (MQ=255)
                                                                                              acaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAAc                                                   >  1:610881/1‑149 (MQ=255)
                                                                                              acaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAAc                                                   >  1:41133/1‑149 (MQ=255)
                                                                                               caaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAACCt                                                  >  2:477042/1‑149 (MQ=255)
                                                                                                    cTGGGCATCGGTTTCAGTGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGAt                                             <  2:566826/149‑1 (MQ=255)
                                                                                                    cTGGGCATCGGTTGCAGCGTACTGCGAGCACTGAACTGCTGGATGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGAt                                             <  2:890806/149‑1 (MQ=255)
                                                                                                         cATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCAGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCATTGACTGACTTTCCATGCTCAAGCGAGAGGCAAAACTGCGCTCCTGt                                        >  2:1241295/1‑149 (MQ=255)
                                                                                                          aTCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                       <  2:982930/149‑1 (MQ=255)
                                                                                                          aTCGGTTGCCGCGTACTGCGAGCACTGACCTGCTGGCTGCGTCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGtt                                       <  2:935712/149‑1 (MQ=255)
                                                                                                                ttCAGTGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCCTGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCa                                 >  2:1258160/1‑149 (MQ=255)
                                                                                                                  cAGCGTACTGCGAGCACTGAACTGCTGTCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCACTGACTGCCTTTCCATGCACACGCGAGAGGCAAAACGTCGATACTGTTGTCTCAtt                               >  1:1064395/1‑149 (MQ=255)
                                                                                                                   aGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTc                              <  1:638617/149‑1 (MQ=255)
                                                                                                                   aGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCGCATTc                              >  1:624145/1‑149 (MQ=255)
                                                                                                                    gCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATTGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTcc                             <  1:582545/149‑1 (MQ=255)
                                                                                                                    gCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTcc                             <  1:1438030/149‑1 (MQ=255)
                                                                                                                       tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCCGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTg                          <  2:1144604/149‑1 (MQ=255)
                                                                                                                       tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTg                          <  1:1488480/149‑1 (MQ=255)
                                                                                                                       tACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTg                          >  1:1334/1‑149 (MQ=255)
                                                                                                                        aCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGt                         <  1:668604/149‑1 (MQ=255)
                                                                                                                        aCTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGt                         >  2:786210/1‑149 (MQ=255)
                                                                                                                         cTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTa                        <  1:547620/149‑1 (MQ=255)
                                                                                                                          tGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTATCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCAAGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTACTTGTaa                       <  2:1447342/149‑1 (MQ=255)
                                                                                                                           gCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCGTGTTGGCTCATTCCTTGTAAt                      >  1:505824/1‑149 (MQ=255)
                                                                                                                            cGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCATTGACTGACTTTCCATGCTCACGCGAGAGGCAAACCTGCGATCCTGTTGGCTCATTCCTTGTAATa                     <  1:899990/149‑1 (MQ=255)
                                                                                                                                 aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCaaa                <  1:268388/149‑1 (MQ=255)
                                                                                                                                 aCTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGAATTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCaaa                >  1:362147/1‑149 (MQ=255)
                                                                                                                                 aCTGAACTGCTGGCTGCCTCAGCCAGCAGTGCATTTTTCCCCAAGCCAGAGCCGACTGAACGATTGCCTGCCGGCGCAATGAATGCCTATCCATGATCACGCGATAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGaaaa                >  2:494796/1‑149 (MQ=255)
                                                                                                                                   tGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGACGACTGAACGATGGCCTGGCGGCGCAATTACTGACTATCCATGATCAAGCTAGAGGCAAAACTGCGCTCGTGTTGGCTCATGCCTTGAAATAGCaaagt              >  2:869131/1‑147 (MQ=255)
                                                                                                                                    gcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAATGATGGCCTGGCGCCGCTATGACTGACTTTCCATGCTCACGCGAGCGGCATCACTGCGAGCCTGTTGGCTCTTTCCTTGTAATAGGAAACTa             >  1:178714/3‑149 (MQ=255)
                                                                                                                                    gcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAAGGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATGCCGTGTAATTGCCAACTa             >  1:184229/3‑149 (MQ=255)
                                                                                                                                    gcaCTGCTGGCAGACGCAGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTGGATATAGCACACTa             >  2:340150/3‑149 (MQ=255)
                                                                                                                                     aaCTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAc            <  1:1190034/149‑1 (MQ=255)
                                                                                                                                     aaCTGCTGGCTGCGTCTGACAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCATTGCCTGACGTTCCATGCTGCCGCGAGAGGCAAAACTGCCATCCTCTAGGCTCATTCCTTGTAATAGCAAACTAc            >  2:303701/1‑149 (MQ=255)
                                                                                                                                      aCTGCTGGCTGCGGCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAcc           <  1:1483906/149‑1 (MQ=255)
                                                                                                                                        tgctgGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTAcacc         >  1:915984/1‑146 (MQ=255)
                                                                                                                                        tgctgGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCa         >  2:613734/1‑149 (MQ=255)
                                                                                                                                        tgctgGCTGCGTCTGCCAGCAGTGCATATATCCCCACGCCCGCGCCGACTGAGCGATGGCCTGGCGGCGCAATGACTGACTTTCCATCCACACGCGAGAGGCAAAACTGCGATCCTGGTGGCTCAATCCTTGTAAAAGCAAACTACCCa         >  1:36597/1‑149 (MQ=255)
                                                                                                                                          gtggctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCaaa       <  2:1214606/148‑1 (MQ=255)
                                                                                                                                          ctggctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCaaa       <  2:1466071/149‑1 (MQ=255)
                                                                                                                                            ggctgcggctgcCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAAc     <  1:440567/149‑1 (MQ=255)
                                                                                                                                             gctgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATCGCCTGGCGGCGCAATGCGTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCACAAcc    >  1:964687/1‑149 (MQ=255)
                                                                                                                                               tgCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACcag  <  1:719837/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CAATGTTTCCATCAATGACTTCGCCTGTTCGCCACAACGAAACGCCTTCATAACCGGCAATCAATAAGGCTTCATTATCTGCCAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAG  >  NC_000913/2962817‑2963108

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: