Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,992,624 T→G 65.6% L491R (CTG→CGG)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,6240TG65.6% 97.1 / 78.6 76L491R (CTG→CGG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base G (21/28);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGT  >  NC_000913/3992354‑3992891
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                           
cgcgTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATTCGCCATCAGCCGCTGGCATATAACCGTTACCTGAATAAACtgggg                                                                                                                                                                                                                                                                         >  2:600349/1‑273 (MQ=255)
cgcgTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                         >  2:82387/1‑275 (MQ=255)
       gaaTAACCTTGGGGTAGCTGCCGGTCCGCAGGATATTGGCGTGCTGACGCGTAAGCTGTATGCCGCGTGTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGAGTTCACCCGATCTCTCGGAACCGAATCTGACCTTGATTTATGTGCCGCCGTGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTggg                                                                                                                                                                                                                                                                  <  4:491376/274‑1 (MQ=255)
       cAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTggg                                                                                                                                                                                                                                                                  <  1:232173/275‑1 (MQ=255)
                 aGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGACATCCGCCGCTGGAATATAACCGTTACCTGAATAACCGGGTGGCGTGGGCATGGTTTa                                                                                                                                                                                                                                                         >  2:183978/1‑274 (MQ=255)
                        ccGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATTTTGAGTCGACCATCGCCTATGCGCCGCTGGAATAAAACCTTACCCCGAAAAAANCGGGGNNGGGGGGCCGGGTTCACGGccc                                                                                                                                                                                                                                                  >  3:20528/1‑273 (MQ=255)
                              gTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGGTGAcc                                                                                                                                                                                                                                           >  2:843272/1‑275 (MQ=255)
                                    aGGATATCTGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGGAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGTAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                                                                                                      <  1:366837/274‑1 (MQ=255)
                                     ggATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  1:20528/275‑1 (MQ=255)
                                        tatCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCACATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTCCCTTAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCCGCTGACCCCGCGCAccc                                                                                                                                                                                                                                 >  1:350957/1‑275 (MQ=255)
                                          tCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                                                                                                               <  2:521852/275‑1 (MQ=255)
                                              cGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCATGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTTACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTGCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTg                                                                                                                                                                                                                            <  1:101878/274‑1 (MQ=255)
                                                tGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCTCGCACCCGTTTGta                                                                                                                                                                                                                          >  1:161005/1‑274 (MQ=255)
                                                    catgcgatatcgatATGCCCGGTTTGAAGAGTTACCACGTAAAGTCACGCTGTTAAACCCGTAGATTTCACCCGATCTCTCGGAAGCGAATCTTACGTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGTTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATAtt                                                                                                                                                                                                                      <  4:138958/261‑1 (MQ=255)
                                                      cgcgGAAGCGGTAGGCCACTTGTGCAANNTTACGAGTGAAAGTGACGTTGGTAAACCCTNTGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGGGCCGCCGTGCCGGGCTAACCGTTCAGGTNGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATAATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTa                                                                                                                                                                                                                     <  3:161005/273‑1 (MQ=255)
                                                                 aaTGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCAGGGTTTAACGGCCTGCTGACCTCGCGCCCCCGGTTGTATATTAAAggctggctct                                                                                                                                                                                                         >  4:37214/2‑266 (MQ=255)
                                                                     ccGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTg                                                                                                                                                                                                                            <  2:37215/251‑1 (MQ=255)
                                                                        cgTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAt                                                                                                                                                                                                 >  2:745083/1‑275 (MQ=255)
                                                                        cgTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAt                                                                                                                                                                                                 >  2:793978/1‑275 (MQ=255)
                                                                                             aaGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTg                                                                                                                                                                                                                                                                    >  2:121685/1‑187 (MQ=255)
                                                                                              aGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaa                                                                                                                                                                                                                    >  2:685450/1‑234 (MQ=255)
                                                                                                tGACGCTGGTAAACCCGCAGATTTAACCCGATCGCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGt                                                                                                                                                                          <  2:407769/274‑1 (MQ=255)
                                                                                                 gACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCa                                                                                                                                                                                                                                                                <  4:121682/187‑1 (MQ=255)
                                                                                                  aCGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGt                                                                                                                                                                                                                <  4:685447/234‑1 (MQ=255)
                                                                                                  aCGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTCATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGCTGGTATCTGTATAACCGCGCGCCAAATGTTGATTCGATCATCAGCCATCAGCCGCAGGAATATAACCGTTCCCTGAATACACGGGTCGGGTGGGCATGGTTTAACGGCCTGCCGACCCTGCGCACCCGTTTGTAGCTTAAAGGTAACGGCATTGTCGTTTTGCCTACATTGCAAGAAATggacc                                                                                                                                                                        >  4:315460/1‑271 (MQ=255)
                                                                                                               cttggATTTCGCCCGATGTCTCGGAACCGAATTGGACCTTTTTTAATGGGCGGCCGGGCCGGGTTAACCGTTCAGGTTTGTATCTGTATAACCGCGCGCCAAATATTTAGTCGATCATCAGCCATAAGCCGCTGGAATATAACCGTTGCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCa                                                                                                                                                           <  4:40630/270‑1 (MQ=255)
                                                                                                                      ttCACCCGACCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGGTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  1:107480/194‑1 (MQ=255)
                                                                                                                         aCCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGc                                                                                                                                                <  2:530012/275‑1 (MQ=255)
                                                                                                                              aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                         <  8:80567/149‑1 (MQ=255)
                                                                                                                              aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGACGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGCTCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                         >  7:62919/1‑149 (MQ=255)
                                                                                                                                   tcGGAAGCGGATATGACCTTTTTTGTTGGGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTTTCTGTGTAACGGGACGCCAAATATTGAGTGGATCCTCAGCCATCAGCCGCTGGAATAGAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGGGTCGCACCATTTCCCGCTGCGCTTAcc                                                                                                                                      <  4:82384/275‑1 (MQ=255)
                                                                                                                                         ccGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                                                                                                                                                              <  6:28528/149‑1 (MQ=255)
                                                                                                                                         ccGAATCATACCTTTATTTGGTTTCCGCCCGGCCGGGTTAACCTTTCAGGTTGTGTTCTGTATAACCGCGCGCCAAATATTGAGTCTATCATCGGCCATCAGCCGCTGGAATATAACCGTTACCTTAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcacc                                                                                                                                <  4:843269/275‑1 (MQ=255)
                                                                                                                                            aaTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGGCCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGCCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGGGTCGCCCCATTTCCCGATGCGCTTACCTGccccga                                                                                                                              >  4:18703/1‑274 (MQ=255)
                                                                                                                                                 gACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGATGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                                                                                                        <  2:18704/275‑1 (MQ=255)
                                                                                                                                                   caggtattgatGTGCCCCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATATCCGCGCGCCAACTATTGAGTCGATCATCAGCCCTCAGCCGCTTGAATATAACCGTTACCGGAGTAAACGGGTGGCGTGAGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAg                                                                                                                      <  2:675641/271‑1 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                                                                                                      >  4:205196/1‑149 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                                                                                                                                          >  3:67375/1‑161 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATNAAAGGTAACGGCCTTGTCGATTTTNCCNAGTTGCa                                                                                                                                                                                   >  4:982045/1‑200 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCAGGGTTTAACGGCCTGCTGACCTCGCGCACCAGTTCGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGAACCATTTACCGCTGCGCGTACCTGCAAGCACAACGGAGgcggcctaaatct                                                                                                         >  4:675638/1‑264 (MQ=255)
                                                                                                                                                                     cgggccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCAcc                                                                                                                                                                                                                                  <  2:205199/149‑1 (MQ=255)
                                                                                                                                                                     cgggccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATAtt                                                                                                                                                                                                                      <  1:67375/161‑1 (MQ=255)
                                                                                                                                                                     cgggccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGgag                                                                                                                                                                               <  2:982048/200‑1 (MQ=255)
                                                                                                                                                                             gcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                                                                                                                                          <  6:15722/149‑1 (MQ=255)
                                                                                                                                                                                       cAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGaa                                                                                                                       >  2:529633/1‑238 (MQ=255)
                                                                                                                                                                                           ttGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgcg                                                                                                                   <  4:529630/238‑1 (MQ=255)
                                                                                                                                                                                            agaaacagtcaacaCCGCGCGCAAATNNTTGAGTCGTCGTAACGCCATAGGCCGATGGNATATACCCGGTGCCTGAAAAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTNTATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGAt                                                                                <  3:152880/259‑1 (MQ=255)
                                                                                                                                                                                              gTATTTGTAAAACCGCGCGCTAACGATTGNNNCGATCATCAGCCGTAGGCCTCTTGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGTATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCg                                                                           <  3:350957/275‑1 (MQ=255)
                                                                                                                                                                                               tATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGt                                                                          <  2:656121/275‑1 (MQ=255)
                                                                                                                                                                                                    gTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                                                                                                                                                          >  1:75654/1‑78 (MQ=255)
                                                                                                                                                                                                     tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCAc                                                                                                                                                                                                                                   >  4:984867/1‑116 (MQ=255)
                                                                                                                                                                                                     tagaaCCGCGCGCGAAATATTGTGTCGGTCATCAGGCAGCAGCCGCTGGACTATAACCGTTGCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCCCCCGTTTGTATATTAAAGGTAACGGCATTGTCGAGTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCt                                                                    <  2:618023/272‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCg                                                                                                                                                                                                                                                                      <  3:75654/78‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCg                                                                           <  2:66743/265‑1 (MQ=255)
                                                                                                                                                                                                         aCCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                                                                                                               <  2:984870/116‑1 (MQ=255)
                                                                                                                                                                                                               cgcCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACCCCGAAGGCGCTCTACAGCCCGTGTGAAATCCCACATCTGGCTGTTATGCTTAACCTGGAATTTGa                                                           >  4:332909/1‑274 (MQ=255)
                                                                                                                                                                                                               cgcCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTgg                                                                                    >  1:152880/1‑249 (MQ=255)
                                                                                                                                                                                                                              cGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTc                                                                                                                                                                         <  6:3845/149‑1 (MQ=255)
                                                                                                                                                                                                                                atcatcAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                                                                                                                                                          >  1:233421/1‑50 (MQ=255)
                                                                                                                                                                                                                                atcatcAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCt                                                                                                                                                                                                                                          >  4:486588/1‑82 (MQ=255)
                                                                                                                                                                                                                                atcatcAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcg                                                                                                                                                                                                                                        >  4:297839/1‑84 (MQ=255)
                                                                                                                                                                                                                                atcatcAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcg                                                                                                                                                                                                                                        >  2:868156/1‑84 (MQ=255)
                                                                                                                                                                                                                                    tcagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCg                                                                                                                                                                                                                                                                      <  3:233421/50‑1 (MQ=255)
                                                                                                                                                                                                                                    tcagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                                                                                                      <  2:486591/82‑1 (MQ=255)
                                                                                                                                                                                                                                    tcagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  2:297842/84‑1 (MQ=255)
                                                                                                                                                                                                                                    tcagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  4:868153/84‑1 (MQ=255)
                                                                                                                                                                                                                                     cagccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCAttt                                                                                                                                                     >  4:286578/1‑162 (MQ=255)
                                                                                                                                                                                                                                          catcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCg                                                                                                                                                 <  2:286581/162‑1 (MQ=255)
                                                                                                                                                                                                                                            tcagccGCTGGAATATAACCGTTACATGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTTGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCTCCATCTGGCGATTATCGTTAACCTGGAATATGGCCAGACCCCAGCGTTCCGCAATCAggag                               >  2:386647/1‑272 (MQ=255)
                                                                                                                                                                                                                                              agccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                                                                                                                                                   <  5:27090/92‑1 (MQ=255)
                                                                                                                                                                                                                                              agccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                                                                                                                                                   >  6:27090/1‑92 (MQ=255)
                                                                                                                                                                                                                                                       aaTATAACGTTTGCCTGAATAAACGGGGGGGGTGGGCAGGGGTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGAATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTGCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAt                   <  4:793975/275‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGaa                                                                <  2:26744/230‑1 (MQ=255)
                                                                                                                                                                                                                                                        atatAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGc                                                                                                                                                >  7:67752/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                gTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGaa                                                                                                                       >  1:247555/1‑166 (MQ=255)
                                                                                                                                                                                                                                                                   aCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgaca                                                                                                                            >  3:244602/1‑158 (MQ=255)
                                                                                                                                                                                                                                                                   aCCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCCCCATTTCCCGCTGCGCTTACCTGcaccgccccc                                                                                                                          >  3:284896/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                    ccTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgcg                                                                                                                   <  3:247555/166‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGca                                                                                                                                  >  8:67483/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       gAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                                                                                                        <  1:244602/158‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAg                                                                                                                      <  1:284896/160‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aaTAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTCAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCGGCACCGACACCGAAGGCGCTCTCCAGCCCGTGTGAGATCCGCCATCTCGCGATTAACGTGAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGTGGGTGCATTTCGATTTCCGTAAGCTGGgtgt  >  4:591555/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                           
CGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGT  >  NC_000913/3992354‑3992891

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