Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,061,235 A→C Y24S (TAT→TCT)  yihM → putative sugar phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,061,2350AC98.9% 326.2 / NA 88Y24S (TAT→TCT) yihMputative sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (48/39);  total (49/39)
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.98e-01

CGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACG  >  NC_000913/4061013‑4061462
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   
cGATATCAGCAAGTATCTTGGCGGTTCAGCGGGGATGCCAGCGATCTTAAATGAAAAGCATAGCGCCGACCAGCAAGGCAATAGTGGCCAAAAAGATATTGAAAATTGGCGATTTGAGGCAGGTGACCTCATAAATAATCTGTTGGTGTTTTATGGCGACAATAAATAACGCAAGAAAGATTCTACAACGTGGCGACACTCTTCCTCTTTACTTACATGCTTCTGTCCTTCATTTAAat                                                                                                                                                                                                                     >  2:793290/1‑239 (MQ=255)
           aGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGa                                                                                                                                                                                                          <  2:507270/239‑1 (MQ=255)
             tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGAAAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGaaa                                                                                                                                                                                                        >  1:626747/1‑239 (MQ=255)
               tCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAaga                                                                                                                                                                                                      <  1:320449/239‑1 (MQ=255)
                         tcttctggTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGGCCAAATAGATATTGAATATTGACGATTTGAGGCAGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctg                                                                                                                                                                                            <  2:682748/233‑1 (MQ=255)
                              gggTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt                                                                                                                                                                                       <  1:672185/239‑1 (MQ=255)
                                                          cATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATc                                                                                                                                                           >  2:366055/1‑239 (MQ=255)
                                                             cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTAc                                                                                                                                                        >  1:744351/1‑239 (MQ=255)
                                                                    aCCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                                                                                                                                                 <  2:158176/223‑1 (MQ=255)
                                                                    aCCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                                                                                                                                                 >  1:158176/1‑223 (MQ=255)
                                                                                 tatTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAata                                                                                                                                                                                                                    <  2:417135/159‑1 (MQ=255)
                                                                                 tatTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAata                                                                                                                                                                                                                    >  1:417135/1‑159 (MQ=255)
                                                                                           aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCg                                                                                                                          >  1:1077661/1‑239 (MQ=255)
                                                                                                tatTGAATATTGGCTATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTCCGTGGCGTCAAAATCCATGTTCTCGATGGCGa                                                                                                                         <  1:793290/235‑1 (MQ=255)
                                                                                                          tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAACATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTg                                                                                                           >  2:1115920/1‑239 (MQ=255)
                                                                                                               aTTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctg                                                                                                                                                                                            <  1:3735/153‑1 (MQ=255)
                                                                                                               aTTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctg                                                                                                                                                                                            >  2:3735/1‑153 (MQ=255)
                                                                                                                 ttGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGa                                                                                                                         >  1:509241/1‑218 (MQ=255)
                                                                                                                 ttGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGa                                                                                                                         <  2:509241/218‑1 (MQ=255)
                                                                                                                               cTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgctg                                                                                                                                                                                         <  1:93835/140‑1 (MQ=255)
                                                                                                                               cTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgctg                                                                                                                                                                                         >  2:93835/1‑140 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGACAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGa                                                                                        >  2:1154096/1‑234 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAg                                                                                                                                                                      <  2:196071/156‑1 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAg                                                                                                                                                                      >  1:196071/1‑156 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTtct                                                                                                               <  1:567075/211‑1 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTtct                                                                                                               >  2:567075/1‑211 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAg                                                                                         >  2:682931/1‑233 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAg                                                                                         <  1:682931/233‑1 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGa                                                                                        >  1:110690/1‑234 (MQ=255)
                                                                                                                                  atcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAActct                                                                                   <  1:1149029/239‑1 (MQ=255)
                                                                                                                                     atTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTg                                                                                <  2:638205/239‑1 (MQ=255)
                                                                                                                                       tAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAActc                                                                                    >  2:587502/1‑233 (MQ=255)
                                                                                                                                         aTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAActc                                                                                    >  2:454053/1‑231 (MQ=255)
                                                                                                                                          tCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGtt                                                                                                                    <  2:863715/198‑1 (MQ=255)
                                                                                                                                          tCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGtt                                                                                                                    >  1:863715/1‑198 (MQ=255)
                                                                                                                                          tCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCttt                                                                           <  1:366055/239‑1 (MQ=255)
                                                                                                                                          tCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCGGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCttt                                                                           <  1:454053/239‑1 (MQ=255)
                                                                                                                                            tGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTTATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTgg                                                                         >  1:954905/1‑239 (MQ=255)
                                                                                                                                              tacgtgcTTTATGGTGACAATAAATAACGCAACAAAGATTCTACAACGTGTCGAGAGTCTTGGTCTTTATTTACATGCTTATGCTTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTCATATGGACGATAAAGCACTCTCTGCCTTTgggg                                                                       >  2:299527/4‑239 (MQ=255)
                                                                                                                                              tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                                                                                                                                                 >  2:3826/1‑149 (MQ=255)
                                                                                                                                              tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                                                                                                                                                                 <  1:3826/149‑1 (MQ=255)
                                                                                                                                              tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTgggg                                                                       >  2:93673/1‑239 (MQ=255)
                                                                                                                                              tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTgggg                                                                       >  2:592904/1‑239 (MQ=255)
                                                                                                                                               aGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                 <  2:175907/212‑1 (MQ=255)
                                                                                                                                               aGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                 >  1:175907/1‑212 (MQ=255)
                                                                                                                                                   gtTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAg                                                                                         >  1:397529/1‑216 (MQ=255)
                                                                                                                                                   gtTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAg                                                                                         <  2:397529/216‑1 (MQ=255)
                                                                                                                                                   gtTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCtt                                                                            >  1:1130810/1‑229 (MQ=255)
                                                                                                                                                   gtTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCtt                                                                            <  2:1130810/229‑1 (MQ=255)
                                                                                                                                                    ttttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGttt                                                                                                                   >  2:250126/1‑189 (MQ=255)
                                                                                                                                                    ttttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGttt                                                                                                                   <  1:250126/189‑1 (MQ=255)
                                                                                                                                                    ttttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATaaa                                                                  >  1:812629/1‑238 (MQ=255)
                                                                                                                                                         tGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCccgc                                                            <  2:812629/239‑1 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCAttt                                                                                                                                                                                                                         >  2:512694/1‑77 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCAttt                                                                                                                                                                                                                         <  1:512694/77‑1 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGa                                                                                                                                                                                                          >  1:119415/1‑92 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGa                                                                                                                                                                                                          <  2:119415/92‑1 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc                                                                                                                                                                                              >  2:829208/1‑104 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc                                                                                                                                                                                              <  1:829208/104‑1 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGccc                                                              >  1:438656/1‑232 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGccc                                                              >  2:804283/1‑232 (MQ=255)
                                                                                                                                                              aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTCTGGGGATAAAGccc                                                              >  1:832226/1‑232 (MQ=255)
                                                                                                                                                                     ataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCtt                                                >  1:648242/1‑239 (MQ=255)
                                                                                                                                                                        aCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGACct                                                                                     >  2:75844/1‑199 (MQ=255)
                                                                                                                                                                        aCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAAct                                                                                     <  1:75844/199‑1 (MQ=255)
                                                                                                                                                                         cGCAAGAAAGATTCTACAACGTTTCGACTCTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCccg                                                             <  2:862822/222‑1 (MQ=255)
                                                                                                                                                                         cGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCccg                                                             >  1:862822/1‑222 (MQ=255)
                                                                                                                                                                          gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                 <  2:138002/185‑1 (MQ=255)
                                                                                                                                                                          gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGa                                                                                                 >  1:138002/1‑185 (MQ=255)
                                                                                                                                                                            aagaaagaTTCTACAACGTGTCGAACCTCTTCCTCTTTATTTTCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTc                                         <  1:596600/239‑1 (MQ=255)
                                                                                                                                                                                   aTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGa                                  >  1:423120/1‑239 (MQ=255)
                                                                                                                                                                                            cGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACTTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTTTGCCTTTGGGGATAATGCCCGCCGTCTGAACCTTGATAGTAATATTGAAACCAgcg                          >  2:966552/1‑238 (MQ=255)
                                                                                                                                                                                             gtgtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAg                            <  1:804283/235‑1 (MQ=255)
                                                                                                                                                                                              tgtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGcc                        >  1:229348/1‑238 (MQ=255)
                                                                                                                                                                                               gtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCt                       <  2:438656/238‑1 (MQ=255)
                                                                                                                                                                                                       tctTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCACGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAg               >  2:961414/1‑238 (MQ=255)
                                                                                                                                                                                                         ttCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAActct                                                                                   <  2:542598/168‑1 (MQ=255)
                                                                                                                                                                                                         ttCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAActct                                                                                   >  1:542598/1‑168 (MQ=255)
                                                                                                                                                                                                          tCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGc                                                                                                                                                                                                              <  2:1047058/44‑1 (MQ=255)
                                                                                                                                                                                                          tCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGc                                                                                                                                                                                                              >  1:1047058/1‑44 (MQ=255)
                                                                                                                                                                                                                 tatttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt                                                                                                                                                                                       >  2:735477/1‑60 (MQ=255)
                                                                                                                                                                                                                 tatttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt                                                                                                                                                                                       <  1:735477/60‑1 (MQ=255)
                                                                                                                                                                                                                   tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTc                                                                                                                                                >  1:948778/1‑97 (MQ=255)
                                                                                                                                                                                                                   tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTc                                                                                                                                                <  2:948778/97‑1 (MQ=255)
                                                                                                                                                                                                                   tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGt                                                                                                                     <  1:942809/124‑1 (MQ=255)
                                                                                                                                                                                                                   tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGt                                                                                                                     >  2:942809/1‑124 (MQ=255)
                                                                                                                                                                                                                   tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGCTAAGGCATCTATTGACg  >  2:813901/1‑239 (MQ=255)
                                                                                                                                                                                                                   tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACg  <  2:423120/239‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   
CGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGTTCTCGATGGCGAGCGTTTTTCTCTTGGTAATATGGACGATAAAGAACTCTCTGCCTTTGGGGATAAAGCCCGCCGTCTGAACCTTGATATTCATATTGAAACCAGCGCCTCAGATAAGGCATCTATTGACG  >  NC_000913/4061013‑4061462

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