Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,175,151 (A)6→5 intergenic (‑125/‑237) coaA ← / → thrU pantothenate kinase/tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,175,1460A.95.2% 255.0 / 8.1 63intergenic (‑120/‑242)coaA/thrUpantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (33/27);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GTCATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGT  >  NC_000913/4174998‑4175288
                                                                                                                                                    |                                                                                                                                              
gTCATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGa                                                                                                                                                <  7:142200/149‑1 (MQ=255)
    ttAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGaaaaa                                                                                                                                            >  1:116351/1‑149 (MQ=255)
    ttAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGaaaaa                                                                                                                                            >  7:129157/1‑149 (MQ=255)
               tctTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAg                                                                                                                                 >  4:87057/1‑148 (MQ=255)
               tctTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAg                                                                                                                                 >  6:56905/1‑148 (MQ=255)
               tctTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAg                                                                                                                                 >  8:290088/1‑148 (MQ=255)
                     ataCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCt                                                                                                                           <  1:36248/148‑1 (MQ=255)
                     ataCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGAGGGGGAGGAATAATAAGAAATATCTCTGCCAGGAAGCTa                                                                                                                          >  2:317280/1‑149 (MQ=255)
                       aCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATc                                                                                                                        >  7:57586/1‑149 (MQ=255)
                       aCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATc                                                                                                                        >  5:286198/1‑149 (MQ=255)
                        cTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATc                                                                                                                        <  8:129157/148‑1 (MQ=255)
                         tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAACTCTGCCAGGAAGCTATCGt                                                                                                                      >  6:135147/1‑149 (MQ=255)
                         tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGt                                                                                                                      >  1:176286/1‑149 (MQ=255)
                         tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAACCTATctt                                                                                                                      >  6:171773/1‑147 (MQ=255)
                           aTAATCCTGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGCATGTGTTAATGGGCAGGATGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGtt                                                                                                                     <  2:147955/148‑1 (MQ=255)
                           aTAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGtt                                                                                                                     >  8:113103/1‑148 (MQ=255)
                             aaGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGa                                                                                                                   >  7:5632/1‑148 (MQ=255)
                                cATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAg                                                                                                               >  2:348951/1‑149 (MQ=255)
                                  tGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGc                                                                                                              >  3:129707/1‑148 (MQ=255)
                                    tttCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAt                                                                                                           <  5:56905/149‑1 (MQ=255)
                                    tttCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAt                                                                                                           <  3:87057/149‑1 (MQ=255)
                                      tCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCATATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAtg                                                                                                          <  4:220460/148‑1 (MQ=255)
                                         ggCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                       >  1:336850/1‑148 (MQ=255)
                                          gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                       >  2:178660/1‑147 (MQ=255)
                                          gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGAAATGTGa                                                                                                       <  5:135147/147‑1 (MQ=255)
                                             gTCATAGCTCTAGCCGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGat                                                                                                  <  2:154317/149‑1 (MQ=255)
                                               cATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGatat                                                                                                >  1:18984/1‑149 (MQ=255)
                                                  aGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGaaaaatctc                                                                                                                                        >  2:162989/1‑103 (MQ=255)
                                                  aGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGaaaaatctc                                                                                                                                        <  1:162989/107‑5 (MQ=255)
                                                         tCGACAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGa                                                                                       >  6:122481/1‑149 (MQ=255)
                                                             cAAATCGATATTTGCCATTGCGCGTAACAGTTCAGTTGAATGGGTTAATGGGCAGGAGGTTAACACCAGAGGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaaa                                                                                  <  1:49606/149‑1 (MQ=255)
                                                              aaaGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaaaa                                                                                 <  7:160828/149‑1 (MQ=255)
                                                                     tatTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTAt                                                                          <  7:113103/149‑1 (MQ=255)
                                                                      atTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                                                           <  7:79731/131‑1 (MQ=255)
                                                                      atTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                                                           >  8:79731/1‑131 (MQ=255)
                                                                       tttGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATg                                                                         <  5:144415/148‑1 (MQ=255)
                                                                            cATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAg                                                                   <  2:18984/149‑1 (MQ=255)
                                                                                    gTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgc                                                            <  5:122481/149‑1 (MQ=255)
                                                                                    gTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                                                           >  8:159612/1‑149 (MQ=255)
                                                                                     tCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                                                           >  6:297246/1‑148 (MQ=255)
                                                                                       acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGtt                                                                      >  2:146613/1‑135 (MQ=255)
                                                                                       acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGtt                                                                      <  1:146613/135‑1 (MQ=255)
                                                                                       acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAt                                                        >  2:248936/1‑149 (MQ=255)
                                                                                        caGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAtt                                                       >  1:90655/1‑149 (MQ=255)
                                                                                         aGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttg                                                      <  2:336850/149‑1 (MQ=255)
                                                                                           ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCACTGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGGAGGAGCGCAttg                                                      <  3:229536/147‑1 (MQ=255)
                                                                                                  gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAg                                                                                                    >  3:106637/1‑94 (MQ=255)
                                                                                                  gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAg                                                                                                    <  4:106637/94‑1 (MQ=255)
                                                                                                            aTGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGa                                   <  1:248936/149‑1 (MQ=255)
                                                                                                               ggCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAATGAATTATGGTTTAGCAGGAGCGCATTTTTGAGCTCAATGATGTTTaaa                                 >  6:220681/1‑148 (MQ=255)
                                                                                                                  aggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgt                             <  1:178660/149‑1 (MQ=255)
                                                                                                                  aggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgt                             >  5:137332/1‑149 (MQ=255)
                                                                                                                    gaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTGTGCCAGGTAGCTATAGTTGAAAAGCAATGTGACAGATATTCCATCTGATTAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtgt                           <  5:220681/149‑1 (MQ=255)
                                                                                                                    gaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtgt                           >  3:163636/1‑149 (MQ=255)
                                                                                                                      gggTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGGTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGc                         <  7:290088/149‑1 (MQ=255)
                                                                                                                        gTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCa                                                                 <  6:265115/107‑1 (MQ=255)
                                                                                                                        gTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCa                                                                 >  5:265115/1‑107 (MQ=255)
                                                                                                                        gTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCt                        >  6:334420/1‑148 (MQ=255)
                                                                                                                           acacCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATc                    <  2:116351/149‑1 (MQ=255)
                                                                                                                           acacCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATc                    <  8:57586/149‑1 (MQ=255)
                                                                                                                           acacCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCAGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATc                    <  6:286198/149‑1 (MQ=255)
                                                                                                                               cAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCATTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGcc                >  1:110583/1‑149 (MQ=255)
                                                                                                                                        aggaATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCg       >  4:232384/1‑149 (MQ=255)
                                                                                                                                        aggaATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGAGTGCTAATCTGCCCTCCGTTCg       >  3:223601/1‑149 (MQ=255)
                                                                                                                                          gaATAATAAGAAAAATCTCTGCCAGGAAGCGATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGc     >  1:218727/1‑149 (MQ=255)
                                                                                                                                             taataaGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  >  8:27504/1‑149 (MQ=255)
                                                                                                                                             taataaGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAGGTGACAGATATTCCCTCTGTTAAAAAGAATTATGGGTTAGCAGGAGCGCATTGTTGAGCTCAATGAGGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  >  1:179962/1‑149 (MQ=255)
                                                                                                                                             taataaGAAAAATCTCTGCCAGGAAGCTATCATTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  <  2:110583/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTCATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGT  >  NC_000913/4174998‑4175288

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