Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,385,208 T→A 5.6% V457E (GTG→GAG) ‡ yjeM → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,385,2080TA 5.6% 269.2 / 4.1 87V457E (GTG→GAG) ‡yjeMputative transporter
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base A (3/2);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTA  >  NC_000913/4385061‑4385353
                                                                                                                                                   |                                                                                                                                                 
tCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCtgtg                                                                                                                                                  >  8:146968/1‑149 (MQ=255)
  gcgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGt                                                                                                                                                <  8:138865/149‑1 (MQ=255)
   cgcTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  7:213152/149‑1 (MQ=255)
      tGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTg                                                                                                                                              >  4:73523/1‑147 (MQ=255)
          tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                        >  7:175813/1‑149 (MQ=255)
          tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCg                                                                                                                                        >  7:84642/1‑149 (MQ=255)
          tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGACGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaa                                                                                                                                            <  2:143086/145‑1 (MQ=255)
          tttCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGACGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCACAGGTTGaa                                                                                                                                            >  1:143086/1‑145 (MQ=255)
              ccGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGgag                                                                                                                                     <  1:64576/148‑1 (MQ=255)
                gTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACt                                                                                                                                  >  2:92454/1‑149 (MQ=255)
                   cTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                <  1:171118/148‑1 (MQ=255)
                   cTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTgg                                                                                                                                <  2:27983/148‑1 (MQ=255)
                    tttAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              >  7:71814/1‑149 (MQ=255)
                      tAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGa                                                                                                                              >  6:204674/1‑147 (MQ=255)
                        aaGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGc                                                                                                                          <  8:48199/149‑1 (MQ=255)
                              gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 >  2:12042/1‑136 (MQ=255)
                              gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTg                                                                                                                                 <  1:12042/136‑1 (MQ=255)
                              gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAg                                                                                                                           >  3:130356/1‑142 (MQ=255)
                              gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAg                                                                                                                           <  4:130356/142‑1 (MQ=255)
                              gTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtg                                                                                                                    <  5:55241/149‑1 (MQ=255)
                                 aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCCACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                  >  7:218804/1‑148 (MQ=255)
                                 aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                  >  1:207937/1‑148 (MQ=255)
                                 aGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATtgtg                                                                                                                  >  8:30423/1‑148 (MQ=255)
                                  ggATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGGTGTCGTACTGGTGGTGACATTTGCGAACGTCTGCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGga                                                                                                                <  4:151397/149‑1 (MQ=255)
                                   gATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgat                                                                                                               <  8:2018/149‑1 (MQ=255)
                                     tctcGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATAATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatga                                                                                                             >  2:132498/1‑149 (MQ=255)
                                      ctcGACAGACCGTTTGTGATTTTCAAAACGCATTAGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCAATCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGgatgat                                                                                                            >  7:81974/1‑149 (MQ=255)
                                              aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               <  2:86112/106‑1 (MQ=255)
                                              aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGtt                                                                                                                                               >  1:86112/1‑106 (MQ=255)
                                              aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                    >  7:4115/1‑149 (MQ=255)
                                              aCCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGcc                                                                                                    >  4:134138/1‑149 (MQ=255)
                                               ccGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGccc                                                                                                   >  8:140551/1‑149 (MQ=255)
                                                    tgtgATTTTAATAACGCATTTGTCGGCAATTATTGCGACAGTTGTTGTCGTACTTGTGGTGACATTTGGGAACGTCTTCACCATCATTCATCCTGTGGTTGAAGCCTGAGACTGGGACAGCACTTTGTGTATGATTGGCGTCCCTGtc                                                                                               <  8:81974/148‑1 (MQ=255)
                                                      tgATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGGTGTCGTTCTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttc                                                                                            <  1:237748/149‑1 (MQ=255)
                                                       gATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGga                                                                                                                                      <  5:152583/106‑1 (MQ=255)
                                                       gATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGga                                                                                                                                      >  6:152583/1‑106 (MQ=255)
                                                          tttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtcttcttct                                                                                        >  1:127320/1‑149 (MQ=255)
                                                           ttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtct                                                                                              >  7:131538/1‑142 (MQ=255)
                                                           ttCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGtct                                                                                              <  8:131538/142‑1 (MQ=255)
                                                            tCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTtctc                                                                                       <  8:134450/148‑1 (MQ=255)
                                                                        tGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGAt                                                                           <  3:73523/148‑1 (MQ=255)
                                                                            ggCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCCCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGa                                                                      >  2:44434/1‑149 (MQ=255)
                                                                             gCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAt                                                                     >  4:101878/1‑149 (MQ=255)
                                                                              cAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGAtt                                                                    <  8:4115/149‑1 (MQ=255)
                                                                                aTGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTa                                                                  >  4:166833/1‑149 (MQ=255)
                                                                                 tGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTAc                                                                 <  6:191429/149‑1 (MQ=255)
                                                                                   aTTGCGACAGTGGTTGTCGTTCTGGTGGTGACATTTGCGAACGTCTTCTCCATGATTGAACCTGTGGTTGAAGCCGGAGTCTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTTGCGATGGCGATTTAcct                                                               <  5:197561/149‑2 (MQ=255)
                                                                                   aTTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCa                                                               <  7:146968/149‑1 (MQ=255)
                                                                                    ttGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAg                                                              <  4:113048/149‑1 (MQ=255)
                                                                                       cGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGaa                                                            <  8:84642/148‑1 (MQ=255)
                                                                                          cAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAt                                                        <  3:134138/149‑1 (MQ=255)
                                                                                            gTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATAATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTg                                                      <  1:132498/149‑1 (MQ=255)
                                                                                             tGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTg                                                      <  2:207937/148‑1 (MQ=255)
                                                                                               gTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    <  1:98983/148‑1 (MQ=255)
                                                                                                   tCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCaga                                               <  1:92454/149‑1 (MQ=255)
                                                                                                    cGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCagag                                              >  6:235189/1‑149 (MQ=255)
                                                                                                     gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgc                                                     <  4:8408/141‑1 (MQ=255)
                                                                                                     gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgc                                                     >  3:8408/1‑141 (MQ=255)
                                                                                                     gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGt                                             >  6:115357/1‑149 (MQ=255)
                                                                                                     gTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGt                                             <  4:164605/149‑1 (MQ=255)
                                                                                                       aCTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgcagca                                                 <  4:69183/143‑1 (MQ=255)
                                                                                                       aCTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgcagca                                                 >  3:69183/1‑143 (MQ=255)
                                                                                                         tggtggTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGc                                          >  3:56443/1‑148 (MQ=255)
                                                                                                                gACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGTCAAAAAATc                                  >  4:34351/1‑149 (MQ=255)
                                                                                                                gACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATc                                  >  6:159197/1‑149 (MQ=255)
                                                                                                                    tttGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGGCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCa                              >  8:63196/1‑149 (MQ=255)
                                                                                                                    attGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGCGACAGCACAGTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGCTGGCGATTTACCAGAACTATTGCAGCAGAGAGGCAAAAAATCCg                                >  3:50663/2‑147 (MQ=255)
                                                                                                                         gAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTggg                         >  6:28483/1‑149 (MQ=255)
                                                                                                                                   aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    >  8:67904/1‑144 (MQ=255)
                                                                                                                                   aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  7:67904/144‑1 (MQ=255)
                                                                                                                                   aCCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAAt                >  3:233746/1‑148 (MQ=255)
                                                                                                                                    ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAtt                                                       <  1:49048/108‑1 (MQ=255)
                                                                                                                                    ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTAtt                                                       >  2:49048/1‑108 (MQ=255)
                                                                                                                                    ccATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATaa              <  3:166833/149‑1 (MQ=255)
                                                                                                                                     catcatTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             >  8:128758/1‑149 (MQ=255)
                                                                                                                                     catcatTCAACATGTGGTTGAAGCCGGAGACTGGGACAGCAAATTGTGGATGTTTGGCGGCCCTGTCTTGTTCTCTATGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAAc             <  4:50663/149‑1 (MQ=255)
                                                                                                                                        gcttCAACCACAGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  1:147985/137‑1 (MQ=255)
                                                                                                                                        gcttCAACCACAGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    >  2:147985/3‑139 (MQ=255)
                                                                                                                                         cttCAACCACAGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGa         >  3:215337/2‑149 (MQ=255)
                                                                                                                                         atTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGACGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGa         >  1:248510/1‑149 (MQ=255)
                                                                                                                                          ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    >  2:122196/1‑137 (MQ=255)
                                                                                                                                          ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTg                    <  1:122196/137‑1 (MQ=255)
                                                                                                                                          ttCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGAt        <  7:228657/149‑1 (MQ=255)
                                                                                                                                             aaCCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGcc     >  1:71570/1‑149 (MQ=255)
                                                                                                                                                cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGCTGGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa  <  8:218804/149‑1 (MQ=255)
                                                                                                                                                cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    >  6:56541/1‑99 (MQ=255)
                                                                                                                                                cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTgca                                                    <  5:56541/99‑1 (MQ=255)
                                                                                                                                                cTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTa  <  7:30423/149‑1 (MQ=255)
                                                                                                                                                  gtgGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGc      <  3:89080/143‑1 (MQ=255)
                                                                                                                                                  gtgGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGc      >  4:89080/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGCGCTGGCTTTCCCGTTCTTTAAAGCACGTCAGGATCTCGACAGACCGTTTGTGATTTTCAAAACGCATTTGTCGGCAATGATTGCGACAGTGGTTGTCGTACTGGTGGTGACATTTGCGAACGTCTTCACCATCATTCAACCTGTGGTTGAAGCCGGAGACTGGGACAGCACATTGTGGATGATTGGCGGCCCTGTCTTCTTCTCGCTGTTAGCGATGGCGATTTACCAGAACTATTGCAGCAGAGTGGCAAAAAATCCGCAGTGGGCGGTGGAATAACCCGATGCCCTA  >  NC_000913/4385061‑4385353

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