Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,334 T→A 11.2% I61N (ATT→AAT) ‡ cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,3340TA11.2% 798.4 / 79.2 252I61N (ATT→AAT) ‡cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (107/116);  new base A (13/15);  total (120/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-05

Only 100 of 262 total aligned reads displayed.

GACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCT  >  NC_000913/3991066‑3991592
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                  
                                              ggAATCAAAGGTATTAAGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGGGTGGATCGCTCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                                                                                                                                <  4:75017/275‑1 (MQ=255)
                                                         cATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACGCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                     <  4:200772/275‑1 (MQ=255)
                                                                         tCAGGCGATACGTTTTGTACCTCTATAGTGAGACTATGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGCTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGAt                                                                                                                                                                                     <  2:402800/275‑1 (MQ=255)
                                                                                      cTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                        <  2:163417/275‑1 (MQ=255)
                                                                                               ctATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGCACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCa                                                                                                                                                                >  2:37776/1‑274 (MQ=255)
                                                                                                  tatTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGTTGAGCTTCCAATTa                                                                                                                                                            >  2:289638/1‑275 (MQ=255)
                                                                                                   attgggagcaggAAACAGAGACTGGATGCCATAAATCAATTGCGGGTGGATGGCGCGCTGGCTGCTAGGGTGCCTGCATGCCAACAGGTATACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCGTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGACGAAAGGTGAGCTTCCAATTa                                                                                                                                                            <  4:37775/263‑1 (MQ=255)
                                                                                                   atgtggaCTCTGAAACCGAGACTGGATGCCATAAATCAATTGCGTGTGGCTCGCGCGCTTGCTGCTAGGGGGCCTGCAATCCAACAGTGCTACAGTCTACTGCCGACATTGTGGCACTATCACCATCCGCGAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACGACCTGAACGAGCTTGAACTGTATCGTGGGATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                            <  1:64406/269‑1 (MQ=255)
                                                                                                    ttGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACt                                                                                                                                                          <  2:326459/275‑1 (MQ=255)
                                                                                                       agaCTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTCCCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCCCTACCTGAACGAGCTTGAACTGTATAGTGACATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgt                                                                                                                                                       >  2:536956/1‑275 (MQ=255)
                                                                                                                         tGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTTCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagc                                                                                                                                     <  1:107499/275‑1 (MQ=255)
                                                                                                                          ggATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAtt                                                                                                                                                                                                                                                                   <  8:367372/148‑1 (MQ=255)
                                                                                                                              gCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGcc                                                                                                                                                                                                                                                               >  7:45389/1‑148 (MQ=255)
                                                                                                                                               tgtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaa                                                                                                                                                                                                                                              <  5:852512/148‑1 (MQ=255)
                                                                                                                                                gtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                                                                                                                                                                                                             >  5:1047664/1‑148 (MQ=255)
                                                                                                                                                gtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAgttg                                                                                                              <  4:495124/275‑1 (MQ=255)
                                                                                                                                                     tggcgggcTTGCTGCTATGGGGCCTGCATTCCAAAAGGTTTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCct                                                                                                         <  4:450163/273‑1 (MQ=255)
                                                                                                                                                       gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACg                                                                                                                                                                                                                                      >  8:784373/1‑148 (MQ=255)
                                                                                                                                                                  gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCttt                                                                                                                                                                                                                                                            >  6:412326/1‑115 (MQ=255)
                                                                                                                                                                       ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                            >  5:183894/1‑110 (MQ=255)
                                                                                                                                                                           ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                                                                                                                                                                                                                                 <  6:521708/117‑1 (MQ=255)
                                                                                                                                                                           ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                                                                                                                                    >  2:546463/1‑146 (MQ=255)
                                                                                                                                                                            cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                          >  5:715159/1‑139 (MQ=255)
                                                                                                                                                                               cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                                                                                                                                  <  5:825947/112‑1 (MQ=255)
                                                                                                                                                                               cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                                                                                                                                  >  6:825947/1‑112 (MQ=255)
                                                                                                                                                                               cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                                                                                                                                <  4:546462/146‑1 (MQ=255)
                                                                                                                                                                                     aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCg                                                                                                                                                                                                        <  5:402944/148‑1 (MQ=255)
                                                                                                                                                                                         ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTgg                                                                      <  2:49001/274‑1 (MQ=255)
                                                                                                                                                                                            cTAAAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATg                                                                                                                                                                                                 >  5:955994/1‑148 (MQ=255)
                                                                                                                                                                                              aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                    <  8:1088170/143‑1 (MQ=255)
                                                                                                                                                                                                   cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                                                                                                                                                                                                                       >  7:1027836/1‑103 (MQ=255)
                                                                                                                                                                                                   cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                                                                                                                                                                                                                       <  8:1027836/103‑1 (MQ=255)
                                                                                                                                                                                                   cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                       >  6:119619/1‑135 (MQ=255)
                                                                                                                                                                                                   cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                       <  5:119619/135‑1 (MQ=255)
                                                                                                                                                                                                    tACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACa                                                                                                                                                                                         <  6:1072227/148‑1 (MQ=255)
                                                                                                                                                                                                     aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                    <  7:853298/136‑1 (MQ=255)
                                                                                                                                                                                                     aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCttt                                                                                                                                                                                                                                                            >  5:1069717/1‑80 (MQ=255)
                                                                                                                                                                                                      cTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                                                                                                                                        <  5:6041/115‑1 (MQ=255)
                                                                                                                                                                                                       tGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAtt                                                                                                                                                                                                                                                                   <  8:1081321/71‑1 (MQ=255)
                                                                                                                                                                                                       tGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAtt                                                                                                                                                                                                                                                                   >  7:1081321/1‑71 (MQ=255)
                                                                                                                                                                                                          cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                                                                                                                                   >  5:239726/1‑84 (MQ=255)
                                                                                                                                                                                                            aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGt                                                                                                                                                                                            <  5:36174/137‑1 (MQ=255)
                                                                                                                                                                                                                   ttaccTATCACCTTTCTTTAATGCCGGTNTNCCTTTATTGTAACGTTCACAAAGGAATTTGCCTTTAACCGCCTGGNGAAACTCAACGCCACTGCCCGAACGGGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAANGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTa                                           <  3:42014/271‑1 (MQ=255)
                                                                                                                                                                                                                   tgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGGGCTTGAACTGTATCGTGGAATGTCAGTACa                                                                                                                                                                                         <  6:1050854/133‑1 (MQ=255)
                                                                                                                                                                                                                    gCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTAc                                          <  2:80218/275‑1 (MQ=255)
                                                                                                                                                                                                                     taCTTTCACCATCAGCTAAGGCCGGGTTACCGTGTTGGAAACGTACCCAAAGGAATTTCCCTTTACATGCCATATGAAACTCAAAGCCACTACCTGAACGAGCTTGAACTGAATCGTGGAAGGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCAAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGAGCAAAGTTGTGCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGGGCGCCAATTGCTAc                                          <  4:289637/273‑1 (MQ=255)
                                                                                                                                                                                                                     cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                          <  5:747458/98‑1 (MQ=255)
                                                                                                                                                                                                                      aCTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                       <  5:898314/116‑1 (MQ=255)
                                                                                                                                                                                                                      aCTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                       >  6:898314/1‑116 (MQ=255)
                                                                                                                                                                                                                       cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                      <  5:414161/148‑1 (MQ=255)
                                                                                                                                                                                                                        tATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAc                                                                                                                                                                                                                               >  5:565798/1‑90 (MQ=255)
                                                                                                                                                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGATCGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                    >  5:1115799/1‑116 (MQ=255)
                                                                                                                                                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                     >  7:951/1‑115 (MQ=255)
                                                                                                                                                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                    <  6:1115799/116‑1 (MQ=255)
                                                                                                                                                                                                                          tCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTa                                    >  2:64419/1‑275 (MQ=255)
                                                                                                                                                                                                                             ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAGTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCa                                                                                                                                                                <  8:918549/148‑1 (MQ=255)
                                                                                                                                                                                                                               aTCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGt                                                                                                                                                                                                            <  7:258305/102‑1 (MQ=255)
                                                                                                                                                                                                                               aTCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGt                                                                                                                                                                                                            >  8:258305/1‑102 (MQ=255)
                                                                                                                                                                                                                               aTCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAt                                                                                                                                                              >  5:50136/1‑148 (MQ=255)
                                                                                                                                                                                                                                tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                                                                                                                                                                                                                            >  1:9476/1‑69 (MQ=255)
                                                                                                                                                                                                                                 ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                            <  6:50136/148‑1 (MQ=255)
                                                                                                                                                                                                                                   gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGa                                                                                                                                                                              >  5:273933/1‑128 (MQ=255)
                                                                                                                                                                                                                                   gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGa                                                                                                                                                                              <  6:273933/128‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAGAGGATcc                                                                                                                                                                                   <  7:1051292/116‑1 (MQ=255)
                                                                                                                                                                                                                                          ccGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                                                                                                                                                   >  8:1051292/1‑116 (MQ=255)
                                                                                                                                                                                                                                           cGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGtata                                                                                                                                                  >  7:598596/1‑148 (MQ=255)
                                                                                                                                                                                                                                             ggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAAAGGATCCGCCGAAAGGTGAGCTTCCAATTCATGGTTTATACACCATGGGCAGCACCTCGTCCGTACGGCAAAGTTGTTCCTCTGACATGGATATCTGGGTCTGTCATCAATCCTGGCCCGATAGCCAACAGCCCCAATTGCTACAACGTAACTGTCGCCTGTCGGaaa                  >  4:172643/1‑274 (MQ=255)
                                                                                                                                                                                                                                                 aCCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                                                                                                                               >  6:471283/1‑97 (MQ=255)
                                                                                                                                                                                                                                                 aCCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                                                                                                                               <  5:471283/97‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                          <  6:742251/67‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCGGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                                                                                                                                         >  8:698836/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                    >  8:334011/1‑120 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                    <  7:334011/120‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgtg                                                                                                                                                      >  6:699681/1‑134 (MQ=255)
                                                                                                                                                                                                                                                           tAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                        >  6:372168/1‑110 (MQ=255)
                                                                                                                                                                                                                                                           tAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAcc                                                                                                                                  <  8:598596/148‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCGTTCCCAAAGGCATGTGCCTTTACACGACTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGGCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAACTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCt  <  1:111597/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                                                                                                                                                                                                                 <  5:321888/66‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                                                                                                                                                                                                                 >  6:321888/1‑66 (MQ=255)
                                                                                                                                                                                                                                                              cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                        <  6:721600/107‑1 (MQ=255)
                                                                                                                                                                                                                                                               gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                     >  6:231821/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACAAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                        >  7:1129716/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                 gtgtaaaGGCAAATGCCTTTGCACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                                             >  7:252925/5‑147 (MQ=255)
                                                                                                                                                                                                                                                                 gtgtaaaGGCAAATGCCTTTGCACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                                             <  8:252925/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gtaaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                                                                                                                                                <  6:335210/76‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gtaaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                            <  7:716489/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                        >  6:97679/2‑101 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                        <  5:97679/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                         >  7:118996/2‑116 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                         <  8:118996/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                    <  7:607205/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                    >  8:607205/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGt                                                                                                                          <  6:573419/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                                                                                                                                        >  8:511572/1‑52 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                                                                                                                                        <  7:511572/52‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgc                                                                                                                                                                                 >  5:737020/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                      aaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGtat                                                                                                                                                   <  6:831242/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                      >  8:799180/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                        ggCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTgg                                                                                                                                                        >  7:794952/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                            tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAc                                                                                                                                                                                          >  6:942055/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                  
GACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCT  >  NC_000913/3991066‑3991592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: