Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,749,326 A→C D32A (GAT→GCT)  sgbU → putative L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,749,3260AC98.7% 244.5 / NA 74D32A (GAT→GCT) sgbUputative L‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (40/33);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑ATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCAT  >  NC_000913/3749094‑3749557
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
ccctTGCGGACTTGCCGATGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑Ctt                                                                                                                                                                                                                                         <  1:288849/235‑1 (MQ=255)
ccccTGCTGACCTGCCGCTGTTGAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑Ctt                                                                                                                                                                                                                                         <  1:705709/235‑1 (MQ=255)
      ctgACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGGGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAACGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑CAAAGCTGCGGTGTTG‑CTTTTGTCGa                                                                                                                                                                                                                                  >  1:357800/1‑236 (MQ=255)
          ccTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCCTGCGCCAATCGACGCTATCTGGGGAGGTGCGCGGGCGTAATCCTCAGTTAGGGATTTATGGAAAAGCGCTGGCGAAAGGTCTCTCCCGGCCGGGGCGGCGTGGGCTGGCACAAACGCTGGGGGTTGG‑CTTTTG‑CGAAATg                                                                                                                                                                                                                              >  2:470734/1‑236 (MQ=255)
              cggCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCg                                                                                                                                                                                                                           <  2:221909/233‑1 (MQ=255)
              ccGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGt                                                                                                                                                                                                                             <  2:501913/233‑1 (MQ=255)
              ccGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGt                                                                                                                                                                                                                             >  1:501913/1‑233 (MQ=255)
                              tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGGTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACcgacga                                                                                                                                                                                                          >  1:60714/1‑236 (MQ=255)
                              tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGGAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGGGCGGCTGGGGCTGGCA‑AAAAACTGCGGGTTTG‑CTTTTGTCGTAAAGTCGGTGGGTGGAAacgacgc                                                                                                                                                                                                          >  2:198147/1‑235 (MQ=255)
                              tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGGGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGAGGCTGGGGCTTGGA‑AAAAGCGGCGGGTTTGCCTTTTGTTGAAAGTTGGGTGGTGGAATCcgcga                                                                                                                                                                                                           >  2:80613/1‑233 (MQ=255)
                                      gaaGTGTTTACTGCCGGGCGGCCACTGGCAGGCGCGGCAACTCCGGCGCAAGTGGCTGGCGATTTCCATGCGAAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGctct                                                                                                                                                                                                  <  2:220265/235‑1 (MQ=255)
                                      aaaGCGTTTATTGCCGGGCGGGCACGGGCAGGCTCGGCAAATCCGGCGGAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGctct                                                                                                                                                                                                  <  2:99205/236‑1 (MQ=255)
                                              tagtgcagagcgggcACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTc                                                                                                                                                                                            <  2:626593/225‑1 (MQ=255)
                                              tATTGCCGGGCGGGCCCTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTc                                                                                                                                                                                            <  2:296274/234‑1 (MQ=255)
                                                           gcACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGcgcc                                                                                                                                                                              <  1:620358/235‑1 (MQ=255)
                                                             aCTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGca                                                                                                                                                                           >  1:596321/1‑236 (MQ=255)
                                                                           cAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCC‑AAAAGCTGCGGGTTCT‑CTTTTGTCGAAATGTCGGTGGGTGAAACCGACGACAGGCTCTCACGTCTTGCCTTGAGCGCCGCGCCAAgcatgttgct                                                                                                                                                             >  2:200597/1‑227 (MQ=255)
                                                                                   gcgcAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGGGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGAGCCGTACAAAGGACTTCGCTGG‑TTGCcg                                                                                                                                                    >  1:716906/1‑236 (MQ=255)
                                                                                   gcgcAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGGTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAAACGACGAACGGCCCTCACGGCTTGACTGGGAGGCCGGCCCAAGGGCTTCGCCGG‑TTGCcg                                                                                                                                                    >  2:72043/1‑236 (MQ=255)
                                                                                            gctggcGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGGGGGTGAAACCGACGGAAGGGCCTTCCGGCTTGACTGGGGTGCCGCACAAAGGACTTGCCTGA‑GTGCcgccgtggttg                                                                                                                                           >  2:505107/1‑228 (MQ=255)
                                                                                                gctaTTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGaaa                                                                                                                                        <  1:113937/232‑1 (MQ=255)
                                                                                                           gcgcAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGTGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGa                                                                                                                            <  1:505107/236‑1 (MQ=255)
                                                                                                                 aTCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGt                                                                                                                      >  1:234327/1‑236 (MQ=255)
                                                                                                                         tATCTGGGGAGGTTCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GAtgtg                                                                                                              <  2:60714/235‑1 (MQ=255)
                                                                                                                         tATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GAtgtg                                                                                                              <  1:569735/235‑1 (MQ=255)
                                                                                                                         tATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GAtgtg                                                                                                              <  1:207628/235‑1 (MQ=255)
                                                                                                                           tCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCtctc                                                                                                                                                                                                 <  1:241520/152‑1 (MQ=255)
                                                                                                                           tCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCtctc                                                                                                                                                                                                 >  2:241520/1‑152 (MQ=255)
                                                                                                                             gggggAGGTGCGCGTGCGTAATCATCAGGTAGGGAGTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GAtgtg                                                                                                              <  1:200597/230‑1 (MQ=255)
                                                                                                                                 gAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGCTGATCGAAACAGGGGTTGGGATCCCGTC‑GATGTGCCTGTCCGc                                                                                                     >  1:235383/1‑236 (MQ=255)
                                                                                                                                  aGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAgcg                                                                                                                                                                                <  2:553337/162‑1 (MQ=255)
                                                                                                                                  aGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAgcg                                                                                                                                                                                >  1:553337/1‑162 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCC‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACCGCTCTCAAGTCTTGACTGGAGCGCCTCACAAAGGACTTCGCTGGCTTCCCGC‑ATGACCGGAACAGGGGTTTGGATCCCATCACATGT‑CCTGTCCGcaca                                                                                                  >  2:338194/1‑235 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTTACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGAACGAAACAGGTGTTGGGATCCCGTG‑CACGTGCCTGTCCGtacac                                                                                                 >  2:634604/1‑235 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCCCAAAGGACTTCGCTGG‑TTGCCGCGATGATCGGAAACGGTGTTCGGATCCCGTC‑GATGTGCCTGTCCGCACAt                                                                                                 >  2:256650/1‑236 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATcc                                                                                                                         >  1:1432/1‑213 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATcc                                                                                                                         <  2:1432/213‑1 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCTTC‑GATGGGCCTGTCCGcaca                                                                                                  >  1:698704/1‑235 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACGTGAGCGCCGCCCAAAGGACTTCGCTGG‑GTGCCGCGATGGTCGAAACAGGGGTTGGGATCCCGTC‑GATGGTCCTggtcgcgca                                                                                                  >  2:643202/1‑227 (MQ=255)
                                                                                                                                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGc                                                                                                             <  1:356456/222‑1 (MQ=255)
                                                                                                                                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGc                                                                                                             >  2:356456/1‑222 (MQ=255)
                                                                                                                                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACAtc                                                                                                <  1:219743/235‑1 (MQ=255)
                                                                                                                                         cgtgcgtAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACAtcgtc                                                                                             >  1:421086/1‑236 (MQ=255)
                                                                                                                                            gcgtAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGc                                                                                           <  2:596321/235‑1 (MQ=255)
                                                                                                                                              gtAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACAtc                                                                                                >  1:118354/1‑228 (MQ=255)
                                                                                                                                               tAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACAtc                                                                                                >  1:736811/1‑227 (MQ=255)
                                                                                                                                                  tcatcaGTTAGGGATTTATGAAAAAGCGATTGCGAAAGATCTCCCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGca                                                                                                    <  2:334821/220‑1 (MQ=255)
                                                                                                                                                      caGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTg                                                                                 >  1:704272/1‑235 (MQ=255)
                                                                                                                                                         ttAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTC‑GATGATCCTGTCCGCACAtcgtc                                                                                             >  1:334821/1‑220 (MQ=255)
                                                                                                                                                             tgATTTATGAATAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGc                                                                            <  1:256650/232‑1 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACCAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAg                                                                   >  1:26637/1‑236 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAg                                                                   >  1:748840/1‑236 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGGCCGCACATCGTGGCTTTCCCTTTGGTATCCCGGACGAg                                                                   >  1:392038/1‑236 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGCGGATGAAACCGACGAACGGCCCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGTAACAGGTGTTGGGATCCCGTC‑GATGGGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGActcg                                                                   >  2:633482/1‑233 (MQ=255)
                                                                                                                                                                        aaaGCGCTGGCGAATGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTTCCTGTTCGCAGCTCGTCGCTTTCCCCTTGGTTGCCGtcaggcgtgcg                                                               >  2:385237/1‑225 (MQ=255)
                                                                                                                                                                        aaaGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAg                                                               >  1:584551/1‑235 (MQ=255)
                                                                                                                                                                        aaaGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCCTGACGAGGCAg                                                               >  1:120701/1‑235 (MQ=255)
                                                                                                                                                                                         tctctcCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCggg                                             <  1:633482/236‑1 (MQ=255)
                                                                                                                                                                                         tctctcCTGGCCGGAGCGGCTGGTGCGGGCA‑AAAAGCTGCGGTTTTG‑CTGTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCggg                                             <  1:584353/236‑1 (MQ=255)
                                                                                                                                                                                           tctcCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGaa                                           >  1:270167/1‑236 (MQ=255)
                                                                                                                                                                                           tctcCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑GTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGGGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGTGCCGGGaa                                           >  2:587947/1‑236 (MQ=255)
                                                                                                                                                                                           ggtcCTGGCCGGAGCAGCTGGTGCTGGCA‑AAAGGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCggg                                             <  2:357800/232‑1 (MQ=255)
                                                                                                                                                                                              ccTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATc                                        <  2:270167/236‑1 (MQ=255)
                                                                                                                                                                                              ccTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATc                                        <  1:587947/236‑1 (MQ=255)
                                                                                                                                                                                                  gCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGa                                    >  1:55959/1‑236 (MQ=255)
                                                                                                                                                                                                    cGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGt                                  >  1:478248/1‑236 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCt                                                                                          >  1:466471/1‑177 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCt                                                                                          <  2:466471/177‑1 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGc                                                                            >  1:491763/1‑191 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGc                                                                            <  2:491763/191‑1 (MQ=255)
                                                                                                                                                                                                                  cTGGCA‑AAAAGCTGCGGTTTTG‑CTTTTGGCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCt                    <  1:565071/236‑1 (MQ=255)
                                                                                                                                                                                                                          gtagttGCGGTTTTG‑CTTTTGTCGAAATGTCGGTGGATGAACCCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgcg              <  2:478248/230‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTGG‑ATTTTGACGGATTGTCGGTGAATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGAtctc       <  2:118354/234‑1 (MQ=255)
                                                                                                                                                                                                                                      ttG‑CTTTTGTCGAGATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGATCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCAt  <  1:611479/235‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
CCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCA‑AAAAGCTGCGGTTTTG‑ATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGG‑TTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTC‑GATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCAT  >  NC_000913/3749094‑3749557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑