Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,390,286 C→A 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,2860CA98.5% 221.5 / NA 68T61N (ACT→AAT) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (33/34);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTT  >  CP009273/4390145‑4390431
                                                                                                                                             |                                                                                                                                                 
catggCGTCGGGAACGTGTTCAAGTTTCTATTTATTTGGTGGATGGTATTAAGCGGTAAGGCCAAATCGAGTCTTTTGATCAGTGCGTTACCCTGTTGAAAATCACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTc                                                                                                                                            <  2:638162/145‑1 (MQ=255)
  cTGCGTCGGGAACGTGTTCCAGTTTCTATTTGTGTGATGAAAGGTAATAAGCTGGAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATAAGTTACAAGCACACGAATCCTAAAGTTGTccc                                                                                                                                          >  2:1160079/1‑149 (MQ=255)
           gcatgtgtTCCAGTCTCTATTTACTGTGTGCATGGTTTTAAGCTGCAAGGGCAAATGAAGTCTTTTGATCAGTTCTTGATCCTGTTCAAAAACACGGTCATCCAGATGGTTTACAAGCACGCGTTTTCTAATGTTGTCCCGTCTCGccc                                                                                                                                 <  2:1004808/145‑1 (MQ=255)
           gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGccc                                                                                                                                 <  2:1169940/149‑1 (MQ=255)
           gAACGTGGTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGccc                                                                                                                                 <  1:499452/149‑1 (MQ=255)
            aaCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCg                                                                                                                                >  1:749367/1‑149 (MQ=255)
            aaCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCg                                                                                                                                <  1:779769/149‑1 (MQ=255)
               gtgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGtt                                                                                                                             <  2:976512/149‑1 (MQ=255)
                 gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtc                                                                                                                           <  2:910439/149‑1 (MQ=255)
                 gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtc                                                                                                                           <  2:1015324/149‑1 (MQ=255)
                  tgccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtct                                                                                                                          <  1:146999/147‑1 (MQ=255)
                   tCCAGTTTCTATTTAGTTGGTGATGGGTATGAGTCCGCCAGGGCAAGTCGAGTCTTTTGATCAGTGCGTGGTCCTGTGGAAAAACACGGTCAGTCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtctc                                                                                                                         <  2:474516/149‑1 (MQ=255)
                   cccAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTtctc                                                                                                                         <  1:1063400/148‑1 (MQ=255)
                     cAGTTTCGTTTTATTGGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGCTGAAAATCACTGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCtcat                                                                                                                       <  1:1172184/149‑1 (MQ=255)
                       gTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCtcatca                                                                                                                     <  2:1318902/149‑1 (MQ=255)
                             atttatgtGGTGAAGTGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTATTCACACTTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                                                                                               <  1:1311466/149‑1 (MQ=255)
                              tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaac                                                                                                              >  1:366397/1‑149 (MQ=255)
                              tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAACCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaac                                                                                                              <  2:411880/149‑1 (MQ=255)
                                tatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGAAACCGAGTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaa                                                                                                            >  1:961929/1‑149 (MQ=255)
                                       tGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCAAGATGGGTTACAAGCACGAGATAAATTAAATCTTCCCGACTCGGCCGGTGACTGATCACAGTAACAACGCCGGt                                                                                                     >  1:632853/1‑149 (MQ=255)
                                        gAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTg                                                                                                    <  1:478295/149‑1 (MQ=255)
                                         aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                                                                                                   >  2:1170115/1‑149 (MQ=255)
                                           tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTATAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                                                                 <  1:815473/149‑1 (MQ=255)
                                            ggTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCgg                                                                                                <  1:182530/149‑1 (MQ=255)
                                             gTATTAAGCTGCAATGGCAAATCGAGTCTTTTCTTAAGTTCGTGATCCTGTTTAAAACCACGGTCAGCCAGATCCTTTTCATTCACGCGCTTTCTAATGTTGTCCCGTCTCGACCGGTTTCTCATCACAGTAACACCGCCGGTGGCGGt                                                                                               <  2:806934/149‑1 (MQ=255)
                                             gTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAAAGATGTCCCGTCTCGCCCGGTTTCTCATTCCAGTAACAACGCCGGTGGCGGt                                                                                               >  2:871330/1‑149 (MQ=255)
                                              tATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTa                                                                                              <  1:969690/149‑1 (MQ=255)
                                              tATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGGTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTa                                                                                              <  1:610307/149‑1 (MQ=255)
                                                 tAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACca                                                                                           <  1:726190/149‑1 (MQ=255)
                                                  aaGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                                                          <  1:865720/149‑1 (MQ=255)
                                                   aGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGCTTTCTCATCCCCGTAACAACGCCGGTGGCGGTACccgc                                                                                         >  2:79500/1‑149 (MQ=255)
                                                     ctgtgaGGGCAAATCGTGTTTTTTGATCAGTTCGTGATCCTGTTGAAAGACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagcag                                                                                       <  1:394776/144‑1 (MQ=255)
                                                         aaGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACGGTACCAACGCCGGTGGCGGTACCAGCAGTAAc                                                                                   >  1:26588/1‑149 (MQ=255)
                                                         aaGGGCAAATCGAGTCTTTTGATCAGTTCGAGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAATAAc                                                                                   >  1:688220/1‑149 (MQ=255)
                                                             gCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTCCCAGCAGTCACTAcc                                                                               >  2:149178/1‑149 (MQ=255)
                                                               aaaTCGTGTCTTTTGATCAGTTCGTGACACTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcat                                                                             <  2:1006885/149‑1 (MQ=255)
                                                                    gAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATga                                                                        >  1:1085149/1‑148 (MQ=255)
                                                                       tCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTag                                                                     >  1:1152015/1‑149 (MQ=255)
                                                                        ctttgtaTCAGTTCGTCGTCCCGCTGAAAACCACGGTCAGCCAGTGGGTTTACAAGCACGCCATTCCTAATGTTGCCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagc                                                                    <  1:4246/143‑1 (MQ=255)
                                                                         ttttGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGGACCAGCAGTAACTACCATCATGGTagca                                                                   >  1:490892/1‑149 (MQ=255)
                                                                               tCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTCCCCTCATGGTAGCAgggcgg                                                             >  2:855972/1‑148 (MQ=255)
                                                                               tCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc                                                             >  1:49445/1‑149 (MQ=255)
                                                                                  gTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGa                                                          >  2:1257310/1‑149 (MQ=255)
                                                                                     cGTGATCCTGTTGAAAAACCCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATa                                                       <  1:233280/149‑1 (MQ=255)
                                                                                          tCCTGTTTAATAACACGGTCAGCTTGATGGTTTACAAGCATGCGATGGCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTcc                                                  <  1:695283/149‑1 (MQ=255)
                                                                                             tGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgcg                                               >  1:1195035/1‑149 (MQ=255)
                                                                                             tGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgcg                                               >  2:470086/1‑149 (MQ=255)
                                                                                               ttGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCa                                             >  2:824731/1‑149 (MQ=255)
                                                                                                tGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa                                            <  2:97827/149‑1 (MQ=255)
                                                                                                         cGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGc                                   <  1:215637/149‑1 (MQ=255)
                                                                                                          ggTCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCg                                  >  1:1315323/1‑149 (MQ=255)
                                                                                                            tCAGCCAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCgaa                                <  2:457376/149‑1 (MQ=255)
                                                                                                                ccAGATGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGaaa                            <  1:410012/149‑1 (MQ=255)
                                                                                                                   gATGGTTGACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCg                         >  2:349143/1‑149 (MQ=255)
                                                                                                                     tGGTTTACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGaa                       >  2:651783/1‑149 (MQ=255)
                                                                                                                        tttACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATaa                    <  1:1201483/149‑1 (MQ=255)
                                                                                                                         ttACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAg                   >  2:766131/1‑149 (MQ=255)
                                                                                                                         ttACAAGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTACGCGCAACAGGACAGCGAAGAAACCGAATAAg                   >  1:1027970/1‑149 (MQ=255)
                                                                                                                             aaGCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGttt               >  1:76905/1‑149 (MQ=255)
                                                                                                                               gCACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAAGAACGTCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGGTTCg             >  1:915347/1‑149 (MQ=255)
                                                                                                                                cACGCGATTTCTAATGTTGTCCCGTCTCGCCCGGATTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTCGCAGCGCGCAGAATACTTCCGCGCAACAGGCCAGCGAAGATACCGAATAAGGTTGCgg            >  2:744349/1‑149 (MQ=255)
                                                                                                                                  cgcgATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGc          <  2:352320/149‑1 (MQ=255)
                                                                                                                                   gcgATTTCTCATGTTGTTACGTCTCGCCCGGTTTCCCACAGCAGTGACAACCCCCGTAGCGGTACCAGCAGCAACCACCATCATGGTAGCAGAGCGCAGAAGACTGCACCGCACCAGGACAGCGAAGAAACCGAATAGGATTTCCGGca         >  2:557279/1‑148 (MQ=255)
                                                                                                                                     gATTTCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGt       >  1:681381/1‑149 (MQ=255)
                                                                                                                                     gATTTCTAATGTTCTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAACAGCAACTACCATCATGGTAGCAGCGCGCAGAATAGTACAGCACAACAGGACAGCGAAGAAACCGAAACAGGGTGCGGGCTGt       >  2:305408/1‑149 (MQ=255)
                                                                                                                                         tCTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTACCTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGCCAGCGAAGAAACCGAATCAGGGTTCGGGCTGttttt   >  1:1174198/1‑149 (MQ=255)
                                                                                                                                          cTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGtttttt  <  1:470086/149‑1 (MQ=255)
                                                                                                                                          cTAATGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGtttttt  >  1:881507/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTT  >  CP009273/4390145‑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: ‑