Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 289,794 A→G 7.0% V239V (GTT→GTC)  yagM ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273289,7940AG 7.0% 190.7 / 2.5 69V239V (GTT→GTC) yagMCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base G (3/2);  total (40/29)
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.88e-01

TACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAA  >  CP009273/289650‑289941
                                                                                                                                                |                                                                                                                                                   
tccccctTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATCGTGTCATCCGGTGCATTAACACGGTTGGAAACATCTGTTTTGCGGAGCGCTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCg                                                                                                                                                 <  2:438919/144‑1 (MQ=255)
 aCCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACGACCgg                                                                                                                                                >  2:1180532/1‑149 (MQ=255)
 aCCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACGACCgg                                                                                                                                                >  2:194092/1‑149 (MQ=255)
 aCCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAAAATCTGTTTTGCGGAGCACTGGCAGCGTAACCACCATTTCATTGACATGGTAAATATCGATCCACCAAACTACCgg                                                                                                                                                >  2:1180537/1‑149 (MQ=255)
      ttAGGTTCTGCTCTTTTGTGAGTATGGGTTTATAAGGATGTCATCCGTTACATTAACACGGTTGGCAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAg                                                                                                                                           <  2:3833/149‑1 (MQ=255)
      ttAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTAATTGACATGGTAAATATTAATCCAACAAACAACAGGATaac                                                                                                                                           >  1:447022/1‑148 (MQ=255)
       tAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAgg                                                                                                                                          >  2:594/1‑149 (MQ=255)
          gTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCCTCAGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGGAAATATCAATCCAACAAACGACCGGCGAAGGCAt                                                                                                                                       >  1:1060876/1‑149 (MQ=255)
          gTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGACACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCACTCCAACAAACAACCGGATAAgaaat                                                                                                                                       >  1:937241/1‑145 (MQ=255)
          gTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCAt                                                                                                                                       >  1:529679/1‑149 (MQ=255)
           ttCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATCTCAATCCAACAAACAACCGGATAAGGCAtc                                                                                                                                      >  2:175014/1‑148 (MQ=255)
               gCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACaa                                                                                                                                  >  1:378231/1‑149 (MQ=255)
                  ctctTGTGAGTATGGGTTTATAATGCTGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATcc                                                                                                                               <  2:220385/149‑1 (MQ=255)
                   tctTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGTAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATccc                                                                                                                              >  1:976180/1‑149 (MQ=255)
                           gTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTttcttc                                                                                                                      <  2:433723/149‑1 (MQ=255)
                             aTGTGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTa                                                                                                                    <  1:91878/149‑1 (MQ=255)
                                ggTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGa                                                                                                                 <  1:944082/149‑1 (MQ=255)
                                 gTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAGTCCAACAAACAGCCGGAGATGGCATAACAATCCCCTTTCTTCTAAGaa                                                                                                                <  2:358813/149‑1 (MQ=255)
                                    tataATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACtt                                                                                                             >  1:55141/1‑149 (MQ=255)
                                    tataATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACtt                                                                                                             >  1:934617/1‑149 (MQ=255)
                                    tataATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACtt                                                                                                             >  1:220520/1‑149 (MQ=255)
                                    tataATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTATTCTAAGAACtc                                                                                                             >  1:219914/1‑148 (MQ=255)
                                             gTCATCCGGTACATTAACACGGTTGGAAACACCTGGTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCATAGAACTCCATTTATcc                                                                                                    >  2:746280/1‑149 (MQ=255)
                                                 tCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGaa                                                                                                >  1:1232169/1‑149 (MQ=255)
                                                          ttAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAGATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTGTCCCGAATATAATGaa                                                                                       >  2:831726/1‑149 (MQ=255)
                                                          ttAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCaa                                                                                       >  2:1200299/1‑149 (MQ=255)
                                                          ttAACACGGTTGGAAACATCTGTATTGCGGAGCACTGGCAACATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCaa                                                                                       >  1:668819/1‑149 (MQ=255)
                                                            aaCACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAAtt                                                                                     >  1:1202570/1‑149 (MQ=255)
                                                              cacGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAACCACCAGGATAAGGCATAACACTCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTgg                                                                                   <  2:876535/149‑1 (MQ=255)
                                                              cacGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTgg                                                                                   <  2:73003/149‑1 (MQ=255)
                                                               acGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTggg                                                                                  >  2:375016/1‑149 (MQ=255)
                                                                cGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGa                                                                                 <  1:789779/149‑1 (MQ=255)
                                                                cGGTTGGAAACATCTGTTTTGCGGAGCACTCGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGa                                                                                 <  2:807162/149‑1 (MQ=255)
                                                                  gTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTggaaag                                                                               >  2:106037/1‑145 (MQ=255)
                                                                   ttGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGt                                                                              >  2:1027745/1‑149 (MQ=255)
                                                                   ttGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATAAATTGGGATGt                                                                              >  1:984459/1‑149 (MQ=255)
                                                                     ggAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTg                                                                            <  2:1133326/149‑1 (MQ=255)
                                                                      gAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACGATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAGGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGa                                                                           >  2:350271/1‑149 (MQ=255)
                                                                      gAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGa                                                                           <  1:213539/149‑1 (MQ=255)
                                                                              tGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTa                                                                   <  1:220592/149‑1 (MQ=255)
                                                                                ttttGCGTAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTaga                                                                 <  2:976180/149‑1 (MQ=255)
                                                                                  ttGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGc                                                               <  2:1014461/149‑1 (MQ=255)
                                                                                  ttGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGc                                                               <  2:4303/149‑1 (MQ=255)
                                                                                            cTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTCCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCa                                                     <  2:1231496/149‑1 (MQ=255)
                                                                                               gcagcaTAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCa                                                  >  1:1119941/1‑149 (MQ=255)
                                                                                               gcagcaTAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCa                                                  >  2:532537/1‑149 (MQ=255)
                                                                                               gcagcaTAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCa                                                  >  2:1147104/1‑149 (MQ=255)
                                                                                                cagcaTAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCaa                                                 <  1:1200299/149‑1 (MQ=255)
                                                                                                 agcaTTACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCaaa                                                <  2:447022/149‑1 (MQ=255)
                                                                                                 agcaTAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCaaa                                                >  1:1040417/1‑149 (MQ=255)
                                                                                                     tAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGtaa                                            >  1:1170797/1‑149 (MQ=255)
                                                                                                       accaccATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAAAGCCTGCCCGGCACCAAAGtaaaa                                          >  2:1206418/1‑147 (MQ=255)
                                                                                                            caTTTCATTGACATGGTAAATATCAATCCAACAAACGACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCt                                     <  1:194092/149‑1 (MQ=255)
                                                                                                            caTTTCATTGACATGGTAAATATCAATCCAACAAACGACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCt                                     <  1:1180532/149‑1 (MQ=255)
                                                                                                            caTTTAATTGACATGGTAACTATCAATCCAACAAACTACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCt                                     <  1:1180537/149‑1 (MQ=255)
                                                                                                             aTTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTc                                    <  2:1202570/149‑1 (MQ=255)
                                                                                                               ttCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGGTGACAAGAGTAGAGCCTGCCCGGCCGCAAAGTAATAACCCTCaa                                  >  2:197471/1‑149 (MQ=255)
                                                                                                                    tGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGgtt                             <  1:984980/149‑1 (MQ=255)
                                                                                                                     gACATGGAAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGgttg                            <  1:191216/149‑1 (MQ=255)
                                                                                                                        aTGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTc                         >  1:977252/1‑149 (MQ=255)
                                                                                                                        aTGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTACAAGGTTGTTc                         >  1:741524/1‑149 (MQ=255)
                                                                                                                        aTGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCATCAAAGGTTGTTc                         >  2:949445/1‑149 (MQ=255)
                                                                                                                        aTGGTAAATACCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTAAAAGGTTGTTc                         >  1:977043/1‑149 (MQ=255)
                                                                                                                          ggTAAATATGAATCCAACAAACAACTGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTa                       <  2:772505/149‑1 (MQ=255)
                                                                                                                           gTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGCATGTTGAAAAGAGTAGAGCCTGCCGGGCATCAAAGTAATAACCCTCCAAgggagtgctag                      >  2:1073961/1‑140 (MQ=255)
                                                                                                                                 ttCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTAGCCCGAATATAATCAATTTGGATGTTGACAAGAGTAGAGCCTGCCCTGCATCAAAGTAATATCCCTCAAAGGGTGTTCTATgaataa                <  2:937241/148‑1 (MQ=255)
                                                                                                                                 atccatccAGCAAACAACCGGATTAGGCATACCAATCCCCTTTCTTCTACGTCGTTCATTTATCCCGAATACGACCAATGGGGTTTTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATgaataa                <  2:1060876/149‑1 (MQ=255)
                                                                                                                                 atCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATgaataa                >  2:213800/1‑149 (MQ=255)
                                                                                                                                     aTCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATAGAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCACAAAGGTTGTATTATgaataagaaa            >  1:754795/1‑148 (MQ=255)
                                                                                                                                         aacaaacaACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATg        >  1:173825/1‑149 (MQ=255)
                                                                                                                                               caacCGGATAAGGCATAACAATCCCCTTTCTGCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGACCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGaa  <  1:831726/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TACCTCTTAGGTTCTGCTCTCTTGTGAGTATGGGTTTATAATGATGTCATCCGGTACATTAACACGGTTGGAAACATCTGTTTTGCGGAGCACTGGCAGCATAACCACCATTTCATTGACATGGTAAATATCAATCCAACAAACAACCGGATAAGGCATAACAATCCCCTTTCTTCTAAGAACTTCATTTATCCCGAATATAATCAATTGGGATGTTGACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAA  >  CP009273/289650‑289941

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