Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,637,683 G→A 6.4% H324Y (CAT→TAT)  bamB ← BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,6830GA 6.4% 200.2 / 4.0 62H324Y (CAT→TAT) bamBBamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/2);  total (30/32)
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.06e-01

TTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCCACATAATCTGA  >  NC_000913/2637543‑2637828
                                                                                                                                            |                                                                                                                                                 
ttCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATACAGCAGAt                                                                                                                                           >  1:114389/1‑149 (MQ=255)
ttCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATACAGCAGAt                                                                                                                                           >  1:335605/1‑149 (MQ=255)
     tCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCtt                                                                                                                                      >  2:32946/1‑149 (MQ=255)
      cTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCttt                                                                                                                                     <  1:14170/149‑1 (MQ=255)
      cTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCttt                                                                                                                                     <  3:140754/149‑1 (MQ=255)
      cTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCttt                                                                                                                                     <  4:341449/149‑1 (MQ=255)
         gAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTtgtg                                                                                                                                  <  3:98880/149‑1 (MQ=255)
           aaCCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTc                                                                                                                                >  3:312063/1‑149 (MQ=255)
            aCCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTcc                                                                                                                               >  1:314684/1‑149 (MQ=255)
                gAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAg                                                                                                                           >  4:245970/1‑149 (MQ=255)
                      tCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCttt                                                                                                                                     >  3:17630/1‑133 (MQ=255)
                      tCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCttt                                                                                                                                     <  4:17630/133‑1 (MQ=255)
                           tttttGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATgc                                                                                                                                                  <  2:232228/115‑1 (MQ=255)
                           tttttGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATgc                                                                                                                                                  >  1:232228/1‑115 (MQ=255)
                             tttGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCa                                                                                                              <  2:470885/149‑1 (MQ=255)
                                  tGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAAt                                                                                                         <  1:165757/149‑1 (MQ=255)
                                    ggCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATgg                                                                                                       >  2:350154/1‑149 (MQ=255)
                                       aacgaaacgaCCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCa                                                                                                    >  2:302815/1‑149 (MQ=255)
                                       aacgaaacgaCCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCa                                                                                                    >  2:301/1‑149 (MQ=255)
                                              cgaCCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCAt                                                                                             <  2:95337/149‑1 (MQ=255)
                                                  cATCTTCGACGTTAATCCAGTGCAGATAACCTTTACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCAcc                                                                                         >  3:301998/1‑149 (MQ=255)
                                                  cATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCAcc                                                                                         >  3:360750/1‑149 (MQ=255)
                                                   aTCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATACAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCAccc                                                                                        <  2:335605/149‑1 (MQ=255)
                                                   aTCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATACAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCAccc                                                                                        <  2:114389/149‑1 (MQ=255)
                                                     cTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCgg                                                                                      <  1:302815/149‑1 (MQ=255)
                                                           cGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAttt                                                                                <  2:383672/149‑1 (MQ=255)
                                                            gTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAtttt                                                                               <  1:263042/149‑1 (MQ=255)
                                                            gTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAtttt                                                                               >  1:242862/1‑149 (MQ=255)
                                                              tAATCCAGTGCAGCTAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGa                                                                             <  4:301998/149‑1 (MQ=255)
                                                               aaTCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGAt                                                                            <  4:418917/149‑1 (MQ=255)
                                                                 tCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCg                                                                          >  3:60360/1‑149 (MQ=255)
                                                                   cAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATc                                                                           <  1:343992/146‑1 (MQ=255)
                                                                   cAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATc                                                                           >  2:343992/1‑146 (MQ=255)
                                                                       gCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCaga                                                                    <  1:301/149‑1 (MQ=255)
                                                                           aTAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCagataga                                                                >  1:284524/1‑149 (MQ=255)
                                                                               ccTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCg                                                            >  2:300657/1‑149 (MQ=255)
                                                                                cTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCCGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGa                                                           <  3:208601/149‑1 (MQ=255)
                                                                                  tCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGAtt                                                         <  1:300657/149‑1 (MQ=255)
                                                                                      tGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCg                                                     >  4:362264/1‑149 (MQ=255)
                                                                                      tGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCg                                                     >  4:174094/1‑149 (MQ=255)
                                                                                        tCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTc                                                   <  1:350154/149‑1 (MQ=255)
                                                                                          aCCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTcga                                                 <  2:242862/149‑1 (MQ=255)
                                                                                          aCCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTcga                                                 <  3:18911/149‑1 (MQ=255)
                                                                                                 accaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAtt                                                                                 <  2:319987/110‑1 (MQ=255)
                                                                                                 accaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAtt                                                                                 >  1:319987/1‑110 (MQ=255)
                                                                                                 accaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGa                                                           >  2:229936/1‑132 (MQ=255)
                                                                                                 accaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGa                                                           <  1:229936/132‑1 (MQ=255)
                                                                                                 accaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGaa                                          <  3:174094/149‑1 (MQ=255)
                                                                                                 accaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGaa                                          >  4:167601/1‑149 (MQ=255)
                                                                                                  ccaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGaaa                                         >  3:121096/1‑149 (MQ=255)
                                                                                                  ccaGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGaaa                                         <  4:360750/149‑1 (MQ=255)
                                                                                                       ttGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCAtt                                    <  3:167601/149‑1 (MQ=255)
                                                                                                         gCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCa                                  >  3:470259/1‑149 (MQ=255)
                                                                                                                taCAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAAcc                           <  4:60360/149‑1 (MQ=255)
                                                                                                                   aGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAtt                                                                                 >  4:103112/1‑92 (MQ=255)
                                                                                                                   aGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCAtt                                                                                 <  3:103112/92‑1 (MQ=255)
                                                                                                                               gTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCca            <  2:314684/149‑1 (MQ=255)
                                                                                                                               gTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAACGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCca            >  4:342175/1‑149 (MQ=255)
                                                                                                                                tCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGGAATCATTCACCGAACCCAGTTCGCGTTTCcac           >  2:469832/1‑149 (MQ=255)
                                                                                                                                   gcagGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATTACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCCACATa        >  2:6338/1‑149 (MQ=255)
                                                                                                                                   gcagGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCCACATa        >  2:158437/1‑149 (MQ=255)
                                                                                                                                     agGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGAc                                                                        >  2:377782/1‑83 (MQ=255)
                                                                                                                                     agGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGAc                                                                        <  1:377782/83‑1 (MQ=255)
                                                                                                                                         gATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCCACATAATCTGa  <  2:284524/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGCAACGAAACGACCATCTTCGACGTTAATCCAGTGCAGATAACCTTCACTGTCACCGACCACCAGGTTGCCATTATACAGCACCGGAGAAGTCAGCAGGCGATGCAGCAGATCGCTTTGTGTCCACAGCGTAACGCCGCCATCAATGGTCAACGCCATCACCCGGTCATTTTGATCGACCAGATAGATGCGATTGCCGTCGACGATGAAATCATTCACCGAACCCAGTTCGCGTTTCCACATAATCTGA  >  NC_000913/2637543‑2637828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑