Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,896 C→T 6.5% pseudogene (1248/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,460,8960CT 6.5% 672.7 / 23.1 208pseudogene (1248/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (111/81);  new base T (7/6);  total (119/89)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

CCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACC  >  CP009273/1460748‑1461043
                                                                                                                                                    |                                                                                                                                                   
ccGGCGACAGCGCTACGGTGGCTAAGAAGGGCACAAAGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                                    <  2:4219195/150‑1 (MQ=18)
ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                                    <  2:1942828/150‑1 (MQ=255)
  ggtgACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                                                                                                                                                  <  1:1308862/147‑1 (MQ=255)
   gCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtg                                                                                                                                                 <  1:2764025/150‑1 (MQ=255)
    cGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATGCGGAGTGCATCGGTATCCAGGTTGAGGGCGACCTGTCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                                <  1:3312965/150‑1 (MQ=18)
     gACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgat                                                                                                                                               >  2:3921967/1‑150 (MQ=255)
         gcgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACg                                                                                                                                           <  1:3460629/150‑1 (MQ=255)
            cAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                        <  2:2366061/150‑1 (MQ=255)
            cAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                        >  1:3555967/1‑150 (MQ=255)
            cAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                        <  2:3017463/150‑1 (MQ=255)
            cAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                        >  2:1329972/1‑150 (MQ=255)
              aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACg                                                                                                                                      <  1:3646080/150‑1 (MQ=255)
              aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACg                                                                                                                                      <  1:1513682/150‑1 (MQ=255)
                       aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGGGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                             <  1:1088480/150‑1 (MQ=255)
                       aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                             >  2:3005824/1‑150 (MQ=255)
                             ggCACCTTGACCGTCACCGTTCCGGAGTCCATCGGTATCGAGGTTGACGTCGACCAGGCGGTTGTTATCAAAGAAGGCGAGAGCGCCATCACCAACGGTGGCACCTGCACGCAGATTAACGGTTATGACGCCACGGCAAACAACATCGGc                                                                                                                       <  2:3438065/150‑1 (MQ=255)
                                    tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                                                                                                                <  2:2728649/150‑1 (MQ=255)
                                    tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAacca                                                                                                                >  1:924674/1‑150 (MQ=255)
                                         gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                           >  1:2366063/1‑150 (MQ=255)
                                         gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                                                           <  1:2893078/150‑1 (MQ=255)
                                          tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                          >  1:1674858/1‑150 (MQ=255)
                                           cACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGACGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAACACCACCGTTGa                                                                                                         >  2:4028246/1‑150 (MQ=255)
                                           cACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                         >  1:259788/1‑150 (MQ=255)
                                           cACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGAGAGAGCGCCATCACCAACGGTGGCACCGGCAAGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                         >  2:2523525/1‑150 (MQ=255)
                                            aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTAAAAAACGAAGGCGAGAGCGCCATCACCCACGTTCGAACCGGCACGCAGATTAACGGTGATGACGCAACGGCAAACAACAACGGCAAAACCACAGATGAc                                                                                                        >  2:4111327/1‑150 (MQ=17)
                                                    ggAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGa                                                                                                >  1:3174598/1‑150 (MQ=255)
                                                      aGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                              >  1:1334768/1‑150 (MQ=255)
                                                       gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                             <  2:575047/150‑1 (MQ=255)
                                                        tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                                                            <  1:2056339/150‑1 (MQ=255)
                                                                  tCCAGGTTGACGGCGACCAGGCGTTAGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                  <  2:1598176/150‑1 (MQ=255)
                                                                    cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAGGACGCCACGGCAAACAACAACGGAAAAACCACCGTTGAAGGCAAGGATTCCACCGGTAcccaa                                                                                >  1:697334/1‑146 (MQ=18)
                                                                          gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                          >  1:1457905/1‑150 (MQ=255)
                                                                           aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                         <  1:3872650/150‑1 (MQ=255)
                                                                                 aCCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAATAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATa                                                                   <  2:859723/150‑1 (MQ=255)
                                                                                 aCCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCAAAAAACAGAGGCAAAACCACAGTTGAGGGAAAGGATACAACCGAAGCGCAAATTGATGGCAAtc                                                                   >  2:4453882/1‑149 (MQ=17)
                                                                                  ccAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATaa                                                                  >  1:1221322/1‑150 (MQ=255)
                                                                                    aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                                                                <  2:4390234/150‑1 (MQ=255)
                                                                                    aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                                                                >  2:4289147/1‑150 (MQ=255)
                                                                                    aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCGAACAACAACGGAGATACCACCGTAGACGGCAAGGTTAAGACCGGGTCGGAAAGTGCTGTCatatcgg                                                                >  2:2646243/1‑144 (MQ=17)
                                                                                    aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTCAAGACGCCACCGCAAACAACACCGGCAAAACCACCGTAGTCGGCAAGGATTCCACCGGTACGGACATTGCTGGCaaatacg                                                                >  2:2646226/1‑145 (MQ=17)
                                                                                     ggCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                               >  2:2451537/1‑150 (MQ=255)
                                                                                     ggCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                               >  2:1976697/1‑150 (MQ=255)
                                                                                       cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGa                                                             <  2:1384933/150‑1 (MQ=255)
                                                                                       cGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGa                                                             >  2:3674591/1‑150 (MQ=255)
                                                                                        ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                            >  1:3015638/1‑150 (MQ=255)
                                                                                         gttgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAg                                                           <  1:3811536/150‑1 (MQ=255)
                                                                                            gttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                                                        <  1:4477072/150‑1 (MQ=255)
                                                                                             ttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGCTGACGCCACGGCAAACACCAACGGCCCAACCACCGTTGACGGCCAGGAGGCCACCGGTACGGACATTGCTGGCACTCCCGGGACGGGGa                                                       >  2:2522889/1‑150 (MQ=14)
                                                                                              tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCCCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                                      <  2:1612189/150‑1 (MQ=255)
                                                                                              tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                                      <  2:219368/150‑1 (MQ=255)
                                                                                               aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                                     <  1:3148913/150‑1 (MQ=255)
                                                                                                acaacGATGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                    <  1:4496789/150‑1 (MQ=255)
                                                                                                acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                    <  1:247371/150‑1 (MQ=255)
                                                                                                acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                    <  2:1019305/150‑1 (MQ=255)
                                                                                                 taacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGcc                                                                                                                                         >  1:2616710/2‑64 (MQ=16)
                                                                                                 caacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                                                   <  2:60158/150‑1 (MQ=255)
                                                                                                  tctgAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                                                  <  1:3052308/147‑1 (MQ=255)
                                                                                                  aacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTCATGACGCCACGGCGAAAAACAAAGGCAAAACCACCGTTGAAGGAAATGATTCCACCGGTACCAAAATAGCTGGCAATAACGGTAAGTTGAtactg                                                  >  2:3379207/1‑146 (MQ=14)
                                                                                                    cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                                <  1:4076296/150‑1 (MQ=255)
                                                                                                    cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                                >  2:4037338/1‑150 (MQ=255)
                                                                                                     gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                                                  >  2:2302729/1‑147 (MQ=255)
                                                                                                      aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATTAACGGTGATGACGCCACGGCAAAAAACATAGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAATGTGATTCTGGACg                                              >  1:4237366/1‑150 (MQ=17)
                                                                                                      aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg                                              >  1:1620481/1‑150 (MQ=255)
                                                                                                      aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg                                              <  1:2872440/150‑1 (MQ=255)
                                                                                                      aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACGGCAAACAACAACGGCAAAACAACCGATGACGGCAACGATTCCACCGGTACGGAAATTCCTTGCAATAACGCGAAGGTGATTCATGACg                                              >  2:2048404/1‑150 (MQ=17)
                                                                                                      aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAGGACGCCACGGCAAACAACAACGTCAAAACCACCGTGGACGGCAAGGATTCCACGGGTACGGAAATTGCTGGCCATAAGGGGAAGGTGATTAATGACg                                              >  2:2048302/1‑150 (MQ=17)
                                                                                                         gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                           <  1:4412491/150‑1 (MQ=255)
                                                                                                             gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTTATGACGCCACGGCATCCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt                                       <  1:2246661/150‑1 (MQ=255)
                                                                                                             gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACCACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt                                       <  1:2247223/150‑1 (MQ=255)
                                                                                                             gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt                                       <  1:3002843/150‑1 (MQ=255)
                                                                                                              agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCCAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGACGGTGATTCAGGACGGCGATCTg                                      >  2:929744/1‑150 (MQ=255)
                                                                                                               gcaccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTgg                                     <  2:2050341/147‑1 (MQ=255)
                                                                                                                cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAGCGGGAAGGTGATTCAGGACTGCGATCTGGa                                    >  1:3692148/1‑150 (MQ=255)
                                                                                                                cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                                    >  1:2180/1‑150 (MQ=255)
                                                                                                                   cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                                 >  2:2824012/1‑150 (MQ=255)
                                                                                                                    aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTc                                <  1:4031643/150‑1 (MQ=255)
                                                                                                                     tCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCa                               >  2:1232427/1‑150 (MQ=255)
                                                                                                                         cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCCAGGATTCCACCGGTACGGAAATTGCTGGCACTAACGGGGAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                           >  2:2130885/1‑150 (MQ=255)
                                                                                                                         cAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAgg                                                          <  1:2332640/119‑1 (MQ=255)
                                                                                                                           aCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCTCGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcg                         <  2:3174912/150‑1 (MQ=255)
                                                                                                                           aCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGc                                                                                                                                    <  1:515175/43‑1 (MQ=255)
                                                                                                                           aCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGc                                                                                                                                    <  1:515468/43‑1 (MQ=255)
                                                                                                                               tGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                     >  1:1950383/1‑150 (MQ=255)
                                                                                                                                ggCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                                                    <  2:1840198/86‑1 (MQ=18)
                                                                                                                                     cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg               >  1:228031/1‑150 (MQ=255)
                                                                                                                                      ggCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt              >  2:4221378/1‑150 (MQ=255)
                                                                                                                                      ggCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt              <  1:49060/150‑1 (MQ=255)
                                                                                                                                         aCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCTGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt           >  1:807564/1‑150 (MQ=255)
                                                                                                                                         aCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt           <  2:3519732/150‑1 (MQ=255)
                                                                                                                                         aCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt           >  1:1785325/1‑150 (MQ=255)
                                                                                                                                             aGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata       >  1:496502/1‑150 (MQ=255)
                                                                                                                                             aGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata       >  2:3911450/1‑150 (MQ=255)
                                                                                                                                               aTTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc     <  1:3987878/150‑1 (MQ=255)
                                                                                                                                               aTTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc     >  2:1087173/1‑150 (MQ=255)
                                                                                                                                               aTTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc     >  1:1848334/1‑150 (MQ=255)
                                                                                                                                                 tAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt              >  2:1552115/1‑139 (MQ=255)
                                                                                                                                                 tAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAc   >  1:1758592/1‑150 (MQ=255)
                                                                                                                                                 tAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAc   >  1:1746467/1‑150 (MQ=255)
                                                                                                                                                  aaCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc  >  2:2605589/1‑150 (MQ=255)
                                                                                                                                                  aaCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc  >  2:2604681/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACC  >  CP009273/1460748‑1461043

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