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

CGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGT  >  NC_000913/3992354‑3992891
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                           
cgcgTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATTC‑GCCATCAGCCGCTGGCATATAACCGTTACCTGAATAAACtgggg                                                                                                                                                                                                                                                                         >  2:600349/1‑273 (MQ=255)
cgcgTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                         >  2:82387/1‑275 (MQ=255)
       gaaTAACCTTGGGGTAGCTGCCGGTCCGCAGGATATTGGCGTGCTGACGCGTAAGCTGTATGCCGCGTGTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGAGTTCACCCGATCTCTCGGAACCGAATCTGACCTTGATTTATGTGCCGCCGTGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTggg                                                                                                                                                                                                                                                                  <  4:491376/274‑1 (MQ=255)
       cAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTggg                                                                                                                                                                                                                                                                  <  1:232173/275‑1 (MQ=255)
                 aGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GACATCCGCCGCTGGAATATAACCGTTACCTGAATAACCGGGTGGCGTGGGCATGGTTTa                                                                                                                                                                                                                                                         >  2:183978/1‑274 (MQ=255)
                        ccGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATTTTGAGTCGACCATCG‑CCTATGCGCCGCTGGAATAAAACCTTACCCCGAAAAAANCGGGGNNGGGGGGCCGGGTTCACGGccc                                                                                                                                                                                                                                                  >  3:20528/1‑273 (MQ=255)
                              gTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGGTGAcc                                                                                                                                                                                                                                           >  2:843272/1‑275 (MQ=255)
                                    aGGATATCTGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGGAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGTAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                                                                                                      <  1:366837/274‑1 (MQ=255)
                                     ggATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  1:20528/275‑1 (MQ=255)
                                        tatCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCACATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTCCCTTAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCCGCTGACCCCGCGCAccc                                                                                                                                                                                                                                 >  1:350957/1‑275 (MQ=255)
                                          tCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                                                                                                               <  2:521852/275‑1 (MQ=255)
                                              cGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCATGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTTACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTGCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTg                                                                                                                                                                                                                            <  1:101878/274‑1 (MQ=255)
                                                tGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCTCGCACCCGTTTGta                                                                                                                                                                                                                          >  1:161005/1‑274 (MQ=255)
                                                    catgcgatatcgatATGCCCGGTTTGAAGAGTTACCACGTAAAGTCACGCTGTTAAACCCGTAGATTTCACCCGATCTCTCGGAAGCGAATCTTACGTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGTTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATAtt                                                                                                                                                                                                                      <  4:138958/261‑1 (MQ=255)
                                                      cgcgGAAGCGGTAGGCCACTTGTGCAANNTTACGAGTGAAAGTGACGTTGGTAAACCCTNTGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGGGCCGCCGTGCCGGGCTAACCGTTCAGGTNGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATAATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTa                                                                                                                                                                                                                     <  3:161005/273‑1 (MQ=255)
                                                                 aaTGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCAGGGTTTAACGGCCTGCTGACCTCGCGCCCCCGGTTGTATATTAAAggctggctct                                                                                                                                                                                                         >  4:37214/2‑266 (MQ=255)
                                                                     ccGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTg                                                                                                                                                                                                                            <  2:37215/251‑1 (MQ=255)
                                                                        cgTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAt                                                                                                                                                                                                 >  2:745083/1‑275 (MQ=255)
                                                                        cgTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAt                                                                                                                                                                                                 >  2:793978/1‑275 (MQ=255)
                                                                                             aaGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTg                                                                                                                                                                                                                                                                    >  2:121685/1‑187 (MQ=255)
                                                                                              aGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaa                                                                                                                                                                                                                    >  2:685450/1‑234 (MQ=255)
                                                                                                tGACGCTGGTAAACCCGCAGATTTAACCCGATCGCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGt                                                                                                                                                                          <  2:407769/274‑1 (MQ=255)
                                                                                                 gACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCa                                                                                                                                                                                                                                                                <  4:121682/187‑1 (MQ=255)
                                                                                                  aCGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGt                                                                                                                                                                                                                <  4:685447/234‑1 (MQ=255)
                                                                                                  aCGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTCATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGCTGGTATCTGTATAACCGCGCGCCAAATGTTGATTCGATCATCA‑GCCATCAGCCGCAGGAATATAACCGTTCCCTGAATACACGGGTCGGGTGGGCATGGTTTAACGGCCTGCCGACCCTGCGCACCCGTTTGTAGCTTAAAGGTAACGGCATTGTCGTTTTGCCTACATTGCAAGAAATggacc                                                                                                                                                                        >  4:315460/1‑271 (MQ=255)
                                                                                                               cttggATTTCGCCCGATGTCTCGGAACCGAATTGGACCTTTTTTAATGGGCGGCCGGGCCGGGTTAACCGTTCAGGTTTGTATCTGTATAACCGCGCGCCAAATATTTAGTCGATCATCA‑GCCATAAGCCGCTGGAATATAACCGTTGCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCa                                                                                                                                                           <  4:40630/270‑1 (MQ=255)
                                                                                                                      ttCACCCGACCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGGTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  1:107480/194‑1 (MQ=255)
                                                                                                                         aCCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGc                                                                                                                                                <  2:530012/275‑1 (MQ=255)
                                                                                                                              aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                         <  8:80567/149‑1 (MQ=255)
                                                                                                                              aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGACGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGCTCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                         >  7:62919/1‑149 (MQ=255)
                                                                                                                                   tcGGAAGCGGATATGACCTTTTTTGTTGGGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTTTCTGTGTAACGGGACGCCAAATATTGAGTGGATCCTCA‑GCCATCAGCCGCTGGAATAGAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGGGTCGCACCATTTCCCGCTGCGCTTAcc                                                                                                                                      <  4:82384/275‑1 (MQ=255)
                                                                                                                                         ccGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                                                                                                                                                              <  6:28528/149‑1 (MQ=255)
                                                                                                                                         ccGAATCATACCTTTATTTGGTTTCCGCCCGGCCGGGTTAACCTTTCAGGTTGTGTTCTGTATAACCGCGCGCCAAATATTGAGTCTATCATCG‑GCCATCAGCCGCTGGAATATAACCGTTACCTTAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcacc                                                                                                                                <  4:843269/275‑1 (MQ=255)
                                                                                                                                            aaTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGGCCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGCCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGGGTCGCCCCATTTCCCGATGCGCTTACCTGccccga                                                                                                                              >  4:18703/1‑274 (MQ=255)
                                                                                                                                                 gACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGATGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                                                                                                        <  2:18704/275‑1 (MQ=255)
                                                                                                                                                   caggtattgatGTGCCCCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATATCCGCGCGCCAACTATTGAGTCGATCATCA‑GCCCTCAGCCGCTTGAATATAACCGTTACCGGAGTAAACGGGTGGCGTGAGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAg                                                                                                                      <  2:675641/271‑1 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                                                                                                      >  4:205196/1‑149 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                                                                                                                                          >  3:67375/1‑161 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATNAAAGGTAACGGCCTTGTCGATTTTNCCNAGTTGCa                                                                                                                                                                                   >  4:982045/1‑200 (MQ=255)
                                                                                                                                                                 ccgccgGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCAGGGTTTAACGGCCTGCTGACCTCGCGCACCAGTTCGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGAACCATTTACCGCTGCGCGTACCTGCAAGCACAACGGAGgcggcctaaatct                                                                                                         >  4:675638/1‑264 (MQ=255)
                                                                                                                                                                     cgggccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCAcc                                                                                                                                                                                                                                  <  2:205199/149‑1 (MQ=255)
                                                                                                                                                                     cgggccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATAtt                                                                                                                                                                                                                      <  1:67375/161‑1 (MQ=255)
                                                                                                                                                                     cgggccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGgag                                                                                                                                                                               <  2:982048/200‑1 (MQ=255)
                                                                                                                                                                             gcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                                                                                                                                          <  6:15722/149‑1 (MQ=255)
                                                                                                                                                                                       cAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGaa                                                                                                                       >  2:529633/1‑238 (MQ=255)
                                                                                                                                                                                           ttGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgcg                                                                                                                   <  4:529630/238‑1 (MQ=255)
                                                                                                                                                                                            agaaacagtcaacaCCGCGCGCAAATNNTTGAGTCG‑TCGTAACGCCATAGGCCGATGGNATATACCCGGTGCCTGAAAAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTNTATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGAt                                                                                <  3:152880/259‑1 (MQ=255)
                                                                                                                                                                                              gTATTTGTAAAACCGCGCGCTAACGATTGNNNCGATCATCA‑GCCGTAGGCCTCTTGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGTATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCg                                                                           <  3:350957/275‑1 (MQ=255)
                                                                                                                                                                                               tATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGt                                                                          <  2:656121/275‑1 (MQ=255)
                                                                                                                                                                                                    gTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                                                                                                                                                          >  1:75654/1‑78 (MQ=255)
                                                                                                                                                                                                     tataACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCAc                                                                                                                                                                                                                                   >  4:984867/1‑116 (MQ=255)
                                                                                                                                                                                                     tagaaCCGCGCGCGAAATATTGTGTCGGTCATCA‑GGCAGCAGCCGCTGGACTATAACCGTTGCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCCCCCGTTTGTATATTAAAGGTAACGGCATTGTCGAGTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCt                                                                    <  2:618023/272‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCg                                                                                                                                                                                                                                                                      <  3:75654/78‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCg                                                                           <  2:66743/265‑1 (MQ=255)
                                                                                                                                                                                                         aCCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                                                                                                               <  2:984870/116‑1 (MQ=255)
                                                                                                                                                                                                               cgcCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACCCCGAAGGCGCTCTACAGCCCGTGTGAAATCCCACATCTGGCTGTTATGCTTAACCTGGAATTTGa                                                           >  4:332909/1‑274 (MQ=255)
                                                                                                                                                                                                               cgcCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTgg                                                                                    >  1:152880/1‑249 (MQ=255)
                                                                                                                                                                                                                              cGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTc                                                                                                                                                                         <  6:3845/149‑1 (MQ=255)
                                                                                                                                                                                                                                atcatcA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                                                                                                                                                          >  1:233421/1‑50 (MQ=255)
                                                                                                                                                                                                                                atcatcA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCt                                                                                                                                                                                                                                          >  4:486588/1‑82 (MQ=255)
                                                                                                                                                                                                                                atcatcA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcg                                                                                                                                                                                                                                        >  4:297839/1‑84 (MQ=255)
                                                                                                                                                                                                                                atcatcA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcg                                                                                                                                                                                                                                        >  2:868156/1‑84 (MQ=255)
                                                                                                                                                                                                                                    tca‑gccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCg                                                                                                                                                                                                                                                                      <  3:233421/50‑1 (MQ=255)
                                                                                                                                                                                                                                    tca‑gccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                                                                                                      <  2:486591/82‑1 (MQ=255)
                                                                                                                                                                                                                                    tca‑gccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  2:297842/84‑1 (MQ=255)
                                                                                                                                                                                                                                    tca‑gccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                    <  4:868153/84‑1 (MQ=255)
                                                                                                                                                                                                                                     ca‑gccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCAttt                                                                                                                                                     >  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)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                           
CGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCA‑GCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGT  >  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: ‑