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

CCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCAT  >  NC_000913/3749094‑3749557
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
ccctTGCGGACTTGCCGATGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCtt                                                                                                                                                                                                                                         <  1:288849/235‑1 (MQ=255)
ccccTGCTGACCTGCCGCTGTTGAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCtt                                                                                                                                                                                                                                         <  1:705709/235‑1 (MQ=255)
      ctgACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGGGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAACGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCACAAAGCTGCGGTGTTGCTTTTGTCGa                                                                                                                                                                                                                                  >  1:357800/1‑236 (MQ=255)
          ccTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCCTGCGCCAATCGACGCTATCTGGGGAGGTGCGCGGGCGTAATCCTCAGTTAGGGATTTATGGAAAAGCGCTGGCGAAAGGTCTCTCCCGGCCGGGGCGGCGTGGGCTGGCACAAACGCTGGGGGTTGGCTTTTGCGAAATg                                                                                                                                                                                                                              >  2:470734/1‑236 (MQ=255)
              cggCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCg                                                                                                                                                                                                                           <  2:221909/233‑1 (MQ=255)
              ccGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGt                                                                                                                                                                                                                             <  2:501913/233‑1 (MQ=255)
              ccGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGt                                                                                                                                                                                                                             >  1:501913/1‑233 (MQ=255)
                              tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGGTTTGCTTTTGTCGAAATGTCGGTGGATGAAACcgacga                                                                                                                                                                                                          >  1:60714/1‑236 (MQ=255)
                              tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGGAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGGGCGGCTGGGGCTGGCAAAAAACTGCGGGTTTGCTTTTGTCGTAAAGTCGGTGGGTGGAAacgacgc                                                                                                                                                                                                          >  2:198147/1‑235 (MQ=255)
                              tCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGGGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGAGGCTGGGGCTTGGAAAAAGCGGCGGGTTTGCCTTTTGTTGAAAGTTGGGTGGTGGAATCcgcga                                                                                                                                                                                                           >  2:80613/1‑233 (MQ=255)
                                      gaaGTGTTTACTGCCGGGCGGCCACTGGCAGGCGCGGCAACTCCGGCGCAAGTGGCTGGCGATTTCCATGCGAAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGctct                                                                                                                                                                                                  <  2:220265/235‑1 (MQ=255)
                                      aaaGCGTTTATTGCCGGGCGGGCACGGGCAGGCTCGGCAAATCCGGCGGAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGctct                                                                                                                                                                                                  <  2:99205/236‑1 (MQ=255)
                                              tagtgcagagcgggcACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTc                                                                                                                                                                                            <  2:626593/225‑1 (MQ=255)
                                              tATTGCCGGGCGGGCCCTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTc                                                                                                                                                                                            <  2:296274/234‑1 (MQ=255)
                                                           gcACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGcgcc                                                                                                                                                                              <  1:620358/235‑1 (MQ=255)
                                                             aCTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGca                                                                                                                                                                           >  1:596321/1‑236 (MQ=255)
                                                                           cAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCCAAAAGCTGCGGGTTCTCTTTTGTCGAAATGTCGGTGGGTGAAACCGACGACAGGCTCTCACGTCTTGCCTTGAGCGCCGCGCCAAgcatgttgct                                                                                                                                                             >  2:200597/1‑227 (MQ=255)
                                                                                   gcgcAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGGGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGAGCCGTACAAAGGACTTCGCTGGTTGCcg                                                                                                                                                    >  1:716906/1‑236 (MQ=255)
                                                                                   gcgcAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGGTTTGCTTTTGTCGAAATGTCGGTGGATGAAAACGACGAACGGCCCTCACGGCTTGACTGGGAGGCCGGCCCAAGGGCTTCGCCGGTTGCcg                                                                                                                                                    >  2:72043/1‑236 (MQ=255)
                                                                                            gctggcGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGGGGGTGAAACCGACGGAAGGGCCTTCCGGCTTGACTGGGGTGCCGCACAAAGGACTTGCCTGAGTGCcgccgtggttg                                                                                                                                           >  2:505107/1‑228 (MQ=255)
                                                                                                gctaTTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGaaa                                                                                                                                        <  1:113937/232‑1 (MQ=255)
                                                                                                           gcgcAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGTGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGa                                                                                                                            <  1:505107/236‑1 (MQ=255)
                                                                                                                 aTCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGt                                                                                                                      >  1:234327/1‑236 (MQ=255)
                                                                                                                         tATCTGGGGAGGTTCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGAtgtg                                                                                                              <  2:60714/235‑1 (MQ=255)
                                                                                                                         tATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGAtgtg                                                                                                              <  1:569735/235‑1 (MQ=255)
                                                                                                                         tATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGAtgtg                                                                                                              <  1:207628/235‑1 (MQ=255)
                                                                                                                           tCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCtctc                                                                                                                                                                                                 <  1:241520/152‑1 (MQ=255)
                                                                                                                           tCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCtctc                                                                                                                                                                                                 >  2:241520/1‑152 (MQ=255)
                                                                                                                             gggggAGGTGCGCGTGCGTAATCATCAGGTAGGGAGTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGAtgtg                                                                                                              <  1:200597/230‑1 (MQ=255)
                                                                                                                                 gAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGCTGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGc                                                                                                     >  1:235383/1‑236 (MQ=255)
                                                                                                                                  aGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAgcg                                                                                                                                                                                <  2:553337/162‑1 (MQ=255)
                                                                                                                                  aGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAgcg                                                                                                                                                                                >  1:553337/1‑162 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCCAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACCGCTCTCAAGTCTTGACTGGAGCGCCTCACAAAGGACTTCGCTGGCTTCCCGCATGACCGGAACAGGGGTTTGGATCCCATCACATGTCCTGTCCGcaca                                                                                                  >  2:338194/1‑235 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTTACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGAACGAAACAGGTGTTGGGATCCCGTGCACGTGCCTGTCCGtacac                                                                                                 >  2:634604/1‑235 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCCCAAAGGACTTCGCTGGTTGCCGCGATGATCGGAAACGGTGTTCGGATCCCGTCGATGTGCCTGTCCGCACAt                                                                                                 >  2:256650/1‑236 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATcc                                                                                                                         >  1:1432/1‑213 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATcc                                                                                                                         <  2:1432/213‑1 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCTTCGATGGGCCTGTCCGcaca                                                                                                  >  1:698704/1‑235 (MQ=255)
                                                                                                                                     tGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACGTGAGCGCCGCCCAAAGGACTTCGCTGGGTGCCGCGATGGTCGAAACAGGGGTTGGGATCCCGTCGATGGTCCTggtcgcgca                                                                                                  >  2:643202/1‑227 (MQ=255)
                                                                                                                                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGc                                                                                                             <  1:356456/222‑1 (MQ=255)
                                                                                                                                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGc                                                                                                             >  2:356456/1‑222 (MQ=255)
                                                                                                                                       cgcgTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACAtc                                                                                                <  1:219743/235‑1 (MQ=255)
                                                                                                                                         cgtgcgtAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACAtcgtc                                                                                             >  1:421086/1‑236 (MQ=255)
                                                                                                                                            gcgtAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGc                                                                                           <  2:596321/235‑1 (MQ=255)
                                                                                                                                              gtAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACAtc                                                                                                >  1:118354/1‑228 (MQ=255)
                                                                                                                                               tAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACAtc                                                                                                >  1:736811/1‑227 (MQ=255)
                                                                                                                                                  tcatcaGTTAGGGATTTATGAAAAAGCGATTGCGAAAGATCTCCCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGca                                                                                                    <  2:334821/220‑1 (MQ=255)
                                                                                                                                                      caGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTg                                                                                 >  1:704272/1‑235 (MQ=255)
                                                                                                                                                         ttAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGATCCTGTCCGCACAtcgtc                                                                                             >  1:334821/1‑220 (MQ=255)
                                                                                                                                                             tgATTTATGAATAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGc                                                                            <  1:256650/232‑1 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACCAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAg                                                                   >  1:26637/1‑236 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAg                                                                   >  1:748840/1‑236 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGGCCGCACATCGTGGCTTTCCCTTTGGTATCCCGGACGAg                                                                   >  1:392038/1‑236 (MQ=255)
                                                                                                                                                                   aTGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGCGGATGAAACCGACGAACGGCCCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGTAACAGGTGTTGGGATCCCGTCGATGGGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGActcg                                                                   >  2:633482/1‑233 (MQ=255)
                                                                                                                                                                        aaaGCGCTGGCGAATGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTTCCTGTTCGCAGCTCGTCGCTTTCCCCTTGGTTGCCGtcaggcgtgcg                                                               >  2:385237/1‑225 (MQ=255)
                                                                                                                                                                        aaaGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAg                                                               >  1:584551/1‑235 (MQ=255)
                                                                                                                                                                        aaaGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGGGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCCTGACGAGGCAg                                                               >  1:120701/1‑235 (MQ=255)
                                                                                                                                                                                         tctctcCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCggg                                             <  1:633482/236‑1 (MQ=255)
                                                                                                                                                                                         tctctcCTGGCCGGAGCGGCTGGTGCGGGCAAAAAGCTGCGGTTTTGCTGTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCggg                                             <  1:584353/236‑1 (MQ=255)
                                                                                                                                                                                           tctcCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGaa                                           >  1:270167/1‑236 (MQ=255)
                                                                                                                                                                                           tctcCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGGTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGGGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGTGCCGGGaa                                           >  2:587947/1‑236 (MQ=255)
                                                                                                                                                                                           ggtcCTGGCCGGAGCAGCTGGTGCTGGCAAAAGGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCggg                                             <  2:357800/232‑1 (MQ=255)
                                                                                                                                                                                              ccTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATc                                        <  2:270167/236‑1 (MQ=255)
                                                                                                                                                                                              ccTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATc                                        <  1:587947/236‑1 (MQ=255)
                                                                                                                                                                                                  gCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGa                                    >  1:55959/1‑236 (MQ=255)
                                                                                                                                                                                                    cGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGt                                  >  1:478248/1‑236 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCt                                                                                          >  1:466471/1‑177 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCt                                                                                          <  2:466471/177‑1 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGc                                                                            >  1:491763/1‑191 (MQ=255)
                                                                                                                                                                                                       aGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGc                                                                            <  2:491763/191‑1 (MQ=255)
                                                                                                                                                                                                                  cTGGCAAAAAGCTGCGGTTTTGCTTTTGGCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCt                    <  1:565071/236‑1 (MQ=255)
                                                                                                                                                                                                                          gtagttGCGGTTTTGCTTTTGTCGAAATGTCGGTGGATGAACCCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGgcgcg              <  2:478248/230‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTGGATTTTGACGGATTGTCGGTGAATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGAtctc       <  2:118354/234‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGCTTTTGTCGAGATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGATCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCAt  <  1:611479/235‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
CCCCTGCTGACCTGCCGCTGTTTAAAGATATCCGCGTGAAAGCGTTTATTGCCGGGCGGGCACTGGCAGGCGCGGCAAATCCGGCGCAAGTGGCTGGCGATTTCCATGCGCAAATCGACGCTATCTGGGGAGGTGCGCGTGCGTAATCATCAGTTAGGGATTTATGAAAAAGCGCTGGCGAAAGATCTCTCCTGGCCGGAGCGGCTGGTGCTGGCAAAAAGCTGCGGTTTTGATTTTGTCGAAATGTCGGTGGATGAAACCGACGAACGGCTCTCACGTCTTGACTGGAGCGCCGCACAAAGGACTTCGCTGGTTGCCGCGATGATCGAAACAGGTGTTGGGATCCCGTCGATGTGCCTGTCCGCACATCGTCGCTTTCCCTTTGGTAGCCGTGACGAGGCAGTGCGCGAACGGGCGCGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCAT  >  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: