Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,041 A→T 100% pseudogene (1393/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,461,0410AT95.6% 874.8 / 21.0 276pseudogene (1393/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base T (177/86);  total (184/92)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 276 total aligned reads displayed.

AACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  CP009273/1460894‑1461187
                                                                                                                                                   |                                                                                                                                                  
aaCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCACTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTcc                                                                                                                                                  >  1:2962109/1‑150 (MQ=255)
aaCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTcc                                                                                                                                                  >  2:189342/1‑150 (MQ=255)
 aCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCg                                                                                                                                                 <  2:3864141/150‑1 (MQ=255)
  ctgTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCgg                                                                                                                                                >  1:922467/3‑150 (MQ=255)
    gTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCg                                                                                                                                              <  2:3522497/150‑1 (MQ=255)
        tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAg                                                                                                                                          >  2:3956042/1‑150 (MQ=255)
         gaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgc                                                                                                                                         <  1:864299/150‑1 (MQ=255)
           cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcgc                                                                                                                                       >  2:1383103/1‑150 (MQ=255)
                 ggCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGt                                                                                                                                 <  2:4967448/150‑1 (MQ=255)
                  gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTg                                                                                                                                >  1:1951353/1‑150 (MQ=255)
                    aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                              >  1:5137754/1‑150 (MQ=255)
                    aaaCAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                              >  1:3275892/1‑150 (MQ=255)
                      acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGTAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATa                                                                                                                            <  1:2648779/150‑1 (MQ=255)
                       caacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATaa                                                                                                                           >  1:1873136/1‑150 (MQ=255)
                          caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                        >  1:4548079/1‑150 (MQ=255)
                          caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                        >  1:2600235/1‑150 (MQ=255)
                           aacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACCGCGCAACGGTGGATAACAAg                                                                                                                       >  2:4560848/1‑150 (MQ=255)
                                  aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCa                                                                                                                <  1:4798744/150‑1 (MQ=255)
                                   aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                               >  2:1614810/1‑150 (MQ=255)
                                   aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                               >  1:4634476/1‑150 (MQ=255)
                                   aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCCCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                               <  2:582459/150‑1 (MQ=255)
                                    accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATg                                                                                                              <  1:3089230/150‑1 (MQ=255)
                                     ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                             >  1:1505609/1‑150 (MQ=255)
                                     ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                             >  1:1342156/1‑150 (MQ=255)
                                     ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                             >  1:3823754/1‑150 (MQ=255)
                                     ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                             >  1:3986236/1‑150 (MQ=255)
                                      caccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAc                                                                                                            >  1:929659/1‑150 (MQ=255)
                                        ccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCg                                                                                                          >  2:2380833/1‑150 (MQ=255)
                                        ccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCg                                                                                                          <  2:2788275/150‑1 (MQ=255)
                                          gTTGACGGCAAGGATTCCACCGGTACTGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAAAAAGGGCACCATGACCGTc                                                                                                        <  2:1109410/150‑1 (MQ=255)
                                          gTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                                                                                        >  1:3152955/1‑150 (MQ=255)
                                            tGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                      >  2:1048219/1‑150 (MQ=255)
                                               cGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGa                                                                                                   >  2:739748/1‑150 (MQ=255)
                                               cGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGa                                                                                                   >  1:359969/1‑150 (MQ=255)
                                                   aaGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                               >  1:3339950/1‑150 (MQ=255)
                                                     ggATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                             >  2:1535070/1‑150 (MQ=255)
                                                     ggATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                             <  2:3287192/150‑1 (MQ=255)
                                                      gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                                                                            <  1:1047127/150‑1 (MQ=255)
                                                         tCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCTACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTcc                                                                                         >  2:1331745/1‑150 (MQ=255)
                                                             ccGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCg                                                                                     >  1:4681802/1‑150 (MQ=255)
                                                               ggTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                                                                                   >  2:2948512/1‑150 (MQ=255)
                                                                  aCGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATc                                                                                >  2:2063223/1‑150 (MQ=255)
                                                                   cGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                                                                               <  2:3503353/150‑1 (MQ=255)
                                                                   cGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc                                                                               <  1:4326084/150‑1 (MQ=255)
                                                                    ggAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCa                                                                              <  2:4414101/150‑1 (MQ=255)
                                                                     gAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAg                                                                             >  1:4786595/1‑150 (MQ=255)
                                                                       aaTTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGAt                                                                           >  2:279106/1‑150 (MQ=255)
                                                                           gCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGAc                                                                       >  2:4593177/1‑150 (MQ=255)
                                                                            cTGTCAATATCGGGAAGGTGATTCCTGACGGCGATCTGGATGTCATCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCATCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACg                                                                      <  2:4828311/150‑1 (MQ=14)
                                                                            cTGGCAATAACGGGAAGGTGATTCAGCACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACg                                                                      >  1:3524433/1‑150 (MQ=255)
                                                                             tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                     >  1:3751268/1‑150 (MQ=255)
                                                                              ggCAAGAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGc                                                                    <  1:131575/150‑1 (MQ=255)
                                                                               gtaaTAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGa                                                                                                                              <  2:981394/89‑1 (MQ=9)
                                                                                cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGa                                                                  >  1:492182/1‑150 (MQ=255)
                                                                                cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGa                                                                  >  1:1831922/1‑150 (MQ=255)
                                                                                cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGa                                                                  >  1:4395175/1‑150 (MQ=255)
                                                                                 aaTAACGGTAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAc                                                                 <  1:212355/150‑1 (MQ=12)
                                                                                 aaTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAc                                                                 >  1:1000389/1‑150 (MQ=255)
                                                                                 aaTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAc                                                                 >  2:3121318/1‑150 (MQ=255)
                                                                                  aTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAcc                                                                >  1:2309839/1‑150 (MQ=255)
                                                                                      cGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc                                                            >  2:1781857/1‑150 (MQ=255)
                                                                                      cGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc                                                            >  1:2152154/1‑150 (MQ=255)
                                                                                           aagTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAgg                                                                                                                      >  1:4083263/3‑87 (MQ=14)
                                                                                           aagTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAgg                                                                                                                      <  2:4083263/85‑1 (MQ=14)
                                                                                           aGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCg                                                       >  2:1339180/1‑150 (MQ=255)
                                                                                            ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGt                                                      <  1:1571505/150‑1 (MQ=255)
                                                                                              tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCa                                                    >  1:1193076/1‑150 (MQ=255)
                                                                                              tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                                     <  2:2263152/101‑1 (MQ=31)
                                                                                                aTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                  >  1:4200664/1‑150 (MQ=255)
                                                                                                 ttCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATa                                                 <  2:2197431/150‑1 (MQ=255)
                                                                                                  tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                >  1:2835121/1‑150 (MQ=255)
                                                                                                    tggACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                              >  1:1416732/2‑150 (MQ=255)
                                                                                                    aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                              >  2:2307904/1‑150 (MQ=255)
                                                                                                    aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                              >  2:538350/1‑150 (MQ=255)
                                                                                                     ggACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGa                                             >  2:214063/1‑150 (MQ=255)
                                                                                                      gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGaa                                            >  2:4217297/1‑150 (MQ=255)
                                                                                                        cGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAgg                                          <  1:615228/150‑1 (MQ=255)
                                                                                                         ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGc                                         >  2:3784383/1‑150 (MQ=255)
                                                                                                         ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGc                                         >  1:5014170/1‑150 (MQ=255)
                                                                                                           cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGTACGTCACCGATCCGGAGTCCATCTGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                       <  1:4947817/150‑1 (MQ=255)
                                                                                                            gATCTGGATGTCAGCGGCGGCGGTCACGGTAGTGATATCTCCGGCGACAGCGCAACGGTGGATAACACGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgag                                      <  2:697597/150‑1 (MQ=255)
                                                                                                             aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgaga                                     >  2:2969395/1‑150 (MQ=255)
                                                                                                                tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                  >  2:4584709/1‑150 (MQ=255)
                                                                                                                tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                  >  2:4527273/1‑150 (MQ=255)
                                                                                                                    tGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                              >  1:2416508/1‑150 (MQ=255)
                                                                                                                     gTCAGCGGCGGTGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATc                             >  1:1632414/1‑150 (MQ=255)
                                                                                                                      tCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                            >  2:3408399/1‑150 (MQ=255)
                                                                                                                        aGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                          >  1:443182/1‑150 (MQ=255)
                                                                                                                          cggcggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaa                        >  2:2028214/1‑150 (MQ=255)
                                                                                                                          cggcggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaa                        <  2:2562404/150‑1 (MQ=255)
                                                                                                                           ggcggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAt                       <  2:2631417/150‑1 (MQ=255)
                                                                                                                             cggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGGATCCAGATTGACGGCGACCAGGCCATCGGCAATACCGACGGCGAGAGCACTATCTCCACTgg                     >  2:3388904/1‑150 (MQ=14)
                                                                                                                             cggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATgg                     <  1:722292/150‑1 (MQ=255)
                                                                                                                               gcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcg                   <  2:1861856/150‑1 (MQ=255)
                                                                                                                                cggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcgg                  >  1:1131120/1‑150 (MQ=255)
                                                                                                                                  gTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                >  2:1500098/1‑150 (MQ=255)
                                                                                                                                   tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac               >  2:3995532/1‑150 (MQ=255)
                                                                                                                                     aCGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg             <  2:363139/150‑1 (MQ=255)
                                                                                                                                     aCGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccg             <  1:4714425/150‑1 (MQ=255)
                                                                                                                                        gTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggca          <  2:323272/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  CP009273/1460894‑1461187

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: