Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,335 T→A 5.5% 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,3350TA 5.5% 968.7 / 31.3 270I61I (ATT→ATA) ‡cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (130/125);  new base A (4/11);  total (134/136)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-04

Only 100 of 275 total aligned reads displayed.

CGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGT  >  NC_000913/3991072‑3991598
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                       
                 cgcgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCCCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCCACAGTCTACTGCCGACATTGTTGCACTACAACCATCCGGCCATGCCGGGTTACCTTTAAGGTAACGGTCCCAAAGGCATTTCCCTTTACGCGCCTGATGAAAcgaca                                                                                                                                                                                                                                              >  4:632318/1‑270 (MQ=255)
                   cgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAAGGTTCCCAAAGGCATTTGCCGTTACACGCCTGCTGAAACTCAACGc                                                                                                                                                                                                                                           >  2:95250/1‑275 (MQ=255)
                               gATGTTGGGGGAACCACAGTCATGACGGGTAGCAAATCAGGGGATACGTCTTGTACCTTGATATTGCGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                                <  4:373936/274‑1 (MQ=255)
                                 tgttgtCCGCATCACGGTACTGACGTGTACAACAGCAGGCGTTACCTCTTGTTCCTCATTATGTAGACTCCGAAACNGGTACTGGTGTCCATAAATCAATTGCGGGAGGATCGCTCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                                                                                                                                              <  3:165598/274‑1 (MQ=255)
                                    ataggagacagcaGTAATGCAGGGTAGCAAATCATGCGATACGTCTTGTACATCTAGATTGAGAATCTGAAACATAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACAGTGTTTCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                                                                                                                                           <  4:510794/263‑1 (MQ=255)
                                                        cGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAtt                                                                                                                                                                                                                                                                         <  4:27725/208‑1 (MQ=255)
                                                          ggTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTACTGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTTATGAAACTCAACGCCACTACCTGAAc                                                                                                                                                                                                                               >  2:510796/1‑248 (MQ=255)
                                                                       gCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATACATCAATTGCGTGTGGATCGCGCGCTTACTGCTATGGGGCCTGCCTTCCAACCGGTCTCCAGTCTACTGCCGACATTGTTGCACTATCACCCTCCGCTAATGCCGGGTTACCTTTCTGGTACCATTCCCAACGGCATTTGCCTTTACACGCCTGATGAAACTCAAAGCCACTACCTGAACGAGCt                                                                                                                                                                                                                          >  2:456066/1‑240 (MQ=255)
                                                                                             atTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTgg                                                                                                                                                                                                           <  2:272817/233‑1 (MQ=255)
                                                                                              ttggtaCTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACt                                                                                                                                                                <  2:632320/270‑1 (MQ=255)
                                                                                               tGAGACGCTTAAACAGAGACTGGATGCCATAAAGCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                               <  1:219997/275‑1 (MQ=255)
                                                                                                   aCTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACTTTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGCTGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTTTTGAAACTCAACGACACTACCTTAACGACCTCTAACTGTAACGTGGAATGCCAATCCAGGCTNNNCCCAAACTGGAGCATCCCATTaacggggt                                                                                                                                                           >  3:227702/1‑268 (MQ=255)
                                                                                                    ctctGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCACTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                                                                                                                                                                                                                                    >  4:138268/1‑201 (MQ=255)
                                                                                                          aaCAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACAcc                                                                                                                                                    <  2:427513/275‑1 (MQ=255)
                                                                                                                   tGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAtt                                                                                                                                                                                                                                                                         >  8:6844/1‑149 (MQ=255)
                                                                                                                    ggATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAtt                                                                                                                                                                                                                                                                         <  5:599908/148‑1 (MQ=255)
                                                                                                                           aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                                  >  6:419376/1‑148 (MQ=255)
                                                                                                                           aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                                  >  6:854163/1‑148 (MQ=255)
                                                                                                                           aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                                                                                                                                                                                                                                  <  5:702000/148‑1 (MQ=255)
                                                                                                                           aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                                                                                                                                                 <  7:187875/149‑2 (MQ=255)
                                                                                                                            tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                                                                                                                                                 >  5:47930/1‑148 (MQ=255)
                                                                                                                            tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgg                                                                                                                                                                                                                                                                <  8:552398/149‑3 (MQ=255)
                                                                                                                            tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                                                                                                                                                 <  6:141772/148‑2 (MQ=255)
                                                                                                                                  aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCt                                                                                                                                                                                                                                                                    >  8:188015/1‑139 (MQ=255)
                                                                                                                                  aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGcc                                                                                                                                                                                                                                                           <  5:292565/148‑1 (MQ=255)
                                                                                                                                  aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCt                                                                                                                                                                                                                                                          >  8:302105/1‑149 (MQ=255)
                                                                                                                                   aTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                                                                                                                                         >  7:758115/1‑149 (MQ=255)
                                                                                                                                   aTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                                                                                                                                         >  7:680341/1‑149 (MQ=255)
                                                                                                                                    ttGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                                                                                                                                         <  5:622990/148‑1 (MQ=255)
                                                                                                                                    ttGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAt                                                                                                                                                                    >  2:602372/1‑233 (MQ=255)
                                                                                                                                         tgtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAACgtt                                                                                                                     >  2:551734/1‑275 (MQ=255)
                                                                                                                                             gATCGCGCGCTTGCTGCTATGGGGCCTACATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                                                                                                                                                 <  5:560251/131‑2 (MQ=255)
                                                                                                                                                  cgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATTTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgt                                                                                                                                                             >  2:119285/1‑226 (MQ=255)
                                                                                                                                                  cgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTg                                                                                                            >  1:73695/1‑275 (MQ=255)
                                                                                                                                                   gggcgtgaggccttGGGGCATTCTTTCCAAAAGGGCGACAGGTGAGTGCCGAAATTTGTGCACAATCACCGGCCGTTAATGCCGGGTTACCTTGATTGTAGAGTTGCCAAAGGCATTGGGCTTGACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACCGGGGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGa                                                                                                           <  4:710549/262‑1 (MQ=255)
                                                                                                                                                   gcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCa                                                                                                                                                                                                                                         >  7:898349/1‑149 (MQ=255)
                                                                                                                                                   gcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                                                                                                                                                 <  8:644517/125‑2 (MQ=255)
                                                                                                                                                    cgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGa                                                                                                           >  2:65659/1‑274 (MQ=255)
                                                                                                                                                          ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGACGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTTCCTTTACACGCCTGATGAAACTCAACGCCACTACCTg                                                                                                                                                                                                                                  >  7:587588/1‑149 (MQ=255)
                                                                                                                                                           tgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                                                                                                                                        >  8:768511/1‑126 (MQ=255)
                                                                                                                                                             ctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                                <  5:100950/148‑1 (MQ=255)
                                                                                                                                                                       tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                                                                                                                                      >  7:9884/1‑148 (MQ=255)
                                                                                                                                                                       tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                                                                                                                                                     >  8:44427/1‑149 (MQ=255)
                                                                                                                                                                        gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTAGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                                                                                                                                                                                                                               <  8:585626/138‑1 (MQ=255)
                                                                                                                                                                         cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                                                                                                                                                                                                                                    <  7:253210/132‑1 (MQ=255)
                                                                                                                                                                         cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt                                                                                                                                                                                                                    >  5:383838/1‑148 (MQ=255)
                                                                                                                                                                         cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAtt                                                                                                                                                                   >  2:343027/1‑197 (MQ=255)
                                                                                                                                                                                aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                                                                                                                                                                                                                      >  5:197752/1‑91 (MQ=255)
                                                                                                                                                                                aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                             <  5:643708/148‑1 (MQ=255)
                                                                                                                                                                                  aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                                                                                                                                                                                                                                <  4:790522/127‑1 (MQ=255)
                                                                                                                                                                                   ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa                                                                                                                                                                                                          <  6:400209/148‑1 (MQ=255)
                                                                                                                                                                                      cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                                                                                                                                                      <  8:903062/149‑1 (MQ=255)
                                                                                                                                                                                             cTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAc                                                                                                                                                                                                >  6:55786/1‑148 (MQ=255)
                                                                                                                                                                                                  gCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                             >  7:481832/1‑130 (MQ=255)
                                                                                                                                                                                                  gCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                                                                                                                                             <  8:481832/130‑1 (MQ=255)
                                                                                                                                                                                                  gCCGACATNNTTGCTCTATCACCATCCGCTAATNCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                                                                                                                                     <  3:54952/138‑1 (MQ=255)
                                                                                                                                                                                                      aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATc                                                                                                                                                                                                               >  2:176275/1‑124 (MQ=255)
                                                                                                                                                                                                      aCATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                                                                                                                                     >  6:648508/1‑134 (MQ=255)
                                                                                                                                                                                                        aTTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGCTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAAt                                                      <  4:481433/275‑1 (MQ=255)
                                                                                                                                                                                                          tgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGACATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaa                                                                                                                                                                                   <  5:861926/148‑1 (MQ=255)
                                                                                                                                                                                                           gttgCACCATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCGGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGc                                                   <  2:497963/275‑1 (MQ=255)
                                                                                                                                                                                                               cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGa                                                                                                                                                                             <  7:302105/149‑1 (MQ=255)
                                                                                                                                                                                                                 cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                                                                                                                                          <  6:227454/102‑1 (MQ=255)
                                                                                                                                                                                                                  tATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                                                                                                                                          >  8:385641/1‑149 (MQ=255)
                                                                                                                                                                                                                     cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTAt                                                                                                                                                                                                                <  1:128300/108‑1 (MQ=255)
                                                                                                                                                                                                                     cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTcc                                                                                                                                                                       <  7:862857/149‑1 (MQ=255)
                                                                                                                                                                                                                           ccGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                                  >  5:497689/1‑148 (MQ=255)
                                                                                                                                                                                                                            cGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTACTTAATCCTGGCTCGTAACTCAAGAGCCACAATTGCTACAACCTAAGTGTaac                                  >  4:598036/1‑273 (MQ=255)
                                                                                                                                                                                                                             gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                               <  7:771570/149‑1 (MQ=255)
                                                                                                                                                                                                                             gCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                               >  8:787424/1‑149 (MQ=255)
                                                                                                                                                                                                                              cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                               <  6:580403/148‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAtt                                                                                                                                                                   >  4:22016/1‑140 (MQ=255)
                                                                                                                                                                                                                                    ccGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGaa                                                                                                                                                                                                         >  1:48751/1‑100 (MQ=255)
                                                                                                                                                                                                                                    ccGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                                  <  4:442346/139‑1 (MQ=255)
                                                                                                                                                                                                                                     aacgTTACCATCAAGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc                                                                                                                                                                                                                     >  7:326622/4‑87 (MQ=255)
                                                                                                                                                                                                                                      gggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                                                                                                                                               <  2:22016/140‑1 (MQ=255)
                                                                                                                                                                                                                                      gggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATac                                                                                                                                                       >  6:655881/1‑148 (MQ=255)
                                                                                                                                                                                                                                      gggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATac                                                                                                                                                       <  5:794333/148‑1 (MQ=255)
                                                                                                                                                                                                                                         ttACCATCAAGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCa                                                                                                                                                   >  7:545149/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCAt                                                                                                                                                  <  7:381758/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ccTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATgg                                                                                                                                                <  7:363432/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ccTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATgg                                                                                                                                                >  8:355549/1‑149 (MQ=255)
                                                                                                                                                                                                                                             cTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                                                                                                                                               >  7:121205/1‑149 (MQ=255)
                                                                                                                                                                                                                                             cTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                                                                                                                                               >  8:372110/1‑149 (MQ=255)
                                                                                                                                                                                                                                              ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgt                                                                                                                                                             <  5:15079/134‑1 (MQ=255)
                                                                                                                                                                                                                                              ttGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgc                                                                                                                                              >  7:163470/1‑149 (MQ=255)
                                                                                                                                                                                                                                               tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTAAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTgg                                                                                                                                                              >  5:907009/1‑132 (MQ=255)
                                                                                                                                                                                                                                               tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTAAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTgg                                                                                                                                                              <  6:907009/132‑1 (MQ=255)
                                                                                                                                                                                                                                                    gTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAcc                                                                                                                                        <  8:121205/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTa                                                                                                                                                                  >  6:163048/1‑121 (MQ=255)
                                                                                                                                                                                                                                                      aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCt                                                                                                                                       >  5:657779/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAgg                                                                                                                                                                                             <  6:768897/92‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGCACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                                                   >  7:642768/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGGAAATGTAGCCTGCTGGAAAACTGGGCCGGCTCGCAGGgggt  >  1:109736/1‑275 (MQ=255)
                                                                                                                                                                                                                                                             ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa                                                                                                                               <  8:881137/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              cAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa                                                                                                                               <  5:773645/148‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCa                                                                                                                                                                      >  2:747481/1‑108 (MQ=255)
                                                                                                                                                                                                                                                               aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAg                                                                                                                              >  5:163318/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                aaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAgg                                                                                                                             >  6:297611/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAg                                                                                                                       <  6:814714/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                       
CGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGAAAACTGGGCCGCCTCGCTGGGTGT  >  NC_000913/3991072‑3991598

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: