Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,286 C→A 77.0% T61N (ACT→AAT)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,390,2860CA77.0% 98.9 / 30.7 61T61N (ACT→AAT) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (12/2);  new base A (18/29);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

TGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTT  >  CP009273/4390148‑4390431
                                                                                                                                          |                                                                                                                                                 
tgtgTCGGGATCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCg                                                                                                                                         <  1:331062/146‑1 (MQ=255)
tGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGAAAGGGCAAATCGAGTCTTTTGAGCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCg                                                                                                                                         <  1:339285/149‑1 (MQ=255)
tGCGTCGGGAACGTGTTCCAGTTTCAATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCg                                                                                                                                         >  1:692834/1‑149 (MQ=255)
   gTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGAGCGACCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGtct                                                                                                                                      <  1:1213606/149‑1 (MQ=255)
     cGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCg                                                                                                                                    <  2:892956/149‑1 (MQ=255)
      gggAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGc                                                                                                                                   >  1:76137/1‑149 (MQ=255)
           cGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt                                                                                                                              >  1:708497/1‑149 (MQ=255)
            gtgtTCGCGTTTCTATTTATGTGGTGAATGGTATTAAGATGCAAGTGCAAATCGAGTCTTTTGAGCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAAGGTTGACCCGTCTCGCCCGGtt                                                                                                                             <  2:217141/149‑1 (MQ=255)
              gtTGCAGTTTCTATTTATTTGGTGAATGGTATGAAGCTGCAAGGGGAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtc                                                                                                                           <  2:476745/149‑1 (MQ=255)
               ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtct                                                                                                                          >  2:1092410/1‑149 (MQ=255)
                ttcAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtctc                                                                                                                         <  2:104685/147‑1 (MQ=255)
                tCCAGTTTCTATTTATTTGGTGAATGGTAGTAAGCTGCAAGCGCACATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGGTGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGACTCGCCCGGTTtctc                                                                                                                         <  2:244017/149‑1 (MQ=255)
                 ccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCtca                                                                                                                        >  1:862370/1‑149 (MQ=255)
                     tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATcac                                                                                                                    >  2:1233888/1‑149 (MQ=255)
                      ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                                                                   <  2:973440/149‑1 (MQ=255)
                         tatttattTGGTTAATGGTATTAAGCGGTAAGGGAAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCGGGTTTCTCATCACAGTa                                                                                                                <  1:647782/149‑1 (MQ=255)
                         tatttattTGGTGAATGGTATTAAGCTGCAAGGACAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCGGTCTCGCCCGGTTTATCATCACAGTa                                                                                                                >  2:97494/1‑149 (MQ=255)
                           tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaac                                                                                                              >  1:818599/1‑149 (MQ=255)
                                   gTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATAAGTTCGTGATCCTGTTGAAAAACACGGTCAGACAGATGGTTTACAAGAACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTTTCATCCAAGTAACAAAGACgg                                                                                                      >  2:807312/1‑149 (MQ=255)
                                       aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                                                                                  <  1:254105/149‑1 (MQ=255)
                                       aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                                                                                  <  2:278094/149‑1 (MQ=255)
                                        tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                                                                 >  1:1162779/1‑149 (MQ=255)
                                        tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                                                                 <  2:842376/149‑1 (MQ=255)
                                        tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCACATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTGTCCCATCACAGTAACACCGCCGGTGGCg                                                                                                 >  1:949272/1‑149 (MQ=255)
                                            agtAAGCTGCAAGGGAAATTCGGGTCTTTTGATCAGTTCGTGATCCTGTTGAAAACGACGGTCAGCCAGGTGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTAc                                                                                             <  2:191371/147‑1 (MQ=255)
                                                  cTGCAAGGGCAAAACGAGGCTTTTGATCAGTTCGAGAACCTGTTGAAAAACACGGTCAGCCAGATGGATTACCAGCACGCGATTTCTAAGGTTGTACAGCCGGGCCCGGTTTCTCATCACAGGAACACCGCCGGTGGCGGTGGcagcgg                                                                                       >  2:123567/1‑147 (MQ=255)
                                                   tGCAAGGGCAAATAGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGACAGCCAGAAGGTGTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGATgcagcggccagcagt                                                                                      >  2:998114/1‑135 (MQ=255)
                                                     cAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTaa                                                                                    <  1:1119801/149‑1 (MQ=255)
                                                      aaGGGCAAATCGAGTTTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                                                                                   <  1:1061473/149‑1 (MQ=255)
                                                       aGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAACAGAACGCGACCTCTAATAGTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACt                                                                                  <  1:600435/149‑1 (MQ=255)
                                                         ggCAAATCGAGTCTTTTGATCAGTTCATGACCCTGTTGAAAAACACGGTCAGCTAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAc                                                                                <  1:474416/149‑1 (MQ=255)
                                                                  agccTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGt                                                                       >  2:1076955/4‑149 (MQ=255)
                                                                   gTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTa                                                                      >  1:1141810/1‑149 (MQ=255)
                                                                    tCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTag                                                                     <  1:383747/149‑1 (MQ=255)
                                                                      ttttGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACCAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagca                                                                   <  1:211535/149‑1 (MQ=255)
                                                                               gTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCGTCATGGGCGCAGCGCGCAGa                                                          >  2:382333/1‑149 (MQ=255)
                                                                                ttCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGaa                                                         <  1:736018/149‑1 (MQ=255)
                                                                                       tCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTcc                                                  <  1:1109489/149‑1 (MQ=255)
                                                                                       tCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTcc                                                  <  1:206873/149‑1 (MQ=255)
                                                                                        ccTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcg                                                 >  2:648540/1‑149 (MQ=255)
                                                                                        ccTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGAAGAATACTTCcg                                                 >  2:30386/1‑149 (MQ=255)
                                                                                         cTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc                                                <  2:1196000/149‑1 (MQ=255)
                                                                                           gTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCgcgc                                              >  1:514211/1‑149 (MQ=255)
                                                                                           gTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGGTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCgcgc                                              <  1:1116305/149‑1 (MQ=255)
                                                                                            ttGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCGCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGGATACTTCCgcgcg                                             >  1:222529/1‑148 (MQ=255)
                                                                                                  aaCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGCAGCAGCGCGCAGAATACTTCCGCGCAACAGGa                                       <  2:692834/149‑1 (MQ=255)
                                                                                                   acacGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTTATCACAGTAACAACGCCGGTGGCGGTTCCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGAc                                      >  1:48571/1‑149 (MQ=255)
                                                                                                       ggTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCg                                  >  1:73338/1‑149 (MQ=255)
                                                                                                       ggTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCg                                  >  1:194051/1‑149 (MQ=255)
                                                                                                       ggTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTAGCGGTACCAGCAGTAACTACCATCATGGTAGCAACGCGCAGAATACTTCCGCGCAACAGGACAGCg                                  >  2:354584/1‑149 (MQ=255)
                                                                                                         tCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACAAGCAGTAACTACCATCATGGTAGCCGCGCGAAGGATACTTCCGCGCAACAGGACGGCgac                                >  1:444929/1‑148 (MQ=255)
                                                                                                           aGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACGGTAACAACGCCGGTGGCGGTACCAGCAATAACTACCCTCAAGGTAGCAGCGCGCAGCATACGTCCGCGCAACAGGACAACgaaga                              >  2:474850/1‑149 (MQ=255)
                                                                                                           aGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAAACCTTCCGCGCAACAGGACAGCgaaga                              >  2:44086/1‑149 (MQ=255)
                                                                                                           aGCCAGATGGTTTACAAGCACGCGATTTCGAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCgaaga                              <  1:471044/149‑1 (MQ=255)
                                                                                                              gtgATGGGTTACAAGCACGCCATTGTCAATGTGGTCCCGTCTCGCCCGGTTTCTCATCAGGCTATCCACGCCGTTGGCGTTACCAGCAGTAACTACCACGATCGCTGAAGCTCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAAc                           <  2:263243/147‑1 (MQ=255)
                                                                                                              cAGATGGGTGACAAGCACGCGATTTCTAATGGTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAAc                           <  2:931872/149‑1 (MQ=255)
                                                                                                                               cgcgATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGc          <  2:76137/149‑1 (MQ=255)
                                                                                                                                gcgATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCt         >  1:931362/1‑149 (MQ=255)
                                                                                                                                gcgATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGATACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAGCAGGGCAGCGAAGAGACCGAATAAGGTGCCGGGCt         >  1:404986/1‑149 (MQ=255)
                                                                                                                                     ttCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGtttt    >  1:777929/1‑149 (MQ=255)
                                                                                                                                       cTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGGGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGtttttt  <  2:249181/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
TGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTT  >  CP009273/4390148‑4390431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑