Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,028,546 G→A A156V (GCC→GTC)  ygfS ← putative electron transport protein YgfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,028,5460GA98.1% 173.2 / ‑4.6 52A156V (GCC→GTC) ygfSputative electron transport protein YgfS
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (29/22);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

ATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGC  >  NC_000913/3028303‑3028785
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
aTGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGAcaacaa                                                                                                                                                                                                                                           <  1:382178/250‑1 (MQ=255)
                 tgCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACt                                                                                                                                                                                                                          <  1:351514/250‑1 (MQ=255)
                  gCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCg                                                                                                                                                                                                                        <  2:488353/251‑1 (MQ=255)
                      tgtTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTg                                                                                                                                                                                                                    >  1:209037/1‑251 (MQ=255)
                       gtTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCGTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTg                                                                                                                                                                                                                    <  1:135396/250‑1 (MQ=255)
                        ttCATCGTTCTGACTTGTAATATCGTTACCGGAATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTACTAGATTCTCCGCGCCTTTCAGCGCGGCGATTTCTTTCAGAGAGGATTCTCACCGCTGACACCAATACGCTGCTGGCGGACTCGCCTg                                                                                                                                                                                                                    >  2:133036/1‑249 (MQ=255)
                         tCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAg                                                                                                                                                                                                                  <  2:559210/250‑1 (MQ=255)
                                        gTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAg                                                                                                                                                                                                   <  1:113330/250‑1 (MQ=255)
                                                     cGGTATTATGCTGGGCTTTGTGACACTGGTCGTCGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCt                                                                                                                                                                                     >  2:495076/1‑251 (MQ=255)
                                                                    cTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCa                                                                                                                                                                           >  1:301683/1‑246 (MQ=255)
                                                                    cTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATGTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTGTCAGCGCGGAGATTTCTTTCAGAGAGGATTCGCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCa                                                                                                                                                                           <  2:301683/246‑1 (MQ=255)
                                                                          gACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGc                                                                                                                                                                <  1:25061/251‑1 (MQ=255)
                                                                            cacTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAg                                                                                                                                                              <  2:730250/251‑1 (MQ=255)
                                                                                 gTCGTAGGCCGCGTCTTTGCACGCGAATGTCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATGTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCGGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTgg                                                                                                                                                         <  2:668063/251‑1 (MQ=255)
                                                                                    gTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGccc                                                                                                                                                      >  1:233810/1‑251 (MQ=255)
                                                                                       ggCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCctc                                                                                                                                                   >  1:15312/1‑251 (MQ=255)
                                                                                                   gCACGCGACTGGCAAGAGCTGCATATTGGCACGGTCATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGACTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTAGCCTGAGTTCTCTTTCAGTAAGCAACTGAAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTc                                                                                                                                       >  2:368926/1‑251 (MQ=255)
                                                                                                        cgAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTcac                                                                                                                                     >  1:91270/1‑248 (MQ=255)
                                                                                                              gCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGc                                                                                                                            >  2:579536/1‑251 (MQ=255)
                                                                                                                aaaaGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTAACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCAt                                                                                                                          <  1:71604/251‑1 (MQ=255)
                                                                                                                               gTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTGCACAATTtgttgt                                                                                                           <  2:145881/251‑1 (MQ=255)
                                                                                                                                  cGGTGATCATTTCTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCg                                                                                                         <  2:565527/250‑1 (MQ=255)
                                                                                                                                     tGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCCCCGCTGACACCAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGt                                                                                                      >  2:730584/1‑250 (MQ=255)
                                                                                                                                     tGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGCGGTTGGGAAATATCATAGATTCTACGCGCATTTCAGCGCGGAGCTTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGCGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCACCGCAGGCTGGCCCCTCTTCACGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGt                                                                                                      >  2:363827/1‑250 (MQ=255)
                                                                                                                                         cATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATccc                                                                                                 <  1:568191/251‑1 (MQ=255)
                                                                                                                                          aTTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAACGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGAGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCCTTTCACAATTTGTTGTCGGGTATcccc                                                                                                >  1:41445/1‑251 (MQ=255)
                                                                                                                                            tACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATccccc                                                                                               <  2:730563/250‑1 (MQ=255)
                                                                                                                                                         gTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCATGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAAt                                                                                 >  2:333690/1‑251 (MQ=255)
                                                                                                                                                            gCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGt                                                                              >  1:338410/1‑251 (MQ=255)
                                                                                                                                                                            gtTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACggg                                                              <  2:233810/251‑1 (MQ=255)
                                                                                                                                                                               aggCAATCTAATAGATTCTCCGCGACTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTAGTTGTCGGGTATCCCCCGGCAACGACTGCATGGTGATCATTCCAACCGGgcac                                                           >  2:24208/2‑251 (MQ=255)
                                                                                                                                                                                 gCAATGTAATAAATTATCAGCGCCTTTCAACGCGGAGATTTCTTTAAGAGCGGATTCTCACCGCTGGCAACAATCCGCTGCTGGCGGACTCGCCTGAGTTATCTTTCAGTACGCAACTGCAGCGCCTGTGTTGGGCATGATTCAACGCAGGCTGGCCACTCTTCCCGCTGTTCACAGAGATCGCATTTCAACATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGTTCATTCCAAACGGgaacgc                                                         >  2:304161/1‑251 (MQ=255)
                                                                                                                                                                                          tAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGc                                                 >  1:73793/1‑250 (MQ=255)
                                                                                                                                                                                           aGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGATAGTATTCTCGCCGCTGACAACAATACGCTGCGGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGCGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCaa                                               <  2:403157/251‑1 (MQ=255)
                                                                                                                                                                                           aGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCa                                                >  1:637804/1‑250 (MQ=255)
                                                                                                                                                                                                ctcCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTg                                          >  1:462922/1‑251 (MQ=255)
                                                                                                                                                                                                ctcCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTg                                          >  2:481353/1‑251 (MQ=255)
                                                                                                                                                                                                  ccGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAActct                                           <  1:189796/248‑1 (MQ=255)
                                                                                                                                                                                                    gcgcCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACATAGATCGCAttt                                                                                                                        <  2:675513/169‑1 (MQ=255)
                                                                                                                                                                                                    gcgcCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCAttt                                                                                                                        >  1:675513/1‑169 (MQ=255)
                                                                                                                                                                                                      gcCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACACTACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGcc                                    >  2:38090/1‑251 (MQ=255)
                                                                                                                                                                                                       ccTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGcc                                    >  1:489188/1‑250 (MQ=255)
                                                                                                                                                                                                       ccTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGcc                                    >  1:405782/1‑250 (MQ=255)
                                                                                                                                                                                                             aGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACaa                             <  2:213538/251‑1 (MQ=255)
                                                                                                                                                                                                              gcgcGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGAATCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGCGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAAc                            >  1:42507/1‑251 (MQ=255)
                                                                                                                                                                                                                  ggAGATTTCTTTCAGAGAGGATTCTCACCGCGGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATGCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACTGTCGCGCACTGACGCAACTTTGACAGCCAATACAACGGGc                        <  2:276326/251‑1 (MQ=255)
                                                                                                                                                                                                                      aTTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGa                     >  1:702499/1‑250 (MQ=255)
                                                                                                                                                                                                                             tcAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTc              <  1:224667/250‑1 (MQ=255)
                                                                                                                                                                                                                              cAGAGAGGCTTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCCCAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCCCTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTg            >  2:619891/1‑251 (MQ=255)
                                                                                                                                                                                                                               agagagGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGc           >  1:410377/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gagGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTACGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGTCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACTGGCAGAATTGGTCTGcacc        >  1:2911/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tctcACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGc  >  1:142296/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
ATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGC  >  NC_000913/3028303‑3028785

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