Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,689 A→G 6.0% pseudogene (1041/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,6890AG 6.0% 467.3 / 12.8 159pseudogene (1041/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (72/72);  new base G (4/5);  total (77/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

GCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCG  >  CP009273/1460540‑1460836
                                                                                                                                                     |                                                                                                                                                   
gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                                                                                                                                     >  2:290855/1‑150 (MQ=255)
  cATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaa                                                                                                                                                   <  1:309840/150‑1 (MQ=255)
   aTCGTCAATAACGACGGCGAGAGCACCAACACCAACGCTGGCACCGGCACCCAGATTACTGGTGATGACGCCACCGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                                                                                                                  <  2:52010/150‑1 (MQ=18)
   aTCGTCAATAACGAAGGCGAGAGCACCAACACCAACGTTGGCACCGTCAACCATATTACTGGTGATGACGCCACCGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                                                                                                                  <  2:1904730/150‑1 (MQ=18)
     cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCTGATTAATGGTGTTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                                                                                                                                                <  1:3049255/150‑1 (MQ=255)
     cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGt                                                                                                                                                <  1:2846615/150‑1 (MQ=255)
      gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                                                                                                                                               <  2:15319/150‑1 (MQ=255)
        cAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTCATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                                                                                                                                             <  2:1993646/150‑1 (MQ=255)
          tttaCGAAGGCGAGAGCACCATCACCAACGGTGGCTCCGTCACCGAGATTAATGGTGATTACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                                                                           <  1:3192431/147‑1 (MQ=17)
          aTACCGAAGGCGAGAGCACCATCACCACCGTTGGCACCGGCACCCAGATTCATGGTGATGACGCCACGGCAAACAACTACGGCAAAACCACCGTTGACTGCTAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                                                                           <  2:4093692/150‑1 (MQ=17)
          aTAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGCTTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                                                                           <  2:4093764/150‑1 (MQ=255)
           tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAGGGCAAACAACAACGGCTAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAGGttgtgca                                                                                                                                          >  2:3132188/1‑145 (MQ=14)
            aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                                                                                                                         >  2:3499593/1‑150 (MQ=255)
            aaCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                                                                                                                         <  2:3160079/150‑1 (MQ=255)
             aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                                                                                                                                        <  1:2927153/150‑1 (MQ=255)
             aCGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCTATCAACAAAGGCAAAACCACAGAAGACGGCAATGATGCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCTgg                                                                                                                                        >  1:2350879/1‑150 (MQ=17)
                aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACg                                                                                                                                     <  2:124332/150‑1 (MQ=255)
                      agagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                                                                                                               <  1:2771508/150‑1 (MQ=255)
                       gagCTCCATCACCAACGGTGGCACCGTCAACCATATTAATGGTGATTACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                                                                                                              <  1:1084508/150‑1 (MQ=17)
                       gagCACCATCACCACCGGTGGCACCGGCACACAGATTAATGGTGATGACGCCACGTCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                                                                                                              <  2:639709/150‑1 (MQ=255)
                        agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                                                                                                                             <  2:2430304/150‑1 (MQ=255)
                        agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                                                                                                                             <  2:2249792/150‑1 (MQ=255)
                        acctgcATCAACAAAGCTGCCTCCGTCACGGATATTACTGGTGATTACGCCACTGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                                                                                                                             <  2:58518/145‑1 (MQ=14)
                            ccAACACCAACGGTGGCTGCGCCACCCAGATTAATGGTGATGACGCCACTGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAGTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                                                                         <  1:3809054/150‑1 (MQ=17)
                                 aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                    >  2:1501075/1‑150 (MQ=255)
                                         tGGCACCGGCACCCAGATTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                            <  1:3544823/150‑1 (MQ=255)
                                          ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                                                                                                           >  2:3512122/1‑150 (MQ=255)
                                          ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                                                                                                           >  2:3512392/1‑150 (MQ=255)
                                           gcaccggcaccgATATTACGGGTGATGACGCCACTGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                                                                                                          <  1:2678909/150‑1 (MQ=17)
                                           gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                                                                                                          >  1:2545631/1‑150 (MQ=255)
                                           gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                                                                                                          >  1:462893/1‑150 (MQ=255)
                                               cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                                                                                                      >  2:361888/1‑150 (MQ=255)
                                               cggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                                                                                                      >  2:362256/1‑150 (MQ=255)
                                                    cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAATGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                 <  1:361888/150‑1 (MQ=255)
                                                     ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa                                                                                                >  2:2296986/1‑150 (MQ=255)
                                                                    tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAg                                                                                 <  1:3529322/150‑1 (MQ=255)
                                                                               cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTgg                                                                      <  2:3206996/150‑1 (MQ=255)
                                                                                 aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                    >  1:178034/1‑150 (MQ=255)
                                                                                 aacaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                    >  1:3430401/1‑150 (MQ=255)
                                                                                            cAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAAAGGTAATAAAGGTAAAGTTATTCAGGAAGGCGATCTGTATGTCAGCGGCGGCGGTCACGGTAGTGATATCACCGGCGACAGCGCAACGGTGGATAGCTAGGGcgc                                                         >  1:3947955/1‑148 (MQ=14)
                                                                                              aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCa                                                       >  2:3898058/1‑150 (MQ=255)
                                                                                                 ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                    >  2:2221400/1‑150 (MQ=255)
                                                                                                 ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                    >  1:946818/1‑150 (MQ=255)
                                                                                                  caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAc                                                   >  1:2400710/1‑150 (MQ=21)
                                                                                                   accGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAcc                                                  <  1:2296986/150‑1 (MQ=255)
                                                                                                   accGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAcc                                                  >  1:1971125/1‑150 (MQ=255)
                                                                                                    ccGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaa                                                               >  1:141777/1‑136 (MQ=255)
                                                                                                      gTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                               <  2:167675/150‑1 (MQ=255)
                                                                                                      gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAATGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATAGATATCACCGGCGACAGCGCAACGGTGGATAACACGGGCACCGGGAGCGTc                                               >  2:3630378/1‑150 (MQ=14)
                                                                                                      gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                      >  1:2259533/1‑143 (MQ=255)
                                                                                                      gTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                               <  2:946818/150‑1 (MQ=255)
                                                                                                        tGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcg                                                                               >  1:4048605/1‑116 (MQ=255)
                                                                                                            ggCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                         >  2:2589002/1‑150 (MQ=255)
                                                                                                            ggCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                         >  2:2588490/1‑150 (MQ=255)
                                                                                                            ggCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                         >  2:1895457/1‑150 (MQ=255)
                                                                                                              cAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                       >  1:3306197/1‑150 (MQ=255)
                                                                                                              cAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                       <  2:3247395/150‑1 (MQ=255)
                                                                                                               aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                      <  1:1594290/150‑1 (MQ=255)
                                                                                                               aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                      <  1:1594242/150‑1 (MQ=255)
                                                                                                               aaGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                      >  1:2831528/1‑150 (MQ=255)
                                                                                                                  gATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                                                                                                            <  2:659567/45‑1 (MQ=255)
                                                                                                                  gATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                   <  2:3306197/150‑1 (MQ=39)
                                                                                                                   aTTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                  >  1:3074295/1‑150 (MQ=39)
                                                                                                                     tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                                                                                                                     <  1:3430399/65‑1 (MQ=40)
                                                                                                                     tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATCACAAGGGCACCATGACCGTCACCGATCCGGAGTcc                                >  2:4106602/1‑150 (MQ=255)
                                                                                                                     tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTcc                                >  1:2246520/1‑150 (MQ=39)
                                                                                                                        aCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATc                             <  2:880134/150‑1 (MQ=11)
                                                                                                                         ccGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCg                            >  1:1680108/1‑150 (MQ=40)
                                                                                                                           ggTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                          >  2:3537220/1‑150 (MQ=255)
                                                                                                                            gTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                         <  1:831619/150‑1 (MQ=255)
                                                                                                                            gTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                         <  1:831624/150‑1 (MQ=21)
                                                                                                                            gTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                         <  1:1499432/150‑1 (MQ=17)
                                                                                                                            gTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                         <  1:1420546/150‑1 (MQ=255)
                                                                                                                            gTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGAGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                         >  1:689904/1‑150 (MQ=255)
                                                                                                                             tACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                        <  2:2400736/150‑1 (MQ=21)
                                                                                                                             tACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                        >  1:2419569/1‑150 (MQ=17)
                                                                                                                             tACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                        <  2:2400710/150‑1 (MQ=255)
                                                                                                                               cGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACTGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCTGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                      >  2:610440/1‑150 (MQ=255)
                                                                                                                                 gAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAg                    >  1:1129377/1‑150 (MQ=255)
                                                                                                                                  aaaTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAgg                   >  2:596224/1‑150 (MQ=17)
                                                                                                                                  aaaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCg                                                                                     >  2:3905918/1‑84 (MQ=40)
                                                                                                                                  aaaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAgg                   >  2:1049803/1‑150 (MQ=255)
                                                                                                                                   aaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                                                                                                              <  2:1394887/42‑1 (MQ=255)
                                                                                                                                   aaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATAGGTATCACCGGCGACAGCGCAACGGTGTATATCAAGGGCTCAATGACCGTCACCGATCCGGAGACCATAGGTATCCAGGt                  >  2:1764644/1‑150 (MQ=14)
                                                                                                                                    aTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt                                                                                                     <  1:3845631/66‑1 (MQ=40)
                                                                                                                                     tCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                    <  1:1051014/98‑1 (MQ=255)
                                                                                                                                      cAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGa               <  1:3437423/150‑1 (MQ=255)
                                                                                                                                       aaCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGAc              <  1:3364039/150‑1 (MQ=255)
                                                                                                                                         cGGTAATAACGGCAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACgg            <  1:2673094/150‑1 (MQ=255)
                                                                                                                                         cGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAAGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACgg            >  1:1809810/1‑150 (MQ=255)
                                                                                                                                          ggTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGc           >  2:2042618/1‑150 (MQ=12)
                                                                                                                                          ggTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACa                                                                >  1:1740533/1‑97 (MQ=255)
                                                                                                                                            taataaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCAGGAGTCCAGCGGTATCCAGGTTGACGGCGa         >  1:2112056/1‑150 (MQ=255)
                                                                                                                                                aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                       <  1:1953108/52‑1 (MQ=255)
                                                                                                                                                aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                       >  2:1953108/1‑52 (MQ=255)
                                                                                                                                                aaCGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAg     >  1:2858525/1‑150 (MQ=39)
                                                                                                                                                   ggAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                                                                                                          >  2:282955/1‑46 (MQ=255)
                                                                                                                                                   ggAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg  >  2:282931/1‑150 (MQ=40)
                                                                                                                                                    gAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                         >  2:275826/1‑110 (MQ=38)
                                                                                                                                                    gAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                         <  1:275826/110‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                                   
GCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCG  >  CP009273/1460540‑1460836

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: