Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,335 T→A 11.1% I61I (ATT→ATA) ‡ cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,3350TA11.1% 804.0 / 78.1 252I61I (ATT→ATA) ‡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 = 3.03e-03

Only 100 of 267 total aligned reads displayed.

GACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCT  >  NC_000913/3991066‑3991609
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                  
                                              ggAATCAAAGGTATTAAGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGGGTGGATCGCTCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                                                                                                                                                 <  4:75017/275‑1 (MQ=255)
                                                         cATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTTTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                                      <  2:202842/275‑1 (MQ=255)
                                                         cATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACGCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                                      <  4:200772/275‑1 (MQ=255)
                                                                         tCAGGCGATACGTTTTGTACCTCTATAGTGAGACTATGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGCTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGAt                                                                                                                                                                                                      <  2:402800/275‑1 (MQ=255)
                                                                                                   attgggagcaggAAACAGAGACTGGATGCCATAAATCAATTGCGGGTGGATGGCGCGCTGGCTGCTAGGGTGCCTGCATGCCAACAGGTATACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCGTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGACGAAAGGTGAGCTTCCAATTa                                                                                                                                                                             <  4:37775/263‑1 (MQ=255)
                                                                                                    ttGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACt                                                                                                                                                                           <  2:326459/275‑1 (MQ=255)
                                                                                                    ttGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACt                                                                                                                                                                           >  2:495125/1‑275 (MQ=255)
                                                                                                     tGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                                          <  2:128476/275‑1 (MQ=255)
                                                                                                                         tGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTTCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagc                                                                                                                                                      <  1:107499/275‑1 (MQ=255)
                                                                                                                                  tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                                             <  6:1116349/147‑1 (MQ=255)
                                                                                                                                      tCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACg                                                                                                                                                                                                                                                                        <  7:493702/148‑1 (MQ=255)
                                                                                                                                       cAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGc                                                                                                                                                                                                                                                                       <  6:480464/148‑1 (MQ=255)
                                                                                                                                          ttGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                                                                                                                                                    <  6:325337/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)
                                                                                                                                                                  gctATGGGGCCGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                                           >  6:414161/1‑148 (MQ=255)
                                                                                                                                                                    tATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCttt                                                                                                                                                                                                                                                                             <  7:148831/113‑1 (MQ=255)
                                                                                                                                                                       ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                                             <  6:183894/110‑1 (MQ=255)
                                                                                                                                                                       ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                                             >  5:183894/1‑110 (MQ=255)
                                                                                                                                                                           ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                                                                                                                                                     >  2:546463/1‑146 (MQ=255)
                                                                                                                                                                            cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGa                                                                                                                                                                                                                                            >  8:1022580/1‑138 (MQ=255)
                                                                                                                                                                              gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                                                                                                                                                                >  6:402944/1‑148 (MQ=255)
                                                                                                                                                                               cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                                                                                                                                                   >  6:825947/1‑112 (MQ=255)
                                                                                                                                                                               cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                                                                                                                                                 <  4:546462/146‑1 (MQ=255)
                                                                                                                                                                                 ttCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                                             <  2:337452/196‑1 (MQ=255)
                                                                                                                                                                                   ccAACAGGTCTAAAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                                                                                                                                                                                                                                               <  7:340912/96‑1 (MQ=255)
                                                                                                                                                                                     aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCg                                                                                                                                                                                                                         <  5:402944/148‑1 (MQ=255)
                                                                                                                                                                                      aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                                        <  6:30837/148‑1 (MQ=255)
                                                                                                                                                                                       cAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                                        <  7:659944/147‑1 (MQ=255)
                                                                                                                                                                                         ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTgg                                                                                       <  2:49001/274‑1 (MQ=255)
                                                                                                                                                                                            cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATg                                                                                                                                                                                                                  >  5:1072227/1‑148 (MQ=255)
                                                                                                                                                                                             tACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                                                                                                                                                                 <  7:784373/148‑1 (MQ=255)
                                                                                                                                                                                              aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                                     <  8:1088170/143‑1 (MQ=255)
                                                                                                                                                                                              aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                                                                                                                                                            >  7:59861/1‑87 (MQ=255)
                                                                                                                                                                                                  tCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTac                                                                                                                                                                                                                                                                           <  6:913344/85‑1 (MQ=255)
                                                                                                                                                                                                   cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTAt                                                                                                                                                                                                                           >  8:359119/1‑132 (MQ=255)
                                                                                                                                                                                                    tACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGaa                                                                           >  2:79391/1‑275 (MQ=255)
                                                                                                                                                                                                     aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCttt                                                                                                                                                                                                                                                                             <  6:1069717/80‑1 (MQ=255)
                                                                                                                                                                                                      cTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                                                                                                                                                         >  6:6041/1‑115 (MQ=255)
                                                                                                                                                                                                          cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                                                                                                                                                    <  6:239726/84‑1 (MQ=255)
                                                                                                                                                                                                           gACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccg                                                                                                                                                                                                   <  6:1047664/148‑1 (MQ=255)
                                                                                                                                                                                                            aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                                                                                                                                                            <  8:1074017/74‑1 (MQ=255)
                                                                                                                                                                                                            aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGt                                                                                                                                                                                                             >  6:36174/1‑137 (MQ=255)
                                                                                                                                                                                                            aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgc                                                                                                                                                                                                  >  7:1084214/1‑148 (MQ=255)
                                                                                                                                                                                                              aTTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                                                                                                                                                                                                                                               <  7:571122/101‑1 (MQ=255)
                                                                                                                                                                                                              aTTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                                                                                                                                                                                                                                               >  8:571122/1‑101 (MQ=255)
                                                                                                                                                                                                                   ttaccTATCACCTTTCTTTAATGCCGGTNTNCCTTTATTGTAACGTTCACAAAGGAATTTGCCTTTAACCGCCTGGNGAAACTCAACGCCACTGCCCGAACGGGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAANGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTa                                                            <  3:42014/271‑1 (MQ=255)
                                                                                                                                                                                                                     taCTTTCACCATCAGCTAAGGCCGGGTTACCGTGTTGGAAACGTACCCAAAGGAATTTCCCTTTACATGCCATATGAAACTCAAAGCCACTACCTGAACGAGCTTGAACTGAATCGTGGAAGGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCAAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGAGCAAAGTTGTGCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGGGCGCCAATTGCTAc                                                           <  4:289637/273‑1 (MQ=255)
                                                                                                                                                                                                                     cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                                           >  6:747458/1‑98 (MQ=255)
                                                                                                                                                                                                                     cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGTTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                                         <  6:955994/148‑1 (MQ=255)
                                                                                                                                                                                                                     cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCa                                                                                                                                                                                                               >  7:254462/1‑126 (MQ=255)
                                                                                                                                                                                                                      aCTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCtt                                                                                                                                                                                                                                                                              >  6:990257/1‑62 (MQ=255)
                                                                                                                                                                                                                       cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                                       <  5:414161/148‑1 (MQ=255)
                                                                                                                                                                                                                        tATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAgg                                                                                                                                                                                           <  8:1048289/143‑1 (MQ=255)
                                                                                                                                                                                                                        tATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAgg                                                                                                                                                                                           >  7:1048289/1‑143 (MQ=255)
                                                                                                                                                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGATCGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                                     >  5:1115799/1‑116 (MQ=255)
                                                                                                                                                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                                     <  6:1115799/116‑1 (MQ=255)
                                                                                                                                                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                                     >  7:918549/1‑148 (MQ=255)
                                                                                                                                                                                                                          tCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTa                                                     >  2:64419/1‑275 (MQ=255)
                                                                                                                                                                                                                             ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAGTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCa                                                                                                                                                                                 <  8:918549/148‑1 (MQ=255)
                                                                                                                                                                                                                               aTCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGt                                                                                                                                                                                                                             >  8:258305/1‑102 (MQ=255)
                                                                                                                                                                                                                                tCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                                                                                                                                                                                                                                             >  1:9476/1‑69 (MQ=255)
                                                                                                                                                                                                                                 ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                                             <  6:50136/148‑1 (MQ=255)
                                                                                                                                                                                                                                   gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGa                                                                                                                                                                                               >  5:273933/1‑128 (MQ=255)
                                                                                                                                                                                                                                    cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                                     <  5:885534/105‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCt                                                                                     <  1:57626/229‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAAAGGATCCGCCGAAAGGTGAGCTTCCAATTCATGGTTTATACACCATGGGCAGCACCTCGTCCGTACGGCAAAGTTGTTCCTCTGACATGGATATCTGGGTCTGTCATCAATCCTGGCCCGATAGCCAACAGCCCCAATTGCTACAACGTAACTGTCGCCTGTCGGaaa                                   >  4:172643/1‑274 (MQ=255)
                                                                                                                                                                                                                                               ttACCATCAAGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACa                                                                                                                                                                                                          <  5:439965/105‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                                           >  5:742251/1‑67 (MQ=255)
                                                                                                                                                                                                                                                    ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCGGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                                                                                                                                                          >  8:698836/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgtg                                                                                                                                                                       >  6:699681/1‑134 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgtg                                                                                                                                                                       <  5:699681/134‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATaca                                                                                                                                                                 >  6:733765/1‑140 (MQ=255)
                                                                                                                                                                                                                                                     tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATaca                                                                                                                                                                 <  5:733765/140‑1 (MQ=255)
                                                                                                                                                                                                                                                           tAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCtt                                                                                                                                                                                                                                    >  5:1119961/1‑67 (MQ=255)
                                                                                                                                                                                                                                                             aCGTTCCCAAAGGCATGTGCCTTTACACGACTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGGCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAACTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCt                   <  1:111597/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa                                                                                                                                                                                                                                  <  5:321888/66‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                                       <  7:328767/109‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                                       >  8:328767/1‑109 (MQ=255)
                                                                                                                                                                                                                                                               gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                                      <  5:231821/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                cgtgtaaaGGCAAATGCCTTTGCACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACt                                                                                                                                                                           <  5:799829/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgtg                                                                                                                                                                       <  8:175147/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gtaaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaaa                                                                                                                                                                                             >  7:1055819/3‑98 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                                                                                                                                                         <  5:97679/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                                          <  6:534758/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                    taaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                                          >  5:534758/2‑116 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                                                                                                                                                                 <  7:1140603/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                                                                                                                                                                 >  8:1140603/1‑77 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                                     >  8:607205/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGt                                                                                                                                           >  5:573419/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                                                                                                                                                         >  8:511572/1‑52 (MQ=255)
                                                                                                                                                                                                                                                                     aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAgg                                                                                                                                                                                           <  5:489981/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAATGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagc                                                                                                                                                      >  8:765001/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                      aaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGtat                                                                                                                                                                    >  5:831242/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                        ggCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                                       <  7:799180/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                                                                                                                                                       >  8:799180/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                            tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAc                                                                                                                                                                                                           <  5:942055/75‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGa                                                                                                                                                                                        <  5:353388/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTGCCAATTACTGGTGTATACACCATGGGCAGCACCTCGGGCGAAGGGCAAAg                                                                                                                                  >  6:866473/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                  
GACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGTGGAAGTCAGCT  >  NC_000913/3991066‑3991609

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