Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,856,021 T→G 25.5% pseudogene (943/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,856,0210TG25.5% 97.7 / 8.3 60pseudogene (943/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (16/28);  new base G (14/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 9.10e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG  >  CP009273/3855874‑3856169
                                                                                                                                                    |                                                                                                                                                      
cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGtt                                                                                                                                                       <  1:820034/149‑1 (MQ=255)
cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGt                                                                                                                                                        <  1:247024/148‑1 (MQ=255)
 gCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGc                                                                                                                                                     >  1:1266778/1‑150 (MQ=255)
  cccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCg                                                                                                                                                    <  1:334024/150‑1 (MQ=255)
  cccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCg                                                                                                                                                    >  1:954254/1‑150 (MQ=255)
           cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATctg                                                                                                                                           >  1:616531/1‑147 (MQ=255)
                cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGGAACAGCAAAcc                                                                                                                                     >  2:176766/1‑150 (MQ=255)
                 acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                    <  1:889687/150‑1 (MQ=255)
                 acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                    <  2:1442200/150‑1 (MQ=255)
                 acaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc                                                                                                                                    <  1:442983/150‑1 (MQ=255)
                  caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCg                                                                                                                                   <  1:65136/150‑1 (MQ=255)
                       aaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGATCAACAAAccccccccc                                                                                                                              >  1:113720/1‑147 (MQ=255)
                       aaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTTCGGGGAACAGGAAACCCGCCCCt                                                                                                                              >  2:640370/1‑150 (MQ=255)
                         gagaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTGGGGGGAACAGCAAGccccccccctt                                                                                                                            >  2:770874/1‑145 (MQ=255)
                            aTACATTAGGGTCGACATTGAGGCCGCCAGCAGGGCGTGTCCCGAAAACAGAAAGGGTGAAATGAAAAGGAGGGTAATTCCCAGCGTTTCGTTATTCCCCACCAGCATGGCGGTTAGGGTTGCGGGAATCAAAAAGCCCGCCACTTTTAc                                                                                                                         <  2:1361939/150‑1 (MQ=255)
                              aCATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGTGGGGAACAGCAAACCCCCCCCCTTTTccc                                                                                                                       >  1:114854/1‑150 (MQ=255)
                               cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGGATCAGCAAGCCCGCCCCCTTTTCCCg                                                                                                                      >  2:495654/1‑150 (MQ=255)
                               cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGGAACAGCAAGCCCGCCCCTTTTTccc                                                                                                                       >  2:1154012/1‑149 (MQ=255)
                                    cgctCGACATTGAGGCCGCCAGCACGGCGTGTCCCGCAAACAGCAACGGTAAAATGAACAGGAAGGTAAATCCCGGCGGTTGGGTAATCCCCACCAGCATGGCGGTTAGGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGAtttt                                                                                                                 <  1:173659/147‑1 (MQ=255)
                                                  gccgccAGCACGGGGTGTCCCGAAAACAGCAACGGTGAAATGAACAGGAAGGTAAATCCCAGGGGTTCGGTAATCCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTg                                                                                                  <  1:1154012/150‑1 (MQ=255)
                                                    cgccAGCGCGGCGTGTACCGCAAACAGAAACGGTAAAATGAACAGGAAGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAAGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGaa                                                                                                <  1:1350447/150‑1 (MQ=255)
                                                     gccAGCGCGGCGTGCACCGAAAACAGCAAGGGTGAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGGTTAAGGTTGGGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAg                                                                                               <  2:304093/150‑1 (MQ=255)
                                                     gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTGCGGGGATCAGCAAACCCGCCACTTTTACCCGATTTTCCGGGGCGGGAGTGGAg                                                                                               >  1:767106/1‑150 (MQ=255)
                                                     gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACAGCAAGCCCCCCCCTTTTACCCGATTTTTCGGGGCGGCAggg                                                                                                  >  2:180532/1‑145 (MQ=255)
                                                      ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCACTTTTACCCCATTTTTCGGGGGGGCAGTGAAg                                                                                               >  1:778869/1‑149 (MQ=255)
                                                       cggCAGGGCGTGTCCCGAAAAAAGAAAGGGTAAAAGGAAAGGGAAGATAAATCCCGGGGGTTGGTTAATCCCCCCCAGCAGGGGGTTAAGGGTGGGGGAAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                             <  2:711216/148‑1 (MQ=255)
                                                           cgggCGTGTCCCGCAAAAAGAAACGGAAAAATGAAAAGGAAGGTAAATTCCAGCGGTTCGGTAATCCCCCCCAGCAGGGCGTTTAGGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATc                                                                                         <  2:1176171/148‑1 (MQ=255)
                                                            gggggtgtCCCGAAAACAGCAAGGGTAAAATGAACAGGAGGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGGGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCg                                                                                        <  1:770874/146‑1 (MQ=255)
                                                                 tgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCCCCACTTTTTCCCGATTTTTCGGTGCGGCAGTGGAGTACATCGCtta                                                                                    >  2:956408/1‑147 (MQ=255)
                                                                 tgtACCGAAAACAGCAACGGTAAAATAAACAGGAAGGAAAATCCCAGCGTTTCGGTAATCCCCACCAGCATGGCGTTTAAGGTTGGGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaaa                                                                                   <  1:67692/150‑1 (MQ=255)
                                                                  gtcccGAAAAAAGAAACGGTAAAATGAACAGGAAGAAAAATCCCGGGGGTTCGATATTCCCCCCCAGCGGGGCGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAg                                                                                  <  1:180532/147‑1 (MQ=255)
                                                                    aCCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAACGTTGTGGGGATCAAGAAGCCCCCCACTTTTTTCCCATTTTTCCGTGCGGCAGTGAAGTACATCGCTAAAGa                                                                                 >  1:725643/1‑149 (MQ=255)
                                                                     ccGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGGGGTTCGGTAATCCCCCCCAGCATGGGGGTTAAGGTTGCGGGAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                               <  1:637852/150‑1 (MQ=255)
                                                                     ccGAAAACAGCAGCGGTAAAATGAACAGGAAGGTAAATCCCGGCGGTTCGGTAATCCCCACCAGCATGGCGTTTAGGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa                                                                               <  1:495654/150‑1 (MQ=255)
                                                                      cccAAACAGCAGGGGTGAAATAAACGGGAAGGAAAATCCGGGGGGTTCGGTAATCCCCCCAACCAGGGCGTTTAAGGTGGCGGGAATCAGAACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt                                                                              <  1:1094321/148‑1 (MQ=255)
                                                                       gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCCCCACCTTTTTCCGATTTTTCCGTGCGGCGAGGGAAGAACACGCTAAAAAAATg                                                                             >  2:96406/1‑150 (MQ=255)
                                                                        cAAACAGCAGCGGTAAAAAGAAAAGGAAGGTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCAGGGGGGTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGc                                                                            <  2:954934/150‑1 (MQ=255)
                                                                                aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGTGGGAGTGAAGGACATCGCTAAAAAAATGCCCACGGCa                                                                    >  2:1060056/1‑150 (MQ=255)
                                                                                    gaaaaaTAAAAAGGAAGGTAAATCCCGGGGTTTGGGTAATCCCCCCCAGCATGGGGGTTAAGGTTGCGGGAATCAAAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaa                                                                <  1:525490/147‑1 (MQ=255)
                                                                                       aaaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACCTTTTCCCGATTTTCCGGGGCGGCAGGGAAGTACATCGCTAAAAAAAAGCCCACGGCACCa                                                                 >  2:124457/1‑146 (MQ=255)
                                                                                         aTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACCGCACCAAAGAtt                                                           >  1:1114081/1‑150 (MQ=255)
                                                                                          tGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCCCCACTTTTTCCCGATTTTTCGGTGCGGCAGGGAAGTACATCGCTAAAAAAAAGCCCACGGCACCAAAGAttt                                                          >  1:553422/1‑150 (MQ=255)
                                                                                                ggAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGGAAGCCCGCCACCTTTACCCGATTTTCCGGGGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCACCAAAGATTTTTGAgg                                                    >  1:675537/1‑149 (MQ=255)
                                                                                                ggAAGGTAAATCCCAGCGGTTCGTTAATCCCCACCAGCATGGCGGTTAAGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGt                                                    <  1:239125/150‑1 (MQ=255)
                                                                                                       aaaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCACCTTTTCCCGATTTTTCGGTGCGGCAGTGGAGGACATCGCCAAAAAAAATCCCACCGCACCCAAGATTTTTTAGTTACCGTg                                             >  2:167605/1‑150 (MQ=255)
                                                                                                        attgccggCGGTTCGGTAATCCCCCCCAGCAGGGGGGTTAGGGTGGCGGGAATCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc                                            <  2:1266778/143‑1 (MQ=255)
                                                                                                        aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAACAAGCCCGCCACTTTTACCCCATTTTTCGGTGCGGCAGGGAAGTACATCGCTAAAAAAATTCCCACGGCACCCAAGATTTTTTAGTTAcccttc                                            >  2:1399478/1‑146 (MQ=255)
                                                                                                              ggCGGTTCGTTAATCCCCCCCAGCAGGGCGGTTAGGGTGGCGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCa                                      <  2:763899/149‑1 (MQ=255)
                                                                                                                 ggTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGGCCGCCCCTTTTACCCGATTTTTCGGGGCGGGAGGGGAGTACATCGCCAAAGAAAAGCCCACCGCACCAAAGATTTTTTAGGTAACGGGCAGGGCaaa                                    >  2:759409/1‑149 (MQ=255)
                                                                                                                        aaaTCCCCCCCAGAAGGGGGGTTAAGGTGGGGGAAACCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCtt                            <  2:1114081/149‑1 (MQ=255)
                                                                                                                               caccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATTCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaa                     >  1:1211548/1‑150 (MQ=255)
                                                                                                                                 ccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaaca                   <  1:956408/150‑1 (MQ=255)
                                                                                                                                      aTGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAt              <  1:1205771/150‑1 (MQ=255)
                                                                                                                                       tGGCGTTTAAGGTGGCGGAAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAtt             <  2:46842/150‑1 (MQ=255)
                                                                                                                                        ggCGGTAAGGGTGGCGGAAATCAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAt              <  2:616531/148‑1 (MQ=255)
                                                                                                                                            ggTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAg         >  2:752085/1‑150 (MQ=255)
                                                                                                                                            ggTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAg         >  2:1278464/1‑150 (MQ=255)
                                                                                                                                             tttAAGGTTGGGGAAACCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGc        <  1:1060056/149‑1 (MQ=255)
                                                                                                                                                 gggTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCt    <  2:553422/149‑1 (MQ=255)
                                                                                                                                                  ggTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc   >  2:986044/1‑150 (MQ=255)
                                                                                                                                                   gTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg  >  2:768683/1‑150 (MQ=255)
                                                                                                                                                   gTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg  >  2:1082686/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG  >  CP009273/3855874‑3856169

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: