Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,144,158 C→A 93.2% intergenic (+253/+19) yghA → / ← exbD putative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,144,1580CA93.2% 691.8 / 38.2 231intergenic (+253/+19)yghA/exbDputative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base C (9/7);  new base A (109/106);  total (118/113)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

TGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGG  >  CP009273/3144013‑3144304
                                                                                                                                                 |                                                                                                                                                  
tGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAt                                                                                                                                                >  1:3303480/1‑150 (MQ=255)
tGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAt                                                                                                                                                >  1:2643300/1‑150 (MQ=255)
tGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAt                                                                                                                                                >  1:1773797/1‑150 (MQ=255)
 gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATc                                                                                                                                               >  1:3454342/1‑150 (MQ=255)
 gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATc                                                                                                                                               <  2:279741/150‑1 (MQ=255)
 gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATc                                                                                                                                               >  1:1983119/1‑150 (MQ=255)
    aTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAgg                                                                                                                                            >  2:2743633/1‑150 (MQ=255)
       gATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATATATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCaa                                                                                                                                         >  1:1289842/1‑150 (MQ=255)
        aTTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAAt                                                                                                                                        >  1:1892331/1‑150 (MQ=255)
        aTTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAt                                                                                                                                        <  1:1065556/150‑1 (MQ=255)
          tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTc                                                                                                                                      >  2:4526830/1‑150 (MQ=255)
                ggCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCtacttt                                                                                                                                >  2:2222100/1‑150 (MQ=255)
                   cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttact                                                                                                                             >  1:1593308/1‑150 (MQ=255)
                      aTAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCg                                                                                                                          <  1:2403568/150‑1 (MQ=255)
                      aTAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCg                                                                                                                          <  1:101249/150‑1 (MQ=255)
                      aTAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCg                                                                                                                          <  1:2089357/150‑1 (MQ=255)
                        aaGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCt                                                                                                                        <  1:4226369/150‑1 (MQ=255)
                          gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCttt                                                                                                                      <  1:3610724/150‑1 (MQ=255)
                          gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCttt                                                                                                                      >  2:3776784/1‑150 (MQ=255)
                           cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTg                                                                                                                     <  2:1133887/150‑1 (MQ=255)
                           cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTg                                                                                                                     >  1:1330456/1‑150 (MQ=255)
                            gTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTgg                                                                                                                    <  2:4140681/150‑1 (MQ=255)
                                 gcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGtt                                                                                                               <  2:4073518/150‑1 (MQ=255)
                                  cATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCGTTGGCGGttt                                                                                                              >  2:4239813/1‑150 (MQ=255)
                                    tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttct                                                                                                            <  2:1372163/150‑1 (MQ=255)
                                      aGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCACTACTACTTTAATTCGCTTTGGCGGTttcttc                                                                                                          >  2:4293292/1‑150 (MQ=255)
                                       ggCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCg                                                                                                         >  2:1414718/1‑150 (MQ=255)
                                          aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCg                                                                                                      >  1:169561/1‑150 (MQ=255)
                                          aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCg                                                                                                      >  2:2996761/1‑150 (MQ=255)
                                           ttGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGa                                                                                                     >  1:3501152/1‑150 (MQ=255)
                                           ttGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGAAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGa                                                                                                     >  2:1039171/1‑150 (MQ=255)
                                            tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAc                                                                                                    <  2:2959824/150‑1 (MQ=255)
                                             gAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAcc                                                                                                   <  2:1453179/150‑1 (MQ=255)
                                              aGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCa                                                                                                  >  2:2895739/1‑150 (MQ=255)
                                                    tAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTa                                                                                            >  2:2121255/1‑150 (MQ=255)
                                                    tAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTa                                                                                            >  1:3203620/1‑150 (MQ=255)
                                                    tAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTa                                                                                            >  1:2558277/1‑150 (MQ=255)
                                                     aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTAt                                                                                           >  2:1310206/1‑150 (MQ=255)
                                                       tGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCt                                                                                         <  2:3957463/150‑1 (MQ=255)
                                                       tGCCGGATGCGTCGTGAACGCCATATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCt                                                                                         <  2:1698998/150‑1 (MQ=255)
                                                         ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTc                                                                                       >  2:1684607/1‑150 (MQ=255)
                                                         ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTc                                                                                       <  1:1438047/150‑1 (MQ=255)
                                                           ggATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCGGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAg                                                                                     <  1:520617/150‑1 (MQ=255)
                                                           ggATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAg                                                                                     <  1:518922/150‑1 (MQ=255)
                                                               gCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTGTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTaa                                                                                 <  2:1851433/150‑1 (MQ=255)
                                                                cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                                <  2:3150676/150‑1 (MQ=255)
                                                                      gAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCt                                                                          >  1:668737/1‑150 (MQ=255)
                                                                         cGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAGTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAt                                                                       >  2:229021/1‑150 (MQ=255)
                                                                          gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATg                                                                      <  2:2179876/150‑1 (MQ=255)
                                                                          gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATg                                                                      <  2:2353187/150‑1 (MQ=255)
                                                                          gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATg                                                                      <  2:2499331/150‑1 (MQ=255)
                                                                           ccTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGc                                                                     >  1:727764/1‑150 (MQ=255)
                                                                            cTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCa                                                                    >  1:1423200/1‑150 (MQ=255)
                                                                            cTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCa                                                                    >  2:3284382/1‑150 (MQ=255)
                                                                             ttATCCGGCCTACAAACGCACGACGATTCATCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAg                                                                   <  2:493320/150‑1 (MQ=255)
                                                                               aTCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTaa                                                                                 <  1:1559701/134‑1 (MQ=255)
                                                                               aTCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTaa                                                                                 >  2:1557948/1‑134 (MQ=255)
                                                                               aTCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTaa                                                                                 >  2:1559701/1‑134 (MQ=255)
                                                                                 ccGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTa                                                               >  1:1644567/1‑150 (MQ=255)
                                                                                  cGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTAt                                                              >  2:837694/1‑150 (MQ=255)
                                                                                  cGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAGTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTAt                                                              >  1:1817355/1‑150 (MQ=255)
                                                                                  cGGCCGACAAACGCACGAAGATTCAACATGGTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTGACTTCGCTGTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTAt                                                              <  2:799986/150‑1 (MQ=37)
                                                                                    gCCTACAATCGCACGAAGATTCAACATGTTGCCATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCCGGTAACCCGCGTGATGCAGCGTAtgc                                                            >  2:3595816/1‑148 (MQ=255)
                                                                                    gCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATcc                                                            >  2:2941328/1‑150 (MQ=255)
                                                                                     ccTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCa                                                           <  2:668737/150‑1 (MQ=255)
                                                                                     ccTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCa                                                           >  1:4264301/1‑150 (MQ=255)
                                                                                     ccTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCa                                                           >  2:2829048/1‑150 (MQ=255)
                                                                                            aCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCtt                                                    >  2:3229059/1‑150 (MQ=255)
                                                                                              gCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCATACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtca                                                  <  1:3690948/150‑1 (MQ=255)
                                                                                              gCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtca                                                  >  2:3755693/1‑150 (MQ=255)
                                                                                              gCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtca                                                  <  1:4058434/150‑1 (MQ=255)
                                                                                               cACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcat                                                 <  1:3907490/150‑1 (MQ=255)
                                                                                                  gAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCaa                                              <  1:3841974/150‑1 (MQ=255)
                                                                                                    aGATTCAACATGTTGCACTAGATGTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTAGTTTACTTCGCGTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTCTCCATTACCTTCATCAACg                                            >  2:2683681/1‑150 (MQ=255)
                                                                                                    aGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACg                                            >  2:2682466/1‑150 (MQ=255)
                                                                                                      aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                          <  2:1995412/150‑1 (MQ=255)
                                                                                                      aTTCAACATGTTGCAATAGATTTTGTAGGCATGATGAGCGTATAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                          >  1:2227833/1‑150 (MQ=255)
                                                                                                       ttCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtct                                         <  2:148741/150‑1 (MQ=255)
                                                                                                        tCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtctc                                        <  1:2860930/150‑1 (MQ=255)
                                                                                                         cAACATGTTGCAATATATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCg                                       <  2:4427094/150‑1 (MQ=255)
                                                                                                          accATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGt                                      >  2:4052072/3‑150 (MQ=255)
                                                                                                              tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                  <  1:2743633/150‑1 (MQ=255)
                                                                                                               gTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCg                                 >  1:15374/1‑150 (MQ=255)
                                                                                                                  gCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACg                              <  1:4219735/150‑1 (MQ=255)
                                                                                                                    aaTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGt                            >  1:517353/1‑150 (MQ=255)
                                                                                                                      tAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGttt                          >  1:915066/1‑150 (MQ=255)
                                                                                                                       aGATTTTGTAGGACTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGTTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACTTCTCGTAATCGACGGtttt                         >  1:2036988/1‑150 (MQ=255)
                                                                                                                        gATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAATACCTTAATCAACGTCTCGTAATCGACGGTTTTg                        >  2:3673629/1‑150 (MQ=255)
                                                                                                                         aTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAATACCTTCATCAACGTCACGTAATCGACGGTTTTGt                       >  1:4421579/1‑150 (MQ=255)
                                                                                                                          ttttGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTc                      <  1:4402644/150‑1 (MQ=255)
                                                                                                                              gTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCt                  <  2:3109512/150‑1 (MQ=255)
                                                                                                                               tAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTc                 >  1:4454238/1‑150 (MQ=255)
                                                                                                                               tAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTc                 >  2:4369230/1‑150 (MQ=255)
                                                                                                                                      gatgaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaag            >  1:3385040/1‑148 (MQ=255)
                                                                                                                                      gatgaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaag            <  2:3385040/148‑1 (MQ=255)
                                                                                                                                         gaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAt       <  2:1774992/150‑1 (MQ=255)
                                                                                                                                            cGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggt    <  1:1684703/150‑1 (MQ=255)
                                                                                                                                            cGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggt    <  1:2865430/150‑1 (MQ=255)
                                                                                                                                             gTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggtg   >  2:2040590/1‑150 (MQ=255)
                                                                                                                                              tAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggtgg  <  1:2613906/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGG  >  CP009273/3144013‑3144304

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: