Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,390,235 T→C I44T (ATC→ACC)  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,2350TC98.4% 198.4 / NA 63I44T (ATC→ACC) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/33);  minor base A (0/1);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATA  >  CP009273/4390091‑4390378
                                                                                                                                                |                                                                                                                                               
aagcagtggggagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTg                                                                                                                                             >  1:553895/11‑149 (MQ=255)
aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTg                                                                                                                                             >  2:195181/1‑149 (MQ=255)
 aGGAAAAGAGAGGAGGGCTAGGGGGCAATGTTTACAAGGTCCGTTCCTGGACGCCCTGCGTCGGGAGCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGt                                                                                                                                            <  1:361080/149‑1 (MQ=255)
      aaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGaaa                                                                                                                                       >  2:805187/1‑149 (MQ=255)
            gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACg                                                                                                                                 <  1:587335/149‑1 (MQ=255)
                 gCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAg                                                                                                                            >  2:56799/1‑149 (MQ=255)
                     aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAg                                                                                                                        >  2:365445/1‑149 (MQ=255)
                       gggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACAGTCAGCCAGAt                                                                                                                      >  2:633245/1‑149 (MQ=255)
                        ggCAATCTTTACAAGTGCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATg                                                                                                                     <  1:135247/149‑1 (MQ=255)
                            aTCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGttt                                                                                                                 <  2:179981/149‑1 (MQ=255)
                              cTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATTGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                                                                                                               <  1:264556/149‑1 (MQ=255)
                               tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTTTTGAAAAACACGGTCAGCCAGATGGTTTACa                                                                                                              >  2:597793/1‑149 (MQ=255)
                                ttACAAGATCCGGTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACaa                                                                                                             <  1:429282/149‑1 (MQ=255)
                                    aaGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCag                                                                                                         >  1:165842/1‑148 (MQ=255)
                                      gATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgc                                                                                                       >  1:730668/1‑149 (MQ=255)
                                        tCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGGTTACAAGCACGCGa                                                                                                     >  2:138789/1‑149 (MQ=255)
                                         gcgcgccgtgccgtaCTGCGGCGGGGCCGTGTTCCTGTTGCTACTGATGTGGTGAATGGCAATAAGCCGGTTGGGCAAATCGAGTCTTTTCTTCCGTTCGAAAAAAAAATAAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAt                                                                                                    <  1:545892/135‑1 (MQ=255)
                                           gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt                                                                                                  >  2:228182/1‑149 (MQ=255)
                                            ttCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc                                                                                                 >  1:622006/1‑149 (MQ=255)
                                              ccTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGGCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                                                                                               >  1:37436/1‑149 (MQ=255)
                                                 gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                            >  1:535584/1‑149 (MQ=255)
                                                      cACCGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTc                                                                                       <  2:1487/149‑1 (MQ=255)
                                                         tGCGTCGGGAACGTGTTCCAGTTTCTATTAATTTGCATGATAGTAGTAAGCTGCAAGGGCAAATCGAGTATTTTGATCAGTTCGTGACTCTGTTGAAAAACACGGGTAGCCAGATGGTTTACAAGCACGCGATTTCTACtgttgttcgg                                                                                    >  2:215189/1‑145 (MQ=255)
                                                           cGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGtc                                                                                  >  2:768365/1‑149 (MQ=255)
                                                                 gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAACACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGccc                                                                            <  1:725981/149‑1 (MQ=255)
                                                                      tgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                       <  2:200489/149‑1 (MQ=255)
                                                                        ttCCAGTTTCTAGTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                     <  1:88839/149‑1 (MQ=255)
                                                                           cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                                                  <  1:768624/149‑1 (MQ=255)
                                                                           cAGTTTCTAGTTATTTGGTGAATGGTATTAAGCGGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                                                  <  1:182218/149‑1 (MQ=255)
                                                                             gTTTCTATTTATTTGGGGAATTGTATTAAGCTGCACGGGCAAATCGAGCCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGCCAGCCAGATGGTTTACAAGCACGCGTTTGCTACTGTTGTCCCGTCTCGCCCGGTGTCtcatct                                                                <  1:864439/149‑2 (MQ=255)
                                                                               ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                              <  2:826577/149‑1 (MQ=255)
                                                                               ttCTATTTATTTGGTGAATGGTATTAAGCGGCAAGGGCAAATCGGGTCTTTTGATCAGTTCATGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                              <  1:12826/149‑1 (MQ=255)
                                                                                  tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTTTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTa                                                           <  1:195181/149‑1 (MQ=255)
                                                                                  tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTa                                                           <  1:218453/149‑1 (MQ=255)
                                                                                    tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGGCTTGTGGTCAGTTCGTGACCCTGTTGAAAAACACGCTCGGCCAGATGGTTTGCAAGGACGCGATTGCTACGGTTCTCCCGTGTCGCCCGGTTTCTCATGCCAGTgac                                                         >  1:813737/1‑149 (MQ=255)
                                                                                          tGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGcc                                                   <  1:138789/149‑1 (MQ=255)
                                                                                           ggTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCg                                                  >  2:412989/1‑149 (MQ=255)
                                                                                             tGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                                                >  1:90241/1‑149 (MQ=255)
                                                                                              gAATGGTATTAAGCTGCAAGGGCGAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCCTCACAGTAACAACGCCGGTg                                               >  2:584976/1‑149 (MQ=255)
                                                                                               aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                                              <  2:97620/149‑1 (MQ=255)
                                                                                                aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTTACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                             <  2:830190/149‑1 (MQ=255)
                                                                                                aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                             <  1:779584/149‑1 (MQ=255)
                                                                                                aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                             <  1:768365/149‑1 (MQ=255)
                                                                                                aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTCCTCATCACAGGAGCAACGCCGGTGGc                                             >  1:848/1‑149 (MQ=255)
                                                                                                      gtaggCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTAcc                                       <  1:205126/145‑1 (MQ=255)
                                                                                                           cTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagcag                                  <  1:365445/149‑1 (MQ=255)
                                                                                                               aaGGGCAAATCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                              <  2:604717/149‑1 (MQ=255)
                                                                                                                     aaaTCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCCGTAACTACcat                        >  2:747692/1‑149 (MQ=255)
                                                                                                                       atagAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatca                      <  2:444086/146‑1 (MQ=255)
                                                                                                                        tCGAGTCTTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCAAGATGGTTTTCAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTTCCAGCAGTAACTACcatcat                     >  2:749724/1‑149 (MQ=255)
                                                                                                                           aGTCCTTTGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGt                  <  1:56799/149‑1 (MQ=255)
                                                                                                                               ttttGATCAGTTCGTGCCCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCTATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagca              <  2:706118/149‑1 (MQ=255)
                                                                                                                               ttttGAGCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCCCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagca              <  1:5165/149‑1 (MQ=255)
                                                                                                                                 ttGATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagcagc            <  1:721113/149‑1 (MQ=255)
                                                                                                                                   gATCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgc          >  2:519963/1‑149 (MQ=255)
                                                                                                                                    aTCAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCACGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcg         >  2:485255/1‑149 (MQ=255)
                                                                                                                                     tCAGGTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAgcgcgc        >  1:458242/1‑149 (MQ=255)
                                                                                                                                      cAGTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCa       <  2:37436/149‑1 (MQ=255)
                                                                                                                                        gTTCGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGACTTCTACTGTTGTCCCGTCCCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGa     >  2:870581/1‑149 (MQ=255)
                                                                                                                                         ttCGTGACCCTATTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTATCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCCGCGCGCAGac    >  1:106746/1‑148 (MQ=255)
                                                                                                                                         ttCGGGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGaa    <  1:425321/149‑1 (MQ=255)
                                                                                                                                           cGTGACCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATa  >  2:641273/1‑149 (MQ=255)
                                                                                                                                           cGTGACCCGGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATa  <  1:357001/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATA  >  CP009273/4390091‑4390378

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