Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,894,614 C→A 15.5% A323S (GCA→TCA)  bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,894,6140CA15.5% 956.2 / 110.5 305A323S (GCA→TCA) bglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (116/141);  new base A (25/23);  total (141/164)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 313 total aligned reads displayed.

AAAGTTACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAA  >  CP009273/3894465‑3894763
                                                                                                                                                     |                                                                                                                                                     
aaaGTTACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATgc                                                                                                                                                       <  2:925210/150‑1 (MQ=255)
aaaGTTACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATgc                                                                                                                                                       <  2:3427181/150‑1 (MQ=255)
aaaGTGTCTGGCGATAGAGTTATTCGCGACCAGGAAGGAGTATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATgc                                                                                                                                                       <  2:641606/150‑1 (MQ=255)
  aGTTACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATgcgc                                                                                                                                                     >  1:476896/1‑150 (MQ=255)
   gTTACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCt                                                                                                                                                    >  1:3778476/1‑150 (MQ=255)
     tACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCttt                                                                                                                                                  <  1:5120793/150‑1 (MQ=255)
       cTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCg                                                                                                                                                <  1:4658080/150‑1 (MQ=255)
       cTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCg                                                                                                                                                <  1:2035452/150‑1 (MQ=255)
            gATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGc                                                                                                                                           >  1:3534953/1‑150 (MQ=255)
            gATAGAGTTAGTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGc                                                                                                                                           <  2:4035099/150‑1 (MQ=255)
               agagTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGt                                                                                                                                        >  2:280755/1‑150 (MQ=255)
                     aTTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAACTTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATa                                                                                                                                  <  2:3878491/150‑1 (MQ=255)
                      ttCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATaa                                                                                                                                 <  1:3252053/150‑1 (MQ=255)
                        cgcgACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACg                                                                                                                               <  2:1281971/150‑1 (MQ=255)
                          cgACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGAc                                                                                                                             <  1:987687/150‑1 (MQ=255)
                             ccAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCAc                                                                                                                          <  1:2561636/150‑1 (MQ=255)
                               aggaaggaGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACtt                                                                                                                        <  2:529236/150‑1 (MQ=255)
                               aggaaggaGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACtt                                                                                                                        <  2:2330127/150‑1 (MQ=255)
                                ggaaggaGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTa                                                                                                                       <  1:4358790/150‑1 (MQ=255)
                                ggaaggaGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTa                                                                                                                       <  1:1803312/150‑1 (MQ=255)
                                   aggaGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCa                                                                                                                    >  2:3321465/1‑150 (MQ=255)
                                        aaTTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAg                                                                                                               <  2:1204872/150‑1 (MQ=255)
                                        aaTTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAg                                                                                                               <  2:4467478/150‑1 (MQ=255)
                                         agtCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAgg                                                                                                              >  1:2012614/3‑150 (MQ=255)
                                         aTTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAgg                                                                                                              >  1:2700885/1‑150 (MQ=255)
                                          ttCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGt                                                                                                             >  2:4531573/1‑150 (MQ=255)
                                          ttCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGt                                                                                                             >  2:2960696/1‑150 (MQ=255)
                                          ttCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGt                                                                                                             >  1:558275/1‑150 (MQ=255)
                                             gTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTACCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGc                                                                                                          >  2:1707150/1‑150 (MQ=255)
                                               tGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCtt                                                                                                        >  2:5085913/1‑150 (MQ=255)
                                                  aGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTAt                                                                                                     >  1:1404863/1‑150 (MQ=255)
                                                   tcaaacTTTATCAAATTTCTGAGTTAAAGATGTGACAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCGTGATTAACTTTTTCAGATGAGCTTTCTTTTTCAGTGACATAACTACCACTAACCATTAATGTGGCTTTATc                                                                                                    >  2:5010725/6‑150 (MQ=255)
                                                     cACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATccc                                                                                                  >  2:3239050/1‑150 (MQ=255)
                                                      aCCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCa                                                                                                 >  1:633056/1‑150 (MQ=255)
                                                       ccTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCAt                                                                                                >  2:3376299/1‑150 (MQ=255)
                                                       ccTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCAt                                                                                                >  2:2436079/1‑150 (MQ=255)
                                                        cTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATa                                                                                               >  1:2891276/1‑150 (MQ=255)
                                                          ttATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAAc                                                                                             >  2:3959119/1‑150 (MQ=255)
                                                          ttATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAAc                                                                                             <  2:74367/150‑1 (MQ=255)
                                                           tATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACg                                                                                            <  1:1150245/150‑1 (MQ=255)
                                                            aTCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACgg                                                                                           <  2:2172817/150‑1 (MQ=255)
                                                             tCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACggg                                                                                          <  1:1610903/150‑1 (MQ=255)
                                                              cAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGa                                                                                         <  1:3216363/150‑1 (MQ=255)
                                                                 aTTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGatat                                                                                      >  2:1894605/1‑150 (MQ=255)
                                                                   ttCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCt                                                                                    <  1:1021650/150‑1 (MQ=255)
                                                                    tCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCtt                                                                                   <  2:3535255/150‑1 (MQ=255)
                                                                    tCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCtt                                                                                   <  1:3122761/150‑1 (MQ=255)
                                                                         gTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAg                                                                              <  2:2584078/150‑1 (MQ=255)
                                                                         gTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAg                                                                              <  1:3718356/150‑1 (MQ=255)
                                                                          ttAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGTCTTTATCCCATAACGGGATATCTTTATAGc                                                                             <  2:4310084/150‑1 (MQ=255)
                                                                             aaGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTa                                                                          >  2:3083103/1‑150 (MQ=255)
                                                                              aGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTaa                                                                         >  2:3658294/1‑150 (MQ=255)
                                                                              aGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTaa                                                                         >  2:4678606/1‑150 (MQ=255)
                                                                                  gtgCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCa                                                                     >  2:1432179/1‑150 (MQ=255)
                                                                                  gtgCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCa                                                                     <  2:2253212/150‑1 (MQ=255)
                                                                                    gCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCTTAAATCAAt                                                                   <  1:3848123/150‑1 (MQ=255)
                                                                                    gCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAt                                                                   <  1:3053329/150‑1 (MQ=255)
                                                                                    gCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAt                                                                   <  2:1904434/150‑1 (MQ=255)
                                                                                       aaaCATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAtggt                                                                >  1:1937304/1‑150 (MQ=255)
                                                                                       aaaCATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAtggt                                                                <  2:2243845/150‑1 (MQ=255)
                                                                                       aaaCATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAtggt                                                                <  1:1154095/150‑1 (MQ=255)
                                                                                           aTCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAtggtattg                                                            <  2:1459544/150‑1 (MQ=255)
                                                                                             ccAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAAtggtattggt                                                          <  1:1504136/150‑1 (MQ=255)
                                                                                                ggTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATggtattggtatt                                                       <  2:1832650/150‑1 (MQ=255)
                                                                                                  tATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTaa                                                     <  2:5127588/150‑1 (MQ=255)
                                                                                                   aTCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAAt                                                    >  1:1947374/1‑150 (MQ=255)
                                                                                                    tCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATc                                                   <  1:2670944/150‑1 (MQ=255)
                                                                                                    tCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATc                                                   >  2:154047/1‑150 (MQ=255)
                                                                                                     cTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATct                                                  <  1:3267511/150‑1 (MQ=255)
                                                                                                      tttCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATctg                                                 >  2:4365817/1‑150 (MQ=255)
                                                                                                           aGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCtgtt                                            <  2:2369947/150‑1 (MQ=255)
                                                                                                            gTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCtgttt                                           <  2:2029133/150‑1 (MQ=255)
                                                                                                             tCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCtgtttg                                          >  2:3122283/1‑150 (MQ=255)
                                                                                                              cATAATACACTTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTGTGCCGTGACATATCGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCtgtttgt                                         <  1:3957458/150‑1 (MQ=255)
                                                                                                               ataagaCCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCtgtttgtt                                        <  1:1221417/150‑1 (MQ=255)
                                                                                                                 aataCCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGtttt                                      >  2:1773211/1‑150 (MQ=255)
                                                                                                                 aataCCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGtttt                                      >  1:1862986/1‑150 (MQ=255)
                                                                                                                 aataCCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGtttt                                      <  2:543344/150‑1 (MQ=255)
                                                                                                                 aataCCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGtttt                                      <  1:2256347/150‑1 (MQ=255)
                                                                                                                  ataCCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTg                                     >  1:2648285/1‑150 (MQ=255)
                                                                                                                     cccGTTATTATCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAg                                  >  1:4180986/1‑150 (MQ=255)
                                                                                                                       cGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGc                                >  1:23350/1‑150 (MQ=255)
                                                                                                                          tattatCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGc                             <  1:3079523/150‑1 (MQ=255)
                                                                                                                          tattatCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGc                             >  1:691179/1‑150 (MQ=255)
                                                                                                                           attatCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCg                            <  1:1579335/150‑1 (MQ=255)
                                                                                                                           attatCCTGATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCg                            >  1:3322635/1‑150 (MQ=255)
                                                                                                                            ttatCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGa                           >  2:671405/1‑150 (MQ=255)
                                                                                                                            ttatCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGa                           <  1:4564626/150‑1 (MQ=255)
                                                                                                                              atCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATc                         <  2:3778516/150‑1 (MQ=255)
                                                                                                                              atCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATc                         >  2:4492412/1‑150 (MQ=255)
                                                                                                                                ccTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGt                       <  1:3698121/150‑1 (MQ=255)
                                                                                                                                  tGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTaa                     >  1:1778672/1‑150 (MQ=255)
                                                                                                                                   gATTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAat                    <  1:4467714/150‑1 (MQ=255)
                                                                                                                                    aTTATCTTTTTCCGATGAGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAatc                   <  1:2397867/150‑1 (MQ=255)
                                                                                                                                                cGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTc       <  1:4954966/150‑1 (MQ=255)
                                                                                                                                                cGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTc       >  1:4072952/1‑150 (MQ=255)
                                                                                                                                                   tGCGCTGTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTcgcg    <  2:3165417/150‑1 (MQ=255)
                                                                                                                                                    gcgcTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGa   <  2:1239152/150‑1 (MQ=255)
                                                                                                                                                    gcgcTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGa   <  1:4012821/150‑1 (MQ=255)
                                                                                                                                                    gcgcTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGa   >  1:1719093/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGTTACTGGCGATAGAGTTATTCGCGACCAGGAAGGAGAATTCGTTGAAGCCACCTTTATCAAATTTCTGCGTTAAAGATGTGCCAAACATCCAGGTATCTTTCCAGTCATAATACCCGTTATTATCCTGATTATCTTTTTCCGATGCGCTTTCGTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAA  >  CP009273/3894465‑3894763

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: