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)
               tctTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAg                                                                                                                                 >  4:87057/1‑148 (MQ=255)
               tctTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAg                                                                                                                                 >  6:56905/1‑148 (MQ=255)
               tctTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAg                                                                                                                                 >  8:290088/1‑148 (MQ=255)
                     ataCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCt                                                                                                                           <  1:36248/148‑1 (MQ=255)
                     ataCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGAGGGGGAGGAATAATAAG‑AAATATCTCTGCCAGGAAGCTa                                                                                                                          >  2:317280/1‑149 (MQ=255)
                       aCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATc                                                                                                                        >  7:57586/1‑149 (MQ=255)
                       aCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATc                                                                                                                        >  5:286198/1‑149 (MQ=255)
                        cTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATc                                                                                                                        <  8:129157/148‑1 (MQ=255)
                         tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAA‑CTCTGCCAGGAAGCTATCGt                                                                                                                      >  6:135147/1‑149 (MQ=255)
                         tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGt                                                                                                                      >  1:176286/1‑149 (MQ=255)
                         tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAACCTATctt                                                                                                                      >  6:171773/1‑147 (MQ=255)
                           aTAATCCTGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGCATGTGTTAATGGGCAGGATGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGtt                                                                                                                     <  2:147955/148‑1 (MQ=255)
                           aTAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGtt                                                                                                                     >  8:113103/1‑148 (MQ=255)
                             aaGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGa                                                                                                                   >  7:5632/1‑148 (MQ=255)
                                cATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAg                                                                                                               >  2:348951/1‑149 (MQ=255)
                                  tGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGc                                                                                                              >  3:129707/1‑148 (MQ=255)
                                    tttCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAt                                                                                                           <  5:56905/149‑1 (MQ=255)
                                    tttCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAt                                                                                                           <  3:87057/149‑1 (MQ=255)
                                      tCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCATATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAtg                                                                                                          <  4:220460/148‑1 (MQ=255)
                                         ggCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                       >  1:336850/1‑148 (MQ=255)
                                          gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                       >  2:178660/1‑147 (MQ=255)
                                          gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGAAATGTGa                                                                                                       <  5:135147/147‑1 (MQ=255)
                                             gTCATAGCTCTAGCCGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGat                                                                                                  <  2:154317/149‑1 (MQ=255)
                                               cATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGatat                                                                                                >  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)
                                                             cAAATCGATATTTGCCATTGCGCGTAACAGTTCAGTTGAATGGGTTAATGGGCAGGAGGTTAACACCAGAGGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaaa                                                                                  <  1:49606/149‑1 (MQ=255)
                                                              aaaGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaaaa                                                                                 <  7:160828/149‑1 (MQ=255)
                                                                     tatTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTAt                                                                          <  7:113103/149‑1 (MQ=255)
                                                                      atTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                                                           <  7:79731/131‑1 (MQ=255)
                                                                      atTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                                                           >  8:79731/1‑131 (MQ=255)
                                                                       tttGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATg                                                                         <  5:144415/148‑1 (MQ=255)
                                                                            cATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAg                                                                   <  2:18984/149‑1 (MQ=255)
                                                                                    gTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgc                                                            <  5:122481/149‑1 (MQ=255)
                                                                                    gTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                                                           >  8:159612/1‑149 (MQ=255)
                                                                                     tCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                                                           >  6:297246/1‑148 (MQ=255)
                                                                                       acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGtt                                                                      >  2:146613/1‑135 (MQ=255)
                                                                                       acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGtt                                                                      <  1:146613/135‑1 (MQ=255)
                                                                                       acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAt                                                        >  2:248936/1‑149 (MQ=255)
                                                                                        caGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAtt                                                       >  1:90655/1‑149 (MQ=255)
                                                                                         aGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttg                                                      <  2:336850/149‑1 (MQ=255)
                                                                                           ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCACTGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGGAGGAGCGCAttg                                                      <  3:229536/147‑1 (MQ=255)
                                                                                                  gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAg                                                                                                    >  3:106637/1‑94 (MQ=255)
                                                                                                  gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAg                                                                                                    <  4:106637/94‑1 (MQ=255)
                                                                                                            aTGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGa                                   <  1:248936/149‑1 (MQ=255)
                                                                                                               ggCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAATGAATTATGGTTTAGCAGGAGCGCATTTTTGAGCTCAATGATGTTTaaa                                 >  6:220681/1‑148 (MQ=255)
                                                                                                                  aggaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgt                             <  1:178660/149‑1 (MQ=255)
                                                                                                                  aggaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgt                             >  5:137332/1‑149 (MQ=255)
                                                                                                                    gaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTGTGCCAGGTAGCTATAGTTGAAAAGCAATGTGACAGATATTCCATCTGATTAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtgt                           <  5:220681/149‑1 (MQ=255)
                                                                                                                    gaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtgt                           >  3:163636/1‑149 (MQ=255)
                                                                                                                      gggTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGGTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGc                         <  7:290088/149‑1 (MQ=255)
                                                                                                                        gTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCa                                                                 <  6:265115/107‑1 (MQ=255)
                                                                                                                        gTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCa                                                                 >  5:265115/1‑107 (MQ=255)
                                                                                                                        gTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCt                        >  6:334420/1‑148 (MQ=255)
                                                                                                                           acacCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATc                    <  2:116351/149‑1 (MQ=255)
                                                                                                                           acacCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATc                    <  8:57586/149‑1 (MQ=255)
                                                                                                                           acacCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCAGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATc                    <  6:286198/149‑1 (MQ=255)
                                                                                                                               cAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCATTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGcc                >  1:110583/1‑149 (MQ=255)
                                                                                                                                        aggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCg       >  4:232384/1‑149 (MQ=255)
                                                                                                                                        aggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGAGTGCTAATCTGCCCTCCGTTCg       >  3:223601/1‑149 (MQ=255)
                                                                                                                                          gaATAATAAG‑AAAAATCTCTGCCAGGAAGCGATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGc     >  1:218727/1‑149 (MQ=255)
                                                                                                                                             taataaG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  >  8:27504/1‑149 (MQ=255)
                                                                                                                                             taataaG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAGGTGACAGATATTCCCTCTGTTAAAAAGAATTATGGGTTAGCAGGAGCGCATTGTTGAGCTCAATGAGGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  >  1:179962/1‑149 (MQ=255)
                                                                                                                                             taataaG‑AAAAATCTCTGCCAGGAAGCTATCATTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  <  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: ‑