Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,058,833 C→T 21.2% G138E (GGA→GAA)  cobU ← bifunctional cobinamide kinase/ cobinamide phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,058,8330CT21.2% 577.3 / 147.8 250G138E (GGA→GAA) cobUbifunctional cobinamide kinase/ cobinamide phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (103/94);  new base T (24/29);  total (127/123)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 262 total aligned reads displayed.

GCCCAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTA  >  CP009273/2058684‑2058982
                                                                                                                                                     |                                                                                                                                                     
gCCCAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTc                                                                                                                                                       >  1:2416110/1‑150 (MQ=255)
gCCCAAAATAATTTACTCAATTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTc                                                                                                                                                       <  2:3634774/150‑1 (MQ=255)
 cccAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTcc                                                                                                                                                      <  1:66977/150‑1 (MQ=255)
  ccAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTtcc                                                                                                                                                     >  1:1063726/1‑150 (MQ=255)
  ccAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTtcc                                                                                                                                                     >  2:4743480/1‑150 (MQ=255)
   cAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCa                                                                                                                                                    <  1:601184/150‑1 (MQ=255)
   cAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCa                                                                                                                                                    >  1:3493048/1‑150 (MQ=255)
     aaaTAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACt                                                                                                                                                  >  1:1407487/1‑150 (MQ=255)
      aataatTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                                                                                                                 <  1:1763240/150‑1 (MQ=255)
       ataatTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTc                                                                                                                                                >  2:1066980/1‑150 (MQ=255)
       ataatTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTc                                                                                                                                                >  1:1309721/1‑150 (MQ=255)
        taatTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATGTCCACTTCg                                                                                                                                               >  2:3819420/1‑150 (MQ=255)
          atTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGtt                                                                                                                                             <  2:1680938/150‑1 (MQ=255)
           tttACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTa                                                                                                                                            <  1:205240/150‑1 (MQ=255)
           tttACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTa                                                                                                                                            >  2:292810/1‑150 (MQ=255)
            ttACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAg                                                                                                                                           <  1:3767918/150‑1 (MQ=255)
                     tGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATa                                                                                                                                  <  1:222028/150‑1 (MQ=255)
                     tGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATa                                                                                                                                  <  1:3291020/150‑1 (MQ=255)
                       aTTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACa                                                                                                                                <  1:1592995/150‑1 (MQ=255)
                       aTTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACa                                                                                                                                <  1:3199399/150‑1 (MQ=255)
                       aTTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACa                                                                                                                                <  1:1325660/150‑1 (MQ=255)
                            tACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCtt                                                                                                                           <  1:4750407/150‑1 (MQ=255)
                              cTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTg                                                                                                                         <  1:2918852/150‑1 (MQ=255)
                               tCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGAGTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCAGTCCCACTTCGTTAGTCACTAATACAACCTTTGc                                                                                                                        >  1:1284639/1‑150 (MQ=255)
                                ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCg                                                                                                                       >  2:3745257/1‑150 (MQ=255)
                                ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCg                                                                                                                       >  2:1447939/1‑150 (MQ=255)
                                ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGAGTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCCCTAATACAACCTTTGCg                                                                                                                       >  2:1309178/1‑150 (MQ=255)
                                  aaTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCggg                                                                                                                     >  2:1683944/1‑150 (MQ=255)
                                  aaTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCggg                                                                                                                     >  2:4342576/1‑150 (MQ=255)
                                           aaCCCCCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTg                                                                                                            <  2:415803/150‑1 (MQ=255)
                                           aaCCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTg                                                                                                            <  2:2849927/150‑1 (MQ=255)
                                            accaccAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                                                                                                           >  1:4358440/1‑150 (MQ=255)
                                             ccaccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggc                                                                                                          <  2:4788923/150‑1 (MQ=255)
                                             ccaccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggc                                                                                                          >  1:4664384/1‑150 (MQ=255)
                                               accaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggcag                                                                                                        >  2:2598409/1‑150 (MQ=255)
                                                   gCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGGCACTAATACAACCTTTGCGGGGCAACGTTGgcaggcag                                                                                                    >  1:2174251/1‑150 (MQ=255)
                                                    ccATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagc                                                                                                   <  1:646324/150‑1 (MQ=255)
                                                      aTACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCaa                                                                                                 >  1:1965696/1‑150 (MQ=255)
                                                         cTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCa                                                                                              <  1:3404869/150‑1 (MQ=255)
                                                          ttCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCaa                                                                                             >  1:4845784/1‑150 (MQ=255)
                                                                tGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTg                                                                                       >  1:3130265/1‑150 (MQ=255)
                                                                tGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTg                                                                                       >  1:1749985/1‑150 (MQ=255)
                                                                    gCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAAtc                                                                                   <  2:1011047/150‑1 (MQ=255)
                                                                     cAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAAtct                                                                                  <  1:965645/150‑1 (MQ=255)
                                                                      aGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAAtctc                                                                                 <  2:4639294/150‑1 (MQ=255)
                                                                       gCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCa                                                                                >  2:1137599/1‑150 (MQ=255)
                                                                         ggCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGACTCTCAGc                                                                              >  1:2654018/1‑150 (MQ=255)
                                                                             aaCTGCTGATTTACCCGCCCGGCAATATCACGAAATTGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTa                                                                          <  2:4658811/150‑1 (MQ=255)
                                                                             aaCTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTa                                                                          <  2:2205877/150‑1 (MQ=255)
                                                                             aaCTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTa                                                                          >  1:1225239/1‑150 (MQ=255)
                                                                              aCTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTaa                                                                         >  2:3859941/1‑150 (MQ=255)
                                                                               ctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAAt                                                                        <  2:3000277/150‑1 (MQ=255)
                                                                               ctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAAt                                                                        <  1:281585/150‑1 (MQ=255)
                                                                               ctgctgATTTACCCGCCCGGAAATATCACGAAATTGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGAAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAAt                                                                        <  1:1838896/150‑1 (MQ=255)
                                                                                 gctgATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCg                                                                      >  2:704645/1‑150 (MQ=255)
                                                                                    gATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCt                                                                   >  2:1167795/1‑150 (MQ=255)
                                                                                     aTTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTg                                                                  >  2:1134758/1‑150 (MQ=255)
                                                                                        tACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTc                                                               >  2:2945850/1‑150 (MQ=255)
                                                                                        tACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTc                                                               <  2:3044436/150‑1 (MQ=255)
                                                                                         aCCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTcc                                                              <  2:2416600/150‑1 (MQ=255)
                                                                                         aCCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTcc                                                              <  2:3130265/150‑1 (MQ=255)
                                                                                         aCCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTcc                                                              <  2:1749985/150‑1 (MQ=255)
                                                                                         aCCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTcc                                                              <  2:1750407/150‑1 (MQ=255)
                                                                                           ccgcccgGCAATATCACTAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCAt                                                            <  1:822681/150‑1 (MQ=255)
                                                                                                cgGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCt                                                       <  2:1483878/150‑1 (MQ=255)
                                                                                                     atatCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATaa                                                  >  1:2146486/1‑150 (MQ=255)
                                                                                                      tataACTAAATTGTCTTGAAAGACTACTCTCCGGAAAAATCCCCATTCCAAATTCGTTTGTCACTAATAAAACCTTTGCGGGGCAAATTTTGCAGGAAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGACTGATAAt                                                 <  1:4588399/150‑1 (MQ=255)
                                                                                                      tatCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAAt                                                 <  2:4358440/150‑1 (MQ=255)
                                                                                                       atCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATc                                                <  1:3005489/150‑1 (MQ=255)
                                                                                                       atCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATc                                                >  1:255765/1‑150 (MQ=255)
                                                                                                         cACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATccc                                              >  1:1199856/1‑150 (MQ=255)
                                                                                                            gAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCAtt                                           >  2:2122235/1‑150 (MQ=255)
                                                                                                            gAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCAtt                                           <  1:356108/150‑1 (MQ=255)
                                                                                                             aaaaTGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAATCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATtc                                          <  2:3405836/150‑1 (MQ=255)
                                                                                                             aaaaTGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATtc                                          <  1:3865747/150‑1 (MQ=255)
                                                                                                             aaaaTGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATtc                                          >  1:1874030/1‑150 (MQ=255)
                                                                                                               aaTGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATtcat                                        <  2:55904/150‑1 (MQ=255)
                                                                                                                 tGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATtcatca                                      >  2:3281978/1‑150 (MQ=255)
                                                                                                                  gTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAg                                     <  2:2380594/150‑1 (MQ=255)
                                                                                                                      tGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTc                                 <  2:1991263/150‑1 (MQ=255)
                                                                                                                       gCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCt                                <  1:3032565/150‑1 (MQ=255)
                                                                                                                          aGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTa                             <  1:1134865/150‑1 (MQ=255)
                                                                                                                            acgacTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATc                           <  1:738750/150‑1 (MQ=255)
                                                                                                                             cgacTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATcg                          >  1:4181877/1‑150 (MQ=255)
                                                                                                                              gacTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATcgc                         >  1:166880/1‑150 (MQ=255)
                                                                                                                              gacTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATcgc                         >  2:2102634/1‑150 (MQ=255)
                                                                                                                                 tctcCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATcgccgc                      >  1:110935/1‑150 (MQ=255)
                                                                                                                                     cGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATa                  <  2:1632965/150‑1 (MQ=255)
                                                                                                                                      ggCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATaa                 <  2:453419/150‑1 (MQ=255)
                                                                                                                                      ggCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCAACGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATaa                 >  2:2753884/1‑150 (MQ=255)
                                                                                                                                        caaaaTCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCTACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTCTCGCCGCCATAATc               <  2:4235042/147‑1 (MQ=255)
                                                                                                                                          caATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCaa             <  2:1225239/150‑1 (MQ=255)
                                                                                                                                          caATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCaa             >  1:4016375/1‑150 (MQ=255)
                                                                                                                                          caATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCaa             <  2:3493048/150‑1 (MQ=255)
                                                                                                                                           aaTCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCaaa            >  1:645713/1‑150 (MQ=255)
                                                                                                                                           aaTCCCCATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCaaa            >  1:1264959/1‑150 (MQ=255)
                                                                                                                                               cccATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAaacaac        <  1:4793749/150‑1 (MQ=255)
                                                                                                                                                ccATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGAAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAacaaca       <  1:3433995/150‑1 (MQ=255)
                                                                                                                                                 cATTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAg      >  2:552017/1‑150 (MQ=255)
                                                                                                                                                  aTTTCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGa     <  1:362352/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCAAAATAATTTACTCATTTGATTTTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGACTGAATCTCAGCATTAATCGCCTGTTCCATCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTA  >  CP009273/2058684‑2058982

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: