Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 743,797 G→T 100% R214R (CGC→CGA)  ybgO ← putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273743,7970GT98.7% 1108.2 / 2.3 345R214R (CGC→CGA) ybgOputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (168/170);  total (172/173)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 351 total aligned reads displayed.

GGTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGACGT  >  CP009273/743650‑743945
                                                                                                                                                   |                                                                                                                                                    
ggTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCActcg                                                                                                                                                    >  2:146556/1‑150 (MQ=255)
ggTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCActcg                                                                                                                                                    >  2:3887735/1‑150 (MQ=255)
  tGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGaa                                                                                                                                                  >  1:1768336/1‑150 (MQ=255)
  tGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCAACTGTGTTTCCGGTGAAATTGAAAAAATGACGCTGCAATCCAGCACTCGaa                                                                                                                                                  >  2:4993838/1‑150 (MQ=255)
    cATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATg                                                                                                                                                <  2:3493081/150‑1 (MQ=255)
     aTTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGt                                                                                                                                               >  2:3733892/1‑150 (MQ=255)
      ttGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtt                                                                                                                                              >  2:5276142/1‑150 (MQ=255)
         aTCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTc                                                                                                                                           <  2:383472/150‑1 (MQ=255)
          tCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTcc                                                                                                                                          <  2:3981790/150‑1 (MQ=255)
            gTTTTGGCTCTTGGTCGTGTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCACCCAATTAAAATTACCAAAAACAAACACCAGTGTTTCCGGGGAAATGGAAAAAATGACGCAGCCAACCCGCACGCGAATGTTTTCCAg                                                                                                                                        >  1:4694093/1‑150 (MQ=255)
               ttGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGAcc                                                                                                                                     <  1:3052310/150‑1 (MQ=255)
                   ctctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTa                                                                                                                                 >  2:3971297/1‑150 (MQ=255)
                     ctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAat                                                                                                                               <  2:1577651/150‑1 (MQ=255)
                     ctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAat                                                                                                                               <  2:1713844/150‑1 (MQ=255)
                     ctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAat                                                                                                                               >  2:2541008/1‑150 (MQ=255)
                     ctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAat                                                                                                                               <  1:1553610/150‑1 (MQ=255)
                      ttGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAAAGGAAAAATTTACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAata                                                                                                                              >  1:1665384/1‑150 (MQ=255)
                      ttGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAACAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAata                                                                                                                              <  1:432219/150‑1 (MQ=255)
                       tGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAatat                                                                                                                             <  1:776033/150‑1 (MQ=255)
                       tGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAatat                                                                                                                             <  2:5218030/150‑1 (MQ=255)
                         gTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGa                                                                                                                           >  2:148278/1‑150 (MQ=255)
                          tCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGat                                                                                                                          <  1:608476/150‑1 (MQ=255)
                           cGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGata                                                                                                                         <  2:2838653/150‑1 (MQ=255)
                            gTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGatat                                                                                                                        >  1:511581/1‑150 (MQ=255)
                             ttttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATAtt                                                                                                                       <  2:1135027/150‑1 (MQ=255)
                              tttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCTGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATAttt                                                                                                                      <  1:2354671/150‑1 (MQ=255)
                              tttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATAttt                                                                                                                      <  2:5186336/150‑1 (MQ=255)
                                    aCTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtt                                                                                                                <  2:4790894/150‑1 (MQ=255)
                                     cTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGttt                                                                                                               >  1:2063369/1‑150 (MQ=255)
                                       gAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                             >  1:1216514/1‑150 (MQ=255)
                                       gAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                             <  1:580394/150‑1 (MQ=255)
                                          cGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGc                                                                                                          >  2:3459692/1‑150 (MQ=255)
                                                  aCATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGt                                                                                                  >  2:1167443/1‑150 (MQ=255)
                                                    aTTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCt                                                                                                >  2:3895513/1‑150 (MQ=255)
                                                    aTTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCt                                                                                                >  1:2796025/1‑150 (MQ=255)
                                                     ttGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCtt                                                                                               <  2:337866/150‑1 (MQ=255)
                                                        tgtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGt                                                                                            >  1:3943596/1‑150 (MQ=255)
                                                         gtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGtt                                                                                           >  1:937505/1‑150 (MQ=255)
                                                         gtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGtt                                                                                           <  1:491811/150‑1 (MQ=255)
                                                           ggCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATAGACGCTGCAATCCAGCACTCTAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTACGTTaa                                                                                         >  1:4925536/1‑150 (MQ=255)
                                                             cGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACg                                                                                       >  1:4359542/1‑150 (MQ=255)
                                                               aCGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCt                                                                                     <  1:2661900/150‑1 (MQ=255)
                                                               aCGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCt                                                                                     >  1:3888232/1‑150 (MQ=255)
                                                               aCGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCt                                                                                     >  1:4955336/1‑150 (MQ=255)
                                                                cGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                    >  2:2229981/1‑150 (MQ=255)
                                                                 ctataTCCCGCCAATAAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATAGACGCTGCAATCCAGCACTCGAATCATATCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCtcc                                                                                   >  2:5266071/2‑148 (MQ=255)
                                                                    tatCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCgg                                                                                >  2:2043428/1‑150 (MQ=255)
                                                                      tCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCa                                                                              >  2:2690015/1‑150 (MQ=255)
                                                                        cAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCAcc                                                                            <  1:4600146/150‑1 (MQ=255)
                                                                         aGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGTGCAGTAGGGTCTTCGTTAACGCTACCGGCACCg                                                                           >  1:3624330/1‑150 (MQ=255)
                                                                            aaattaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCa                                                                        <  1:455517/150‑1 (MQ=255)
                                                                             aattaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCaa                                                                       >  2:440095/1‑150 (MQ=255)
                                                                                taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAAtt                                                                    >  2:2917419/1‑150 (MQ=255)
                                                                                 aaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTg                                                                   <  1:186520/150‑1 (MQ=255)
                                                                                 aaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTg                                                                   <  1:3509488/150‑1 (MQ=255)
                                                                                  aattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGc                                                                  >  2:1857875/1‑150 (MQ=255)
                                                                                  aattTACC‑AAATCAATCACCTGTGTTTCCGGGGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTATTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACTCTACCGGCACCGTCAAATTGCa                                                                 >  1:1591243/1‑150 (MQ=255)
                                                                                   attTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCa                                                                 >  1:1137951/1‑150 (MQ=255)
                                                                                     ttACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacg                                                               >  1:4971223/1‑150 (MQ=255)
                                                                                          aaaaTCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGat                                                          <  1:4209793/150‑1 (MQ=255)
                                                                                                aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTc                                                    <  2:3648046/150‑1 (MQ=255)
                                                                                                aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTc                                                    >  2:4744638/1‑150 (MQ=255)
                                                                                                      cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAg                                              <  1:1479151/150‑1 (MQ=255)
                                                                                                      cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAg                                              >  1:5202890/1‑150 (MQ=255)
                                                                                                       tgtgTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGc                                             <  1:1010549/150‑1 (MQ=255)
                                                                                                        gtgtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCg                                            <  2:4702303/150‑1 (MQ=255)
                                                                                                          gtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTa                                          >  2:2291209/1‑150 (MQ=255)
                                                                                                            gtCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACt                                        <  2:126914/149‑1 (MQ=255)
                                                                                                             tCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACtt                                       <  1:4110312/150‑1 (MQ=255)
                                                                                                              ccGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACttt                                      >  2:4687918/1‑150 (MQ=255)
                                                                                                                 gTGAAATGGAAAAATTGACGCTGCAATCCAGAACGCGAATGATTTACAGACCAGTAATATGATATGTCAGGTGTTGGTCAGAAGGGTCTTCTATAACGCTACCGGCACCGTCAAATTGCACGAAGATATAGTCGCTAAGCTTAcaaagaa                                   >  1:4397529/1‑144 (MQ=255)
                                                                                                                  tGAAATTGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATa                                  <  2:4260563/150‑1 (MQ=255)
                                                                                                                  tGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATa                                  >  1:3227477/1‑150 (MQ=255)
                                                                                                                  tGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATa                                  >  2:4622692/1‑150 (MQ=255)
                                                                                                                  tGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATa                                  <  1:3834517/150‑1 (MQ=255)
                                                                                                                   gAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATaa                                 <  1:828252/150‑1 (MQ=255)
                                                                                                                       tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCAt                             <  1:3642161/150‑1 (MQ=255)
                                                                                                                       tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCAt                             <  2:4199803/150‑1 (MQ=255)
                                                                                                                        ggAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATg                            <  1:4510435/150‑1 (MQ=255)
                                                                                                                         gAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgc                           <  2:213218/150‑1 (MQ=255)
                                                                                                                          aaaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcg                          <  1:2457720/150‑1 (MQ=255)
                                                                                                                             aaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcggcg                       >  1:1423440/1‑150 (MQ=255)
                                                                                                                             aaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcggcg                       <  2:1387897/150‑1 (MQ=255)
                                                                                                                              aTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGcggcgg                      >  1:2230957/1‑150 (MQ=255)
                                                                                                                              aTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGcggcgg                      <  2:4943410/150‑1 (MQ=255)
                                                                                                                                 gACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                   >  1:516066/1‑150 (MQ=255)
                                                                                                                                  aCGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAtt                  <  2:31725/150‑1 (MQ=255)
                                                                                                                                   cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAttt                 >  2:853028/1‑150 (MQ=255)
                                                                                                                                    gCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtc                <  1:4085013/150‑1 (MQ=255)
                                                                                                                                     cTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtct               <  2:195889/150‑1 (MQ=255)
                                                                                                                                         aaTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa           >  1:2937650/1‑150 (MQ=255)
                                                                                                                                         aaTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa           >  1:34440/1‑150 (MQ=255)
                                                                                                                                         aaTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa           >  2:2945008/1‑150 (MQ=255)
                                                                                                                                         aaTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa           <  1:5276142/150‑1 (MQ=255)
                                                                                                                                          aTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAt          >  1:3562605/1‑150 (MQ=255)
                                                                                                                                             cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtttt       <  1:3733892/150‑1 (MQ=255)
                                                                                                                                             cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtttt       <  2:3240998/150‑1 (MQ=255)
                                                                                                                                              aGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTg      <  1:2558314/150‑1 (MQ=255)
                                                                                                                                              aGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTg      <  2:3045514/150‑1 (MQ=255)
                                                                                                                                                cACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGAc    <  1:3328564/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGACGT  >  CP009273/743650‑743945

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