Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,633,006 C→A 19.6% intergenic (‑82/+281) cspI ← / ← rzpQ Qin prophage; cold shock protein/Rz‑like protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,633,0060CA19.6% 627.2 / 102.4 234intergenic (‑82/+281)cspI/rzpQQin prophage; cold shock protein/Rz‑like protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (108/80);  new base A (14/32);  total (122/112)
Fisher's exact test for biased strand distribution p-value = 1.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 236 total aligned reads displayed.

GATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATT  >  CP009273/1632862‑1633152
                                                                                                                                                |                                                                                                                                                  
gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCa                                                                                                                                                 <  2:1565266/148‑1 (MQ=255)
 atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCaa                                                                                                                                                >  1:2522318/1‑148 (MQ=255)
  tgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAag                                                                                                                                               <  2:958224/148‑1 (MQ=255)
   gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAaga                                                                                                                                              >  2:2003490/1‑148 (MQ=255)
    aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAagag                                                                                                                                             <  2:2454404/148‑1 (MQ=255)
    aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAagag                                                                                                                                             <  1:505036/148‑1 (MQ=255)
     aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAgg                                                                                                                                            >  2:2683566/1‑148 (MQ=255)
     aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAgg                                                                                                                                            <  1:1981978/148‑1 (MQ=255)
          aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACaa                                                                                                                                       <  2:3216521/148‑1 (MQ=255)
            aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAaat                                                                                                                                     <  2:1304106/148‑1 (MQ=255)
            aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAaat                                                                                                                                     <  2:2654892/148‑1 (MQ=255)
              cTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataa                                                                                                                                   >  1:1732340/1‑148 (MQ=255)
                     aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAtt                                                                                                                            <  1:2594331/148‑1 (MQ=255)
                     aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAAtt                                                                                                                            <  1:248813/148‑1 (MQ=255)
                     aGGGTTAAACCATGTCACTACCCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAtt                                                                                                                            <  2:2101541/148‑1 (MQ=255)
                        gTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAACTCTCAAGAGGTACAAATAATAAAATtata                                                                                                                         >  1:1618923/1‑148 (MQ=255)
                        gTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATtata                                                                                                                         <  1:2663354/148‑1 (MQ=255)
                            aaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACaa                                                                                                                     >  1:1821354/1‑148 (MQ=255)
                             aCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAt                                                                                                                    >  2:1529127/1‑148 (MQ=255)
                             aCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAAt                                                                                                                    >  2:3139630/1‑148 (MQ=255)
                              ccATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATg                                                                                                                   >  1:3102888/1‑148 (MQ=255)
                                aTTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAc                                                                                                                 <  1:335621/148‑1 (MQ=255)
                                aTTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAc                                                                                                                 >  1:1281541/1‑148 (MQ=255)
                                aTTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAc                                                                                                                 >  1:2396477/1‑148 (MQ=255)
                                aTTTCACTAAACCAGGCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGTTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAc                                                                                                                 >  1:3158734/1‑148 (MQ=255)
                                 tttCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACt                                                                                                                >  1:1794444/1‑148 (MQ=255)
                                    cACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCt                                                                                                             <  1:1225901/148‑1 (MQ=255)
                                      cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTc                                                                                                           >  2:2919652/1‑148 (MQ=255)
                                      cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTc                                                                                                           >  2:1610620/1‑148 (MQ=255)
                                      cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTc                                                                                                           <  1:3063540/148‑1 (MQ=255)
                                         aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGa                                                                                                        <  2:1266422/148‑1 (MQ=255)
                                         aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGa                                                                                                        <  1:814402/148‑1 (MQ=255)
                                            cAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATaaa                                                                                                                               >  1:1496003/1‑122 (MQ=255)
                                                 aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTg                                                                                                <  2:477064/148‑1 (MQ=255)
                                                 aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTg                                                                                                <  2:2088236/148‑1 (MQ=255)
                                                 aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTg                                                                                                >  1:138572/1‑148 (MQ=255)
                                                   tttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTa                                                                                              >  1:549771/1‑148 (MQ=255)
                                                    ttGTTATACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTaa                                                                                             <  2:2629304/148‑1 (MQ=255)
                                                    ttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAATAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTaa                                                                                             >  2:1840286/1‑148 (MQ=255)
                                                    ttGTTAGACATAATTATTAACTTTTGAAGAAATTAGCCCTTGGGCAGAATGTTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTaa                                                                                             <  1:3199332/148‑1 (MQ=255)
                                                        tAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACaaa                                                                                         >  1:1123459/1‑148 (MQ=255)
                                                        tAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACaaa                                                                                         >  1:871422/1‑148 (MQ=255)
                                                         aGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                                                                                        <  1:1566790/148‑1 (MQ=255)
                                                         aGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                                                                                        >  1:2419292/1‑148 (MQ=255)
                                                          gACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAcc                                                                                       <  2:386966/148‑1 (MQ=255)
                                                           aCATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCa                                                                                      >  2:1198/1‑148 (MQ=255)
                                                            cATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAg                                                                                     >  1:2180700/1‑148 (MQ=255)
                                                                attattACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAaca                                                                                 <  1:2683566/148‑1 (MQ=255)
                                                                 ttattaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAacac                                                                                >  2:1141468/1‑148 (MQ=255)
                                                                  tattaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACAcc                                                                               <  1:1212342/148‑1 (MQ=255)
                                                                  tattaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACAcc                                                                               <  2:2506323/148‑1 (MQ=255)
                                                                   attaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCa                                                                              <  1:2096869/148‑1 (MQ=255)
                                                                   attaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCa                                                                              >  2:2186205/1‑148 (MQ=255)
                                                                   attaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCa                                                                              <  2:231228/148‑1 (MQ=255)
                                                                    ttaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCAt                                                                             >  2:190447/1‑148 (MQ=255)
                                                                     taCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCAtt                                                                            >  2:1026451/1‑148 (MQ=255)
                                                                      aCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTa                                                                           <  1:2899892/148‑1 (MQ=255)
                                                                      aCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTa                                                                           >  1:719409/1‑148 (MQ=255)
                                                                      aCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTa                                                                           <  1:2761416/148‑1 (MQ=255)
                                                                       ccTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTaa                                                                          >  2:1386046/1‑148 (MQ=255)
                                                                         ttttGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACg                                                                        >  2:54377/1‑148 (MQ=255)
                                                                          tttGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGc                                                                       >  2:1459462/1‑148 (MQ=255)
                                                                          tttGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGc                                                                       <  1:1381763/148‑1 (MQ=255)
                                                                              aagaaATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGa                                                                   <  2:611929/148‑1 (MQ=255)
                                                                                   aTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAac                                                              >  2:130599/1‑148 (MQ=255)
                                                                                   aTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAac                                                              >  1:992262/1‑148 (MQ=255)
                                                                                    ttAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAacc                                                             <  2:844272/148‑1 (MQ=255)
                                                                                     tAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAacca                                                            >  1:1006992/1‑148 (MQ=255)
                                                                                      aGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAaccac                                                           >  2:632256/1‑148 (MQ=255)
                                                                                       gCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAaccacc                                                          >  1:433244/1‑148 (MQ=255)
                                                                                       gCCCTTGGGCAGAATGGTCCGAAAAAAAATATAAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAAGACTGCTTAAGATAAAATTGTAACAAACCAGAACACCATTAACGCATGATTAaccacc                                                          >  2:1423002/1‑148 (MQ=255)
                                                                                        cccTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCAccc                                                         >  2:1464109/1‑148 (MQ=255)
                                                                                         ccTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCa                                                        <  1:1617353/148‑1 (MQ=255)
                                                                                          cTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCAt                                                       <  1:663803/148‑1 (MQ=255)
                                                                                            tGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAg                                                     >  2:1618589/1‑148 (MQ=255)
                                                                                              ggCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCa                                                   >  1:110461/1‑148 (MQ=255)
                                                                                               gCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCaa                                                  <  1:826321/148‑1 (MQ=255)
                                                                                                   aaTTGTCCTAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGAt                                              <  2:2779911/148‑1 (MQ=255)
                                                                                                   aaTGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGAt                                              >  1:2551515/1‑148 (MQ=255)
                                                                                                     tGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTa                                            >  1:2844940/1‑148 (MQ=255)
                                                                                                       gTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACAATTAACGCATGATTAACCACCCATAGCAAGGATTACt                                          >  2:109824/1‑148 (MQ=255)
                                                                                                        tCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACtt                                         >  1:2043924/1‑148 (MQ=255)
                                                                                                         ccGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACttt                                        >  2:1227357/1‑148 (MQ=255)
                                                                                                         ccGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACttt                                        <  2:3195857/148‑1 (MQ=255)
                                                                                                            aaaaaaaaTATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGt                                     >  1:1666574/1‑148 (MQ=255)
                                                                                                                  aaTATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGaa                               >  2:2394693/1‑148 (MQ=255)
                                                                                                                    tatCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGaaaa                             >  2:576942/1‑148 (MQ=255)
                                                                                                                    tatCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGaaaa                             >  2:576777/1‑148 (MQ=255)
                                                                                                                         gagagaAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAAcaca                        >  1:920132/1‑148 (MQ=255)
                                                                                                                               aaaaCCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATg                  <  1:1467950/148‑1 (MQ=255)
                                                                                                                                aaaCCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGa                 >  1:3042977/1‑148 (MQ=255)
                                                                                                                                aaaCCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGa                 >  1:1254951/1‑148 (MQ=255)
                                                                                                                                  aCCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGaaa               >  2:1344957/1‑148 (MQ=255)
                                                                                                                                  aCCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGaaa               <  2:3102888/148‑1 (MQ=255)
                                                                                                                                    caacaaGGAAATCTCAAGATGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGa             <  2:3073511/148‑1 (MQ=255)
                                                                                                                                      acaaGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAAt           >  2:289470/1‑148 (MQ=255)
                                                                                                                                       caaGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATa          >  1:1710085/1‑148 (MQ=255)
                                                                                                                                         aGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGc        >  2:579994/1‑148 (MQ=255)
                                                                                                                                          ggAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCt       <  2:685873/148‑1 (MQ=255)
                                                                                                                                             aaTCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCttta    <  1:2983116/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATGATTAACCACCCATAGCAAGGATTACTTTTGTAAAGAAAAACACAGCAATGAAAGAATAGCTTTATT  >  CP009273/1632862‑1633152

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