Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,858 G→T 72.4% V71V (GTC→GTA)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,8580GT72.4% 365.6 / 252.4 283V71V (GTC→GTA) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base T (106/99);  total (140/143)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 285 total aligned reads displayed.

GCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGG  >  CP009273/4219711‑4220007
                                                                                                                                                   |                                                                                                                                                     
gcTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTAc                                                                                                                                                     <  1:2344950/150‑1 (MQ=255)
 cTCCAGCCACTATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACa                                                                                                                                                    <  2:3024192/150‑1 (MQ=255)
   ccAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAAc                                                                                                                                                  <  1:107221/150‑1 (MQ=255)
         aaTATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATg                                                                                                                                            <  2:4972376/150‑1 (MQ=255)
          atatACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGa                                                                                                                                           <  2:3059790/150‑1 (MQ=255)
           tataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGaa                                                                                                                                          >  1:77228/1‑150 (MQ=255)
           tataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGaa                                                                                                                                          >  1:3549791/1‑150 (MQ=255)
           tataCAGAGCGCCACGTTTAACCAAGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAGCGCGGCACCGGGATCTACAACGCTATGaa                                                                                                                                          >  2:4945672/1‑150 (MQ=255)
            ataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAAt                                                                                                                                         <  1:4928857/150‑1 (MQ=255)
             taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATa                                                                                                                                        <  2:2447421/150‑1 (MQ=255)
              aCAGACCGCCCCGTTTCACCACGTGCGTAAGAGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTCCTGCAACTGGAAAGTATTCCCGCCGCCTACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATAc                                                                                                                                       <  2:2676238/150‑1 (MQ=255)
                agagCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCg                                                                                                                                     <  1:2703451/150‑1 (MQ=255)
                 gagCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCgg                                                                                                                                    <  2:2820669/150‑1 (MQ=255)
                  agCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGt                                                                                                                                   <  1:1419878/150‑1 (MQ=255)
                   gcgcCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACAGGTg                                                                                                                                  >  2:1051159/1‑150 (MQ=255)
                     gcCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGAc                                                                                                                                >  1:598174/1‑150 (MQ=255)
                      ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACa                                                                                                                               <  1:3760549/150‑1 (MQ=255)
                      ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACa                                                                                                                               >  1:2521251/1‑150 (MQ=255)
                         cGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGaa                                                                                                                            <  2:189060/150‑1 (MQ=255)
                                 cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAg                                                                                                                    >  1:3227081/1‑150 (MQ=255)
                                   cGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCg                                                                                                                  <  2:743884/150‑1 (MQ=255)
                                   cGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCg                                                                                                                  >  2:3569508/1‑150 (MQ=255)
                                     tCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGa                                                                                                                <  1:5131411/150‑1 (MQ=255)
                                      cAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAg                                                                                                               <  2:4649772/150‑1 (MQ=255)
                                      cAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAg                                                                                                               >  1:4655530/1‑150 (MQ=255)
                                       aGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGc                                                                                                              <  2:4414059/150‑1 (MQ=255)
                                          aaTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgag                                                                                                           >  1:1953101/1‑150 (MQ=255)
                                          aaTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgag                                                                                                           >  2:292296/1‑150 (MQ=255)
                                           aTTGTTGACAGCAGAACGCGCTAGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGACGACAATCACGATCTCAGCATTTTCAATCGCTGCAACGGGATCTACAACGATATGAATAACGGAGACAGAAACACCCAGCGGAGCgcga                                                                                                          >  2:1657258/1‑150 (MQ=255)
                                                tGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGc                                                                                                     >  1:86941/1‑150 (MQ=255)
                                                 gCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCa                                                                                                    >  2:2756381/1‑150 (MQ=255)
                                                  ccAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTCTGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAg                                                                                                   >  2:3935567/1‑150 (MQ=255)
                                                   cagcagCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGc                                                                                                  <  2:2266506/150‑1 (MQ=255)
                                                    agcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                                                                                                 <  1:2054603/150‑1 (MQ=255)
                                                    agcagcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                                                                                                 >  2:5026703/1‑150 (MQ=255)
                                                    agcagcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                                                                                                 <  2:930978/150‑1 (MQ=255)
                                                     gcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCg                                                                                                <  1:1140940/150‑1 (MQ=255)
                                                      cagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGt                                                                                               >  1:4736556/1‑150 (MQ=255)
                                                       agcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                                                                                              <  2:1488047/150‑1 (MQ=255)
                                                       agcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                                                                                              >  1:3662503/1‑150 (MQ=255)
                                                        gcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGttt                                                                                             <  2:3127005/150‑1 (MQ=255)
                                                        gcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGttt                                                                                             >  2:1460220/1‑150 (MQ=255)
                                                             cgcTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCt                                                                                        <  1:3384461/150‑1 (MQ=255)
                                                              gcTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtg                                                                                       >  1:4545724/1‑150 (MQ=255)
                                                                 cgcgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTgtgt                                                                                    <  2:1567443/150‑1 (MQ=255)
                                                                   cgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTaa                                                                                  <  2:3786001/150‑1 (MQ=255)
                                                                    ggCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAat                                                                                 >  2:1324400/1‑150 (MQ=255)
                                                                    ggCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAat                                                                                 >  1:4287476/1‑150 (MQ=255)
                                                                      cACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatca                                                                               <  1:1281735/150‑1 (MQ=255)
                                                                               tgCAACTGGAAAGTATTCCCGCCTCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGACCCGCAGCCGTTTTCGCTGTGTAATCATACCAGGTc                                                                      >  1:2321223/2‑150 (MQ=255)
                                                                                gCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCt                                                                     <  1:983888/150‑1 (MQ=255)
                                                                                   aCTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCg                                                                  >  1:521308/1‑150 (MQ=255)
                                                                                      ggAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTa                                                               <  1:408643/150‑1 (MQ=255)
                                                                                       gAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTAGCTGTGTAATCATCCCAGGTCTGCGTTAc                                                              >  1:2254188/1‑150 (MQ=255)
                                                                                          aGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGcc                                                           >  2:64210/1‑150 (MQ=255)
                                                                                               tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGa                                                      >  1:5132749/1‑150 (MQ=255)
                                                                                               tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGa                                                      >  1:5156530/1‑150 (MQ=255)
                                                                                                cccGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaa                                                     >  2:2738905/1‑150 (MQ=255)
                                                                                                  cgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAg                                                   >  2:3258226/1‑150 (MQ=255)
                                                                                                  cgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCAAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCAACCCAGGTCTGCGTTACGCCAGCGAAAg                                                   >  2:4850242/1‑150 (MQ=255)
                                                                                                  cgccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAg                                                   >  1:3130169/1‑150 (MQ=255)
                                                                                                   gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAgg                                                  >  2:1069280/1‑150 (MQ=255)
                                                                                                    ccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAggg                                                 <  2:3742627/150‑1 (MQ=255)
                                                                                                    ccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAggg                                                 >  2:4709189/1‑150 (MQ=255)
                                                                                                     cgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGa                                                >  1:2160561/1‑150 (MQ=255)
                                                                                                     cgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGa                                                >  2:3281747/1‑150 (MQ=255)
                                                                                                      gccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGAt                                               >  2:3855100/1‑150 (MQ=255)
                                                                                                        cgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaa                                             <  2:2128021/150‑1 (MQ=255)
                                                                                                         gACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaaa                                            >  1:2633257/1‑150 (MQ=255)
                                                                                                          aCAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAac                                           <  1:5099676/150‑1 (MQ=255)
                                                                                                           cAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCTGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAaca                                          <  1:2114949/150‑1 (MQ=255)
                                                                                                            aaTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAacac                                         >  2:1300698/1‑150 (MQ=255)
                                                                                                             aTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCTGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACAcc                                        <  2:1813146/150‑1 (MQ=255)
                                                                                                               cACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGc                                      >  1:2991203/1‑150 (MQ=255)
                                                                                                                 cGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTg                                    >  1:1825673/1‑150 (MQ=255)
                                                                                                                 cGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTg                                    <  1:1894302/150‑1 (MQ=255)
                                                                                                                  gATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGGGTTACGCCAGCGAAAGGGATAAACACCGCTGa                                   >  1:2866611/1‑150 (MQ=255)
                                                                                                                   aTCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAg                                  <  2:2340000/150‑1 (MQ=255)
                                                                                                                   aTCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAg                                  >  2:4066200/1‑150 (MQ=255)
                                                                                                                   aTCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAg                                  >  2:4085356/1‑150 (MQ=255)
                                                                                                                    tctcAGCATTTTCAATCGCGGCAACGGGATCTACAAAGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAAACGCAGCCGTATTCGCAGTGTAATCAAACCAGGTCATCGTTACGCCAGCGAAATGGATAAACACCCCTGAgc                                 >  1:4100619/1‑150 (MQ=255)
                                                                                                                     ctcAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcg                                >  1:3620973/1‑150 (MQ=255)
                                                                                                                     ctcAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcg                                >  2:4117268/1‑150 (MQ=255)
                                                                                                                       cAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggc                              >  1:4291829/1‑150 (MQ=255)
                                                                                                                        aGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggcg                             >  1:4456929/1‑150 (MQ=255)
                                                                                                                        aGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggcg                             <  1:3588039/150‑1 (MQ=255)
                                                                                                                                cAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCa                     >  2:5112249/1‑150 (MQ=255)
                                                                                                                                cAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCa                     >  1:1733349/1‑150 (MQ=255)
                                                                                                                                 aaTCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                    <  1:719680/150‑1 (MQ=255)
                                                                                                                                    cgcgGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTg                 >  2:4280481/1‑150 (MQ=255)
                                                                                                                                    cgcgGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTg                 >  1:2821169/1‑150 (MQ=255)
                                                                                                                                    cgcgGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTg                 >  1:4228119/1‑150 (MQ=255)
                                                                                                                                      cgGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGtt               >  1:4867028/1‑150 (MQ=255)
                                                                                                                                      cgGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGtt               >  1:669475/1‑150 (MQ=255)
                                                                                                                                         cAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAg            >  2:2904052/1‑150 (MQ=255)
                                                                                                                                         cAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAg            <  1:1578026/150‑1 (MQ=255)
                                                                                                                                              ggATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt       <  1:2756381/150‑1 (MQ=255)
                                                                                                                                              ggATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt       >  2:2853195/1‑150 (MQ=255)
                                                                                                                                               gATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTa      <  1:398673/150‑1 (MQ=255)
                                                                                                                                                aTCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAg     <  1:2450952/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGG  >  CP009273/4219711‑4220007

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: ‑