Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,907,998 A→T 8.8% intergenic (‑57/+31) eno ← / ← pyrG enolase/CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,907,9980AT 8.8% 336.5 / 18.6 102intergenic (‑57/+31)eno/pyrGenolase/CTP synthetase
Reads supporting (aligned to +/- strand):  ref base A (55/37);  new base T (6/3);  total (61/41)
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.33e-01

Only 100 of 104 total aligned reads displayed.

AACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACC  >  NC_000913/2907851‑2908144
                                                                                                                                                   |                                                                                                                                                  
aaCCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAg                                                                                                                                                   >  3:674624/1‑149 (MQ=255)
aaCCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAg                                                                                                                                                   <  4:743013/149‑1 (MQ=255)
    accCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGt                                                                                                                                                >  2:16293/1‑148 (MQ=255)
    accCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGt                                                                                                                                                >  2:428919/1‑148 (MQ=255)
    accCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTa                                                                                                                                               >  4:202496/1‑149 (MQ=255)
     cccTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTa                                                                                                                                               <  2:28583/148‑1 (MQ=255)
         ccAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgcg                                                                                                                                           >  2:609752/1‑148 (MQ=255)
         ccAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgcg                                                                                                                                           <  2:216030/148‑1 (MQ=255)
          cAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGtt                                                                                                                                         >  4:515519/1‑149 (MQ=255)
            gATGTACTGCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGc                                                                                                                                       >  4:76637/1‑149 (MQ=255)
             aTGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGc                                                                                                                                       <  1:289108/148‑1 (MQ=255)
              tGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGcc                                                                                                                                      >  2:794143/1‑148 (MQ=255)
               gTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCg                                                                                                                                     <  1:794143/148‑1 (MQ=255)
                    tCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGtt                                                                                                                                         <  3:137608/139‑1 (MQ=255)
                    tCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                                                                >  1:157209/1‑148 (MQ=255)
                    tCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                                                                >  2:836530/1‑148 (MQ=255)
                    tCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                                                                >  2:291904/1‑148 (MQ=255)
                    tCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                                                                <  1:291904/148‑1 (MQ=255)
                         ttCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACtttt                                                                                                                          >  3:627339/1‑149 (MQ=255)
                         ttCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACtttt                                                                                                                          <  3:515519/149‑1 (MQ=255)
                            aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACtttttt                                                                                                                        <  1:428919/148‑1 (MQ=255)
                            aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTtta                                                                                                                       <  3:888592/149‑1 (MQ=255)
                            aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTtta                                                                                                                       <  3:76637/149‑1 (MQ=255)
                            aaCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTtta                                                                                                                       >  4:205593/1‑149 (MQ=255)
                             aCAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                                                                      >  4:355635/1‑149 (MQ=255)
                               aGTCGGGTTACCACGGGAGTCGATGATTTCACGATCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTGTttactt                                                                                                                    >  4:646853/1‑149 (MQ=255)
                                gTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTTACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                   >  4:590857/1‑149 (MQ=255)
                                gTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                   >  4:424701/1‑149 (MQ=255)
                                    ggTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTacttactt                                                                                                                >  2:894247/1‑148 (MQ=255)
                                       tACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGcc                                                                                                            <  3:355635/149‑1 (MQ=255)
                                       tACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGcc                                                                                                            <  3:670346/149‑1 (MQ=255)
                                        aCCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCGGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGcc                                                                                                            <  2:766231/148‑1 (MQ=255)
                                             gggAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACg                                                                                                       >  2:80503/1‑148 (MQ=255)
                                                     atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                                                                >  1:329358/1‑115 (MQ=255)
                                                     atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTa                                                                                                                                <  2:329358/115‑1 (MQ=255)
                                                       gatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACt                                                                                            >  4:114112/1‑149 (MQ=255)
                                                            cACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCt                                                                                        <  1:80503/148‑1 (MQ=255)
                                                                aCCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTggcg                                                                                    >  2:536788/1‑148 (MQ=255)
                                                                   gatgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGaa                                                                                              >  4:352889/1‑135 (MQ=255)
                                                                   gatgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGaa                                                                                              <  3:352889/135‑1 (MQ=255)
                                                                    atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCttt                                                                               <  4:627339/149‑1 (MQ=255)
                                                                    atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCttt                                                                               >  4:488529/1‑149 (MQ=255)
                                                                    atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCtt                                                                                >  2:28791/1‑148 (MQ=255)
                                                                     tgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                              >  3:893089/1‑149 (MQ=255)
                                                                      gatTTTTACTATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTca                                                                             <  3:590857/149‑1 (MQ=255)
                                                                      gatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                              <  1:609752/148‑1 (MQ=255)
                                                                         ttttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTcaca                                                                           >  1:531145/1‑148 (MQ=255)
                                                                          tttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaaa                                                                         <  3:205593/149‑1 (MQ=255)
                                                                          tttACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaa                                                                          <  1:536788/148‑1 (MQ=255)
                                                                            tACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGc                                                                       >  4:310170/1‑149 (MQ=255)
                                                                             aCGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCg                                                                                          <  1:422377/129‑1 (MQ=255)
                                                                             aCGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCg                                                                                          >  2:422377/1‑129 (MQ=255)
                                                                             aCGATTTTGAACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCGGGAACTCGCTGGCGGCTTTCACAAAGcc                                                                      >  4:579205/1‑149 (MQ=255)
                                                                              cGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCt                                                                     <  4:674624/149‑1 (MQ=255)
                                                                                         aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCgg                                                          >  3:891957/1‑149 (MQ=255)
                                                                                         aTTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCg                                                           <  1:377876/148‑1 (MQ=255)
                                                                                          ttAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCggg                                                         >  3:279215/1‑149 (MQ=255)
                                                                                              gTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                            <  3:508734/62‑1 (MQ=255)
                                                                                              gTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                            >  4:508734/1‑62 (MQ=255)
                                                                                                 ttCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa                                                   <  1:16293/148‑1 (MQ=255)
                                                                                                   ccTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCa                                                <  3:310170/149‑1 (MQ=255)
                                                                                                    cTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGt                                                                                                                                                <  2:476840/52‑1 (MQ=255)
                                                                                                    cTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCAc                                               >  3:844697/1‑149 (MQ=255)
                                                                                                    cTCAAGTCACGAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGt                                                                                                                                                >  1:476840/1‑52 (MQ=255)
                                                                                                     tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                                                                      >  3:856353/1‑77 (MQ=255)
                                                                                                     tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTCCACAAAGCCTGCAAACAGCGGGTGACCATCACg                                              >  3:405519/1‑149 (MQ=255)
                                                                                                      cAAGTCATTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGt                                             <  3:579205/149‑1 (MQ=255)
                                                                                                              tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTa                                                                                                                                               >  2:473061/1‑43 (MQ=255)
                                                                                                              tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTa                                                                                                                                               <  1:473061/43‑1 (MQ=255)
                                                                                                               aGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                              >  1:581717/1‑107 (MQ=255)
                                                                                                               aGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTc                                                                              <  2:581717/107‑1 (MQ=255)
                                                                                                                 ttAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGt                                  >  4:216077/1‑149 (MQ=255)
                                                                                                                  tAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTcaca                                                                           >  1:378663/1‑107 (MQ=255)
                                                                                                                  tAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTcaca                                                                           <  2:378663/107‑1 (MQ=255)
                                                                                                                           aCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCgg                         <  1:894247/148‑1 (MQ=255)
                                                                                                                             tCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTTCTTCACCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGat                                                     <  3:721064/118‑2 (MQ=255)
                                                                                                                             tCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCACCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGAc                                                     >  4:721064/1‑118 (MQ=255)
                                                                                                                              ccAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATg                      >  1:686243/1‑148 (MQ=255)
                                                                                                                              ccAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATg                      >  1:129161/1‑148 (MQ=255)
                                                                                                                                 gACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGaa                   >  2:245868/1‑148 (MQ=255)
                                                                                                                                 gACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAAc                  >  3:568047/1‑149 (MQ=255)
                                                                                                                                    aaacaaCGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTgg               >  4:203333/1‑149 (MQ=255)
                                                                                                                                      acaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGc              >  1:177315/1‑148 (MQ=255)
                                                                                                                                      acaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCa             >  4:57529/1‑149 (MQ=255)
                                                                                                                                       caacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCaa            <  3:488529/149‑1 (MQ=255)
                                                                                                                                         gcgcgTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGa                        >  4:581412/2‑135 (MQ=255)
                                                                                                                                         acGCGTACCCTGGGTACGCGTTGCCGCTCTATCTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGa                        <  3:581412/135‑1 (MQ=255)
                                                                                                                                         acGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc          >  3:790052/1‑149 (MQ=255)
                                                                                                                                          cgcgTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc          >  2:644690/1‑148 (MQ=255)
                                                                                                                                          cgcgTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc          >  1:523286/1‑148 (MQ=255)
                                                                                                                                           gcgTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc         >  1:164073/1‑148 (MQ=255)
                                                                                                                                           gcgTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCa        >  3:662765/1‑149 (MQ=255)
                                                                                                                                               aCCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGaa    <  4:810779/149‑1 (MQ=255)
                                                                                                                                                cccTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTgg                                           >  4:299556/1‑109 (MQ=255)
                                                                                                                                                cccTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAAc   >  4:595565/1‑149 (MQ=255)
                                                                                                                                                cccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAAc   >  4:402650/1‑149 (MQ=255)
                                                                                                                                                cccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAAc   >  3:798159/1‑149 (MQ=255)
                                                                                                                                                  cAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAAcc  >  2:128355/1‑148 (MQ=255)
                                                                                                                                                   aGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc         <  1:636508/140‑1 (MQ=255)
                                                                                                                                                   aGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc         >  2:636508/1‑140 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACC  >  NC_000913/2907851‑2908144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: