Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,653 A→T 5.6% L133H (CTC→CAC)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,6530AT 5.6% 226.2 / 3.6 71L133H (CTC→CAC) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base A (36/31);  new base T (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGC  >  NC_000913/1964511‑1964799
                                                                                                                                              |                                                                                                                                                  
cTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCAACTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACACACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAAt                                                                                                                                                >  3:64561/1‑147 (MQ=255)
  cATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCACGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGa                                                                                                                                              >  5:243379/1‑147 (MQ=255)
   aTCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGa                                                                                                                                              <  6:156381/146‑1 (MQ=255)
       tCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAaca                                                                                                                                        >  5:170727/1‑148 (MQ=255)
        cACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACaaa                                                                                                                                          <  4:64561/145‑1 (MQ=255)
          ccAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAAcacca                                                                                                                                     <  7:148673/148‑1 (MQ=255)
            aatcaatcCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAAcaccac                                                                                                                                    >  5:186734/1‑147 (MQ=255)
            aatcaatcCGGCGTTAAGGTCGGCGTCAATCGCCATATGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAAcaccac                                                                                                                                    >  1:253798/1‑147 (MQ=255)
             atcaatcCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGATAATGACAAACACCACGa                                                                                                                                  >  5:10282/1‑148 (MQ=255)
             atcaatcCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGa                                                                                                                                  >  7:108118/1‑148 (MQ=255)
                  tcCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGACATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATAcc                                                                                                                              >  7:315213/1‑147 (MQ=255)
                            aGGTCGGCGTCAATCGCCATCTGATTACCCGGCATACCATCGAGAACAACGCGCGCACCAACTTCTGCGATACTCCCGGCACCTTTTGTAATGACCATAAAGTTGATGATCACGAGAATGACAATCACCACGATACCGATAGCGAAAt                                                                                                                   <  4:223546/148‑1 (MQ=255)
                                 cgcctCAACCGCCATCAGCATACCCCGCACACCACCCCGCACACCCCGCCCACCCACATCCGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATCATCACGAGAATGACAAACACCACGCTACCGAAAGCGAAATTGcca                                                                                                              <  8:205958/144‑1 (MQ=255)
                                      cAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAAc                                                                                                         >  2:305040/1‑148 (MQ=255)
                                       aaTCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACg                                                                                                        <  8:26377/148‑1 (MQ=255)
                                        aTCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGa                                                                                                       <  2:253798/148‑1 (MQ=255)
                                           gCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAgg                                                                                                     >  7:145981/1‑147 (MQ=255)
                                              aTCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGt                                                                                                 >  5:93619/1‑148 (MQ=255)
                                                cTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCg                                                                                                                       <  4:177236/124‑1 (MQ=255)
                                                cTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCg                                                                                                                       >  3:177236/1‑124 (MQ=255)
                                                  gCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACAAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGAcc                                                                                             >  1:274766/1‑148 (MQ=255)
                                                   cTTACCCGGAATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACAAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGAc                                                                                              <  2:274766/146‑1 (MQ=255)
                                                    ttACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCg                                                                                            >  6:24686/1‑147 (MQ=255)
                                                      aCCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGa                                                                                           <  6:186734/146‑1 (MQ=255)
                                                          ggCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCtt                                                                                     <  1:305040/148‑1 (MQ=255)
                                                          ggCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGGAATGACCATAAAGTTGATGATCACGAGAATGCCAACCACCACGAGACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCtt                                                                                     >  4:300987/1‑148 (MQ=255)
                                                               aCCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAg                                                                                                  <  5:33543/130‑1 (MQ=255)
                                                               aCCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAg                                                                                                  >  6:33543/1‑130 (MQ=255)
                                                               aCCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGac                                                                                 >  3:183604/1‑147 (MQ=255)
                                                                  aTCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAACCACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGaccacc                                                                             >  1:180373/1‑148 (MQ=255)
                                                                     gagaACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCtt                                                                           <  1:49221/147‑1 (MQ=255)
                                                                      agaACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTcc                                                                         <  8:276264/148‑1 (MQ=255)
                                                                          cAAAGCGCGCACCCACTTCTTCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgc                                                                     <  2:389438/148‑1 (MQ=255)
                                                                             aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgc                                                                  >  1:170226/1‑148 (MQ=255)
                                                                             aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgc                                                                  >  4:12316/1‑148 (MQ=255)
                                                                                gcgcACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccgc                                                               <  6:93619/148‑1 (MQ=255)
                                                                                       cACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcgccg                                                        <  8:145981/148‑1 (MQ=255)
                                                                                       cACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcgccg                                                        >  7:123233/1‑148 (MQ=255)
                                                                                         cTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGt                                                      <  6:170727/148‑1 (MQ=255)
                                                                                         cTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGt                                                      <  6:243379/148‑1 (MQ=255)
                                                                                           tCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTAt                                                    >  2:168166/1‑148 (MQ=255)
                                                                                                gATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCt                                               <  8:315213/148‑1 (MQ=255)
                                                                                                  tACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCAAGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGACCTTc                                             >  2:85141/1‑148 (MQ=255)
                                                                                                    cGCCCGTCACCTTTGTTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCCCCAACGAGGAAGTGACCGACCGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCTGTATGCCCTTcc                                            <  6:10282/147‑1 (MQ=255)
                                                                                                     gCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCAt                                          <  5:24686/148‑1 (MQ=255)
                                                                                                        cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGt                                                      >  2:22632/1‑133 (MQ=255)
                                                                                                        cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGt                                                      <  1:22632/133‑1 (MQ=255)
                                                                                                        cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTa                                        >  7:306968/1‑147 (MQ=255)
                                                                                                                ttGGTAATGACCATAAAGTTGATGATCTCGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTcc                                                                         <  3:159412/106‑1 (MQ=255)
                                                                                                                ttGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTcc                                                                         >  4:159412/1‑106 (MQ=255)
                                                                                                                ttGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                               >  2:343020/1‑148 (MQ=255)
                                                                                                                      aTGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGtt                          <  3:128573/147‑1 (MQ=255)
                                                                                                                       tGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTg                         >  2:26688/1‑147 (MQ=255)
                                                                                                                       tGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGa                        <  2:359412/148‑1 (MQ=255)
                                                                                                                        gACCATAAAGTTGATGATCACGTGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTg                         >  5:115383/1‑146 (MQ=255)
                                                                                                                        gACCATAAAGTTGATGATCACGTGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTg                         >  5:110482/1‑146 (MQ=255)
                                                                                                                        gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGAt                                 <  7:182302/138‑1 (MQ=255)
                                                                                                                        gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGAt                                 >  8:182302/1‑138 (MQ=255)
                                                                                                                        gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGa                        <  5:50320/147‑1 (MQ=255)
                                                                                                                         aCCATAAAGTTGATGATCACGTGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGaa                       <  6:110482/147‑1 (MQ=255)
                                                                                                                         aCCATAAAGTTGATGATCACGTGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGaa                       <  6:115383/147‑1 (MQ=255)
                                                                                                                         aCCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTg                         >  6:50320/1‑145 (MQ=255)
                                                                                                                           cATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGc                     >  8:303894/1‑147 (MQ=255)
                                                                                                                           cATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGc                     >  8:88517/1‑147 (MQ=255)
                                                                                                                            aTAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGcc                    <  8:97896/147‑1 (MQ=255)
                                                                                                                              aaaGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACg                 >  5:251872/1‑148 (MQ=255)
                                                                                                                              aaaGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACg                 >  3:185866/1‑148 (MQ=255)
                                                                                                                                 gTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTa              <  4:183604/148‑1 (MQ=255)
                                                                                                                                 gTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTa              <  4:185866/148‑1 (MQ=255)
                                                                                                                                         tCACGAGAAAGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGc         >  7:79260/1‑145 (MQ=255)
                                                                                                                                          cACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAAAGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGa     >  5:247811/1‑148 (MQ=255)
                                                                                                                                              agaATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGc  >  8:376810/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CTCATCTTCACCAATCAATCCGGCGTTAAGGTCGGCGTCAATCGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGC  >  NC_000913/1964511‑1964799

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