Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,853 C→T 6.3% P471P (CCG→CCA)  eptC ← LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,8530CT 6.3% 193.4 / 2.7 63P471P (CCG→CCA) eptCLPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base T (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

TCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGTTATATGAC  >  NC_000913/4148705‑4148994
                                                                                                                                                    |                                                                                                                                             
tCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgc                                                                                                                                               >  2:116986/1‑149 (MQ=255)
  aGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcgg                                                                                                                                             >  2:183384/1‑149 (MQ=255)
   gTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTGATTACGCCCCTGTGTTTGGTgcggc                                                                                                                                            >  6:45391/1‑149 (MQ=255)
   gTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggc                                                                                                                                            >  3:222576/1‑149 (MQ=255)
       cGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggcggc                                                                                                                                         <  3:2304/148‑1 (MQ=255)
          cAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCgtgt                                                                                                                                     >  6:44370/1‑149 (MQ=255)
             gCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCg                                                                                                                                   >  1:162423/1‑148 (MQ=255)
             gCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCg                                                                                                                                   >  8:92449/1‑148 (MQ=255)
             gCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGt                                                                                                                                  <  1:183384/149‑1 (MQ=255)
                            aaCATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGGTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCg                                                                                                                    <  3:164323/148‑1 (MQ=255)
                              cATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTgg                                                                                                                 <  3:49140/149‑1 (MQ=255)
                                tCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGt                                                                                                                >  8:42885/1‑148 (MQ=255)
                                tCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTc                                                                                                               <  1:116986/149‑1 (MQ=255)
                                  cTGCGAGAAATCACGGGGACGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAATCCTCTTCACCGTGGt                                                                                                                >  5:79940/1‑146 (MQ=255)
                                      gagaAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAgaga                                                                                                          >  1:77093/1‑148 (MQ=255)
                                      gagaAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGaa                                                                                                         <  2:8793/149‑1 (MQ=255)
                                       agaAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAATCCTCTTCACCGTGGTCAGAGaaa                                                                                                        <  7:706/149‑1 (MQ=255)
                                          aaTCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATaa                                                                                                     <  2:162423/149‑1 (MQ=255)
                                           aTCACGGGGAGGAGGCGCTTGCCATTTTTCTGACGTCCACAGCAGGATCGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATaaa                                                                                                    <  4:173853/149‑1 (MQ=255)
                                              aCGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAAcc                                                                                                  <  7:92449/148‑1 (MQ=255)
                                                      gAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCa                                                                                                              <  4:129034/128‑1 (MQ=255)
                                                      gAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCa                                                                                                              >  3:129034/1‑128 (MQ=255)
                                                             ttGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGt                                                                                  >  7:147592/1‑149 (MQ=255)
                                                               gacATTTTTCTGACGTCCACAGGAGGAACGGAATGGTGTACATATGACTCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACAGTTGTCAGAGAAATAAACCAGGAAACCGTTCGGGtt                                                                                 <  4:181394/146‑1 (MQ=255)
                                                               gCCATTTTTCTGACGTGCACAGCAGGAACGGATTGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCATGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGt                                                                                  <  7:80050/147‑1 (MQ=255)
                                                                  aTTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACACCCCTGAGTTTTGTGTGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTtgct                                                                             >  4:185419/1‑149 (MQ=255)
                                                                     tttCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTtgctgc                                                                           <  5:44370/148‑1 (MQ=255)
                                                                           aCGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTaaa                                                                     <  2:77093/148‑1 (MQ=255)
                                                                           aCGTCCACAGCAGGAACGGAATCGTGTACATATGACGCCTCCGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAg                                                                    >  8:69312/1‑149 (MQ=255)
                                                                           aCGACCACAGCAGGAACGGAATGGTGTACATATGACTCGTCGGGTTGTCCTCATTACGCCCCTGAGATTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGCAATAACCCAGGAAACCGTTCGGGTTTGCTGCTTTaa                                                                      <  5:45391/147‑1 (MQ=255)
                                                                                  cagcatGAACGGAATGGTGTACATATGACGCGTCTGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGTCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTaa                                                             <  3:16731/149‑1 (MQ=255)
                                                                                     cagGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAAt                                                            >  5:201210/1‑147 (MQ=255)
                                                                                     cagGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCa                                                          <  8:147592/149‑1 (MQ=255)
                                                                                        gAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGc                                                       >  3:194885/1‑149 (MQ=255)
                                                                                          aCGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTg                                                     >  7:124324/1‑149 (MQ=255)
                                                                                           cGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTgg                                                    >  8:99099/1‑149 (MQ=255)
                                                                                              aaTGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCaa                                                 >  6:160674/1‑149 (MQ=255)
                                                                                                   tgtACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACAGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACcaca                                            >  1:82284/1‑149 (MQ=255)
                                                                                                    gtACATATGACGCGTCGGGTTGTCCTCATTACGCCCATGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACAt                                           >  8:68100/1‑149 (MQ=255)
                                                                                                    gtACATATGACGCGTCGGGTTGTCCTCATTACGCCACTGAGTTTTGTGCGGAGGCGTGTCGTAAACCTCTTAACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGATGCTTTAAAGTCTTTAATCAGGCTGTCAAACACAt                                           >  5:182788/1‑149 (MQ=255)
                                                                                                        atatGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATc                                       <  5:160674/149‑1 (MQ=255)
                                                                                                         tatGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCg                                      >  5:228298/1‑149 (MQ=255)
                                                                                                          atGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGt                                     >  7:212887/1‑149 (MQ=255)
                                                                                                           tGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGt                                     <  7:99099/148‑1 (MQ=255)
                                                                                                               gcgTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATAc                                <  8:212887/149‑1 (MQ=255)
                                                                                                                    gggTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGt                            <  8:124324/148‑1 (MQ=255)
                                                                                                                     ggTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTa                          <  7:69312/149‑1 (MQ=255)
                                                                                                                          tCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGtt                     >  3:16418/1‑149 (MQ=255)
                                                                                                                            cTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCTCATGATCGTTATACAGGGTAGCGTTg                    <  4:194885/148‑1 (MQ=255)
                                                                                                                             tCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCCCCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGt                   >  3:126268/1‑148 (MQ=255)
                                                                                                                             tCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGt                   >  7:251913/1‑148 (MQ=255)
                                                                                                                             tCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGt                   >  8:34052/1‑148 (MQ=255)
                                                                                                                              cATTACGCCCCTGAGGTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTAAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCa                 <  8:25102/149‑1 (MQ=255)
                                                                                                                               aTTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCCCATGATCGTTATACAGGTTAGCGTTGTCa                 >  5:141599/1‑148 (MQ=255)
                                                                                                                                    gCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCg           >  4:206417/1‑149 (MQ=255)
                                                                                                                                      cccTGAGTTTTGTGTGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGtt         <  3:185419/149‑1 (MQ=255)
                                                                                                                                      cccTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCg           >  7:38985/1‑147 (MQ=255)
                                                                                                                                      cccTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCg           <  8:38985/147‑1 (MQ=255)
                                                                                                                                      cccTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCg           >  7:239032/1‑147 (MQ=255)
                                                                                                                                      cccTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCg           <  8:239032/147‑1 (MQ=255)
                                                                                                                                       ccTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGTta        >  1:236113/1‑149 (MQ=255)
                                                                                                                                        cTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTaaa                                                                     <  6:30499/87‑1 (MQ=255)
                                                                                                                                        cTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTaaa                                                                     >  5:30499/1‑87 (MQ=255)
                                                                                                                                           aGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGTGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGTTatat     <  3:7854/148‑1 (MQ=255)
                                                                                                                                           aGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGTtata      <  6:201210/147‑1 (MQ=255)
                                                                                                                                             ttttGTGTGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCa                 <  8:134179/134‑1 (MQ=255)
                                                                                                                                             ttttGTGTGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCa                 >  7:134179/1‑134 (MQ=255)
                                                                                                                                             ttttGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGTTATATGAc  >  3:211456/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAAACCGTTCGGGTTTGCTGCTTTAAAGTCTTTAATCAGGCTGGCAACCACATGATCGTTATACAGGTTAGCGTTGTCATAATCGTTATATGAC  >  NC_000913/4148705‑4148994

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: ‑