Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,154,051 G→C G162A (GGA→GCA)  btuB → vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,154,0510GC99.1% 482.6 / ‑4.1 111G162A (GGA→GCA) btuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (57/53);  minor base T (1/0);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

Only 100 of 121 total aligned reads displayed.

TATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAT  >  CP009273/4153953‑4154148
                                                                                                  |                                                                                                 
tatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcagcag                                                                                                 >  1:2155930/1‑101 (MQ=255)
tatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGCGGTGAATATCATCACGCCGCGCGATGAACCCGGAACGGAAATTTcagcag                                                                                                 >  1:244969/1‑101 (MQ=255)
tagccGTGGGCCGCGCTGCGCTGTTTATGGTTTCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCGGGAACGGAAATTTcagcag                                                                                                 <  2:303398/98‑1 (MQ=255)
 atCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGAAAATTTCagcagg                                                                                                >  2:3564423/1‑101 (MQ=255)
   ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGt                                                                                              >  2:1133273/1‑101 (MQ=255)
   ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGt                                                                                              <  1:788671/101‑1 (MQ=255)
    cGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTg                                                                                             >  2:1753312/1‑101 (MQ=255)
    cGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTg                                                                                             >  1:2951040/1‑101 (MQ=255)
     gTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAgggggg                                                                                            >  2:126491/1‑101 (MQ=255)
       gggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGAAAATTTCAGCAGGGTgggg                                                                                          <  2:352178/101‑1 (MQ=255)
        ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGa                                                                                         <  2:1929819/101‑1 (MQ=255)
         gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGaa                                                                                        >  1:3550362/1‑101 (MQ=255)
         gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGaa                                                                                        <  2:2855480/101‑1 (MQ=255)
          ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAg                                                                                       >  1:3905417/1‑101 (MQ=255)
          ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAg                                                                                       <  1:2191132/101‑1 (MQ=255)
           cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGc                                                                                      >  1:549392/1‑101 (MQ=255)
           cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGc                                                                                      <  2:1252735/101‑1 (MQ=255)
           cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGc                                                                                      <  1:342625/101‑1 (MQ=255)
           cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGc                                                                                      >  1:2342289/1‑101 (MQ=255)
             cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCaa                                                                                    >  2:2520318/1‑101 (MQ=255)
             cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCaa                                                                                    >  1:1044953/1‑101 (MQ=255)
               cTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATa                                                                                  <  2:1053232/101‑1 (MQ=255)
                tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAg                                                                                 <  1:1940349/101‑1 (MQ=255)
                tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGGGGGGAAGCAATAg                                                                                 >  1:681608/1‑101 (MQ=255)
                  cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGtt                                                                               >  1:691667/1‑101 (MQ=255)
                   gCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTa                                                                              <  1:1810301/101‑1 (MQ=255)
                     tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATc                                                                            >  1:2024446/1‑101 (MQ=255)
                      gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCa                                                                           >  1:2892083/1‑101 (MQ=255)
                       tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAg                                                                          <  2:3136759/101‑1 (MQ=255)
                        ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGa                                                                         >  1:412451/1‑101 (MQ=255)
                        ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGa                                                                         <  1:1936796/101‑1 (MQ=255)
                          aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAAc                                                                       <  2:2745174/101‑1 (MQ=255)
                          aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGGGGGGAAGCAATAGTTATCAGAAc                                                                       >  2:1568742/1‑101 (MQ=255)
                           tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACt                                                                      <  2:1068235/101‑1 (MQ=255)
                             gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                                    <  1:336560/101‑1 (MQ=255)
                              ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                                   >  1:213340/1‑101 (MQ=255)
                                ccGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                                 >  1:2117819/1‑101 (MQ=255)
                                 cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                                >  2:3231874/1‑101 (MQ=255)
                                 cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                                >  1:1148104/1‑101 (MQ=255)
                                  gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                                               <  2:1850356/101‑1 (MQ=255)
                                   aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                              >  1:456830/1‑101 (MQ=255)
                                    tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtct                                                             >  2:700384/1‑101 (MQ=255)
                                     gCAATAGGCGGGNGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGCGGGGANGCNATAGCTATCAGAACTATGATGtctc                                                            >  2:1430351/1‑101 (MQ=255)
                                     gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                                            <  2:2558570/101‑1 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          >  1:636330/1‑101 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          <  2:3848981/101‑1 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          <  2:1416390/101‑1 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          >  1:2323849/1‑101 (MQ=255)
                                         tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                        <  1:30787/101‑1 (MQ=255)
                                          aGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                                       >  1:181538/1‑101 (MQ=255)
                                          aGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                                       <  2:1342110/101‑1 (MQ=255)
                                           ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgca                                                      <  2:1743943/101‑1 (MQ=255)
                                             cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                    <  2:3236013/101‑1 (MQ=255)
                                                 gtggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                                >  1:2580336/1‑101 (MQ=255)
                                                  tggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                               >  2:3687016/1‑101 (MQ=255)
                                                   ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                              >  2:1917868/1‑101 (MQ=255)
                                                     tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                            >  1:1471147/1‑101 (MQ=255)
                                                     tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                            <  1:1981262/101‑1 (MQ=255)
                                                        atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                         <  2:563225/101‑1 (MQ=255)
                                                         tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                        >  1:3108370/1‑101 (MQ=255)
                                                         tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                        <  2:2956359/101‑1 (MQ=255)
                                                         tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                        >  1:3258227/1‑101 (MQ=255)
                                                          atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                                       >  2:730144/1‑101 (MQ=255)
                                                          atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                                       <  2:142870/101‑1 (MQ=255)
                                                            catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAATCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                                     <  1:1130353/101‑1 (MQ=255)
                                                            catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                                     >  2:2950496/1‑101 (MQ=255)
                                                            catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                                     <  1:720078/101‑1 (MQ=255)
                                                                 cgacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca                                >  2:1129621/1‑101 (MQ=255)
                                                                 cgacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca                                >  1:3207026/1‑101 (MQ=255)
                                                                 cgacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca                                <  2:3158012/101‑1 (MQ=255)
                                                                  gacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                               >  1:3369080/1‑101 (MQ=255)
                                                                  gacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                               >  1:3279882/1‑101 (MQ=255)
                                                                  gacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                               >  2:2459993/1‑101 (MQ=255)
                                                                   gcgcgcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                              <  2:2690592/100‑1 (MQ=255)
                                                                    cgcgcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                             <  1:2060959/101‑1 (MQ=255)
                                                                       gcgATGAACCCGGAACGGAAATTTCAGTAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                          >  1:61553/1‑101 (MQ=255)
                                                                       gcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                          <  1:1734192/101‑1 (MQ=255)
                                                                       gcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                          >  1:157904/1‑101 (MQ=255)
                                                                       gcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                          <  1:3370042/101‑1 (MQ=255)
                                                                        cgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                         <  2:3026529/101‑1 (MQ=255)
                                                                         gATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                        <  1:285088/101‑1 (MQ=255)
                                                                           tGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc                      >  1:747652/1‑101 (MQ=255)
                                                                           tGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc                      <  2:2760281/101‑1 (MQ=255)
                                                                           tGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc                      >  1:1066165/1‑101 (MQ=255)
                                                                               cccGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt                  <  2:2138855/101‑1 (MQ=255)
                                                                               cccGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt                  <  1:3373971/101‑1 (MQ=255)
                                                                                ccGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg                 >  1:603297/1‑101 (MQ=255)
                                                                                 cggaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTgg                <  1:1090253/101‑1 (MQ=255)
                                                                                 cggaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTgg                <  1:781162/101‑1 (MQ=255)
                                                                                    aacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg             >  2:2703471/1‑101 (MQ=255)
                                                                                     acggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGa            >  1:920742/1‑101 (MQ=255)
                                                                                     acggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGa            >  1:3320114/1‑101 (MQ=255)
                                                                                     acggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGa            >  2:555977/1‑101 (MQ=255)
                                                                                      cggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt           <  2:1200247/101‑1 (MQ=255)
                                                                                       ggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt          <  2:3238239/101‑1 (MQ=255)
                                                                                         aaaTTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAt        <  2:3499198/101‑1 (MQ=255)
                                                                                         aaaTTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAt        <  1:1969991/101‑1 (MQ=255)
                                                                                         aaaTTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAt        >  2:1562377/1‑101 (MQ=255)
                                                                                          aaTTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATg       >  1:607636/1‑101 (MQ=255)
                                                                                               cagcagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAt  <  2:2307737/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAT  >  CP009273/4153953‑4154148

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