Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 305,655 G→T 6.5% S28* (TCG→TAG)  ecpB ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273305,6550GT 6.5% 1186.3 / 35.7 307S28* (TCG→TAG) ecpBECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (157/130);  new base T (9/11);  total (166/141)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 318 total aligned reads displayed.

GTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTC  >  CP009273/305506‑305804
                                                                                                                                                     |                                                                                                                                                     
gtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACg                                                                                                                                                       <  2:3736388/150‑1 (MQ=255)
  aGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGAGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACgat                                                                                                                                                     >  2:4481393/1‑150 (MQ=255)
      aCTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGata                                                                                                                                                 >  1:1029110/1‑150 (MQ=255)
        tCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGatata                                                                                                                                               >  2:4058605/1‑150 (MQ=255)
         cATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGatatat                                                                                                                                              >  1:3451838/1‑150 (MQ=255)
         cATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGatatat                                                                                                                                              <  2:3834753/150‑1 (MQ=255)
          aTCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATc                                                                                                                                             <  2:4785579/150‑1 (MQ=255)
             cGGCTTGTCCATTGAGATTACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGcc                                                                                                                                          <  1:4581585/150‑1 (MQ=255)
               gCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGa                                                                                                                                        >  2:2053203/1‑150 (MQ=255)
               gCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGa                                                                                                                                        >  1:132743/1‑150 (MQ=255)
                       aTTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACTATGATATATCGCCGACATCCAgc                                                                                                                                >  1:2425611/1‑150 (MQ=255)
                        ttGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAgcg                                                                                                                               <  1:3134353/150‑1 (MQ=255)
                          gagaTAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGcc                                                                                                                             >  1:604919/1‑150 (MQ=255)
                           agaTAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACTATGATATATCGCCGACATCCAGCGCCt                                                                                                                            <  2:1947810/150‑1 (MQ=255)
                            gaTAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTg                                                                                                                           >  1:2135820/1‑150 (MQ=255)
                              tAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGGTCATAAACGATGATATATCGCCGACATCCAGCGCCTggg                                                                                                                         >  2:5043076/1‑150 (MQ=255)
                                      ccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGat                                                                                                                 >  2:3203806/1‑150 (MQ=255)
                                      ccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGat                                                                                                                 <  1:2217432/150‑1 (MQ=255)
                                          cgtcAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTc                                                                                                             <  2:4547095/150‑1 (MQ=255)
                                           gtcAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTcc                                                                                                            >  1:1552792/1‑150 (MQ=255)
                                           gtcAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTcc                                                                                                            >  1:3549996/1‑150 (MQ=255)
                                             cAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCgg                                                                                                          <  2:3452096/150‑1 (MQ=255)
                                                   gTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAAca                                                                                                    <  1:270170/150‑1 (MQ=255)
                                                    tGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAAcag                                                                                                   >  2:1956476/1‑150 (MQ=255)
                                                    tGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAAcag                                                                                                   <  1:2092131/150‑1 (MQ=255)
                                                        gagagCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCaga                                                                                               <  1:5121175/150‑1 (MQ=255)
                                                        gagagCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCaga                                                                                               >  1:1270495/1‑150 (MQ=255)
                                                         agagCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCagag                                                                                              >  2:4330373/1‑150 (MQ=255)
                                                          gagCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGc                                                                                             <  2:3666189/150‑1 (MQ=255)
                                                               gTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGagag                                                                                        >  1:670262/1‑150 (MQ=255)
                                                                  cGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCa                                                                                     <  1:4874426/150‑1 (MQ=255)
                                                                     gaCGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGaa                                                                                  >  1:2511478/1‑150 (MQ=255)
                                                                     gaCGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGaa                                                                                  >  1:2898634/1‑150 (MQ=255)
                                                                       cGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAgg                                                                                <  2:1741547/150‑1 (MQ=255)
                                                                       cGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAgg                                                                                >  1:1181460/1‑150 (MQ=255)
                                                                       cGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAgg                                                                                >  1:237638/1‑150 (MQ=255)
                                                                         gATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTg                                                                              >  1:4184570/1‑150 (MQ=255)
                                                                         gATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTg                                                                              >  1:3908457/1‑150 (MQ=255)
                                                                         gATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGAGATTCCGGAAAACAGCAGAGCGCGAGGCAGGAGGTg                                                                              >  2:4807389/1‑150 (MQ=255)
                                                                          atatTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGc                                                                             >  1:1937984/1‑150 (MQ=255)
                                                                           tatTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCt                                                                            <  2:4455252/150‑1 (MQ=255)
                                                                            atTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCtt                                                                           <  1:1902335/150‑1 (MQ=255)
                                                                            atTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCtt                                                                           <  2:2536444/150‑1 (MQ=255)
                                                                              tGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCATTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCtttt                                                                         <  2:3948679/150‑1 (MQ=255)
                                                                                aTAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTc                                                                       <  1:4266240/150‑1 (MQ=255)
                                                                                  aaGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCAt                                                                     <  2:2751326/150‑1 (MQ=255)
                                                                                   aGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATa                                                                    >  2:3593660/1‑150 (MQ=255)
                                                                                     gCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAAc                                                                  >  2:3926728/1‑150 (MQ=255)
                                                                                     gCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAAc                                                                  <  1:2673394/150‑1 (MQ=255)
                                                                                      cGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACTATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACa                                                                 <  1:3783986/150‑1 (MQ=255)
                                                                                       gACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACaa                                                                >  1:3359718/1‑150 (MQ=255)
                                                                                        aCCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACTATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAAt                                                               <  1:139105/150‑1 (MQ=255)
                                                                                          cACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTc                                                             <  2:1494447/150‑1 (MQ=255)
                                                                                           aCTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTcc                                                            <  2:1297068/150‑1 (MQ=255)
                                                                                            cTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCg                                                           <  1:5004775/150‑1 (MQ=255)
                                                                                             tGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGt                                                          >  1:720460/1‑150 (MQ=255)
                                                                                              gTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGCGCTTTTTCATAACAATTCCGTc                                                         <  2:2846616/150‑1 (MQ=255)
                                                                                                cGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGACGGTGCTTTTTCATAACAATTCCGTCCa                                                       >  2:1643990/1‑150 (MQ=255)
                                                                                                      tGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATa                                                 >  1:2149858/1‑150 (MQ=255)
                                                                                                       gTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaa                                                <  2:1712374/150‑1 (MQ=255)
                                                                                                        tttttGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaaa                                               >  1:1474229/1‑150 (MQ=255)
                                                                                                        tttttGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaaa                                               <  2:1230782/150‑1 (MQ=255)
                                                                                                        tttttGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaaa                                               <  1:3065275/150‑1 (MQ=255)
                                                                                                            tGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAg                                           >  2:4830368/1‑150 (MQ=255)
                                                                                                             gATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgg                                          >  1:474158/1‑150 (MQ=255)
                                                                                                              aTCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggg                                         >  1:2251522/1‑150 (MQ=255)
                                                                                                               tCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggg                                        <  2:2180555/150‑1 (MQ=255)
                                                                                                               tCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggg                                        <  1:3591939/150‑1 (MQ=255)
                                                                                                               tCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggg                                        >  1:15098/1‑150 (MQ=255)
                                                                                                               tCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGAGGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAgggg                                        <  2:5123397/150‑1 (MQ=255)
                                                                                                                cGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATATAGCAggggg                                       <  2:1548868/150‑1 (MQ=255)
                                                                                                                cGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggggg                                       >  2:538797/1‑150 (MQ=255)
                                                                                                                 gTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGACTAAAGCAgggggg                                      >  2:4706373/1‑150 (MQ=255)
                                                                                                                   tttGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCt                                    <  2:2672410/150‑1 (MQ=255)
                                                                                                                   tttGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCt                                    <  1:2976704/150‑1 (MQ=255)
                                                                                                                    ttGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTa                                   >  2:1049111/1‑150 (MQ=255)
                                                                                                                     tGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCCGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTAc                                  >  1:2053666/1‑150 (MQ=255)
                                                                                                                       cTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTAccc                                >  2:793149/1‑150 (MQ=255)
                                                                                                                         cAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCct                              >  1:1603861/1‑150 (MQ=255)
                                                                                                                         cAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCct                              >  1:815308/1‑150 (MQ=255)
                                                                                                                          aGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCctg                             >  2:1014834/1‑150 (MQ=255)
                                                                                                                           gCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCctgc                            <  1:4112938/150‑1 (MQ=255)
                                                                                                                             gTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCctgctg                          >  1:2824211/1‑150 (MQ=255)
                                                                                                                              tgctgcTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTgg                         >  2:3882521/1‑150 (MQ=255)
                                                                                                                                 tgctgTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGTTACCCCTGCTGGTAc                      >  2:4119216/1‑150 (MQ=255)
                                                                                                                                 tgctgTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTAc                      >  2:2977918/1‑150 (MQ=255)
                                                                                                                                 tgctgTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTAc                      <  1:2328888/150‑1 (MQ=255)
                                                                                                                                   ctgTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACAt                    >  2:1821666/1‑150 (MQ=255)
                                                                                                                                      tCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCag                 <  2:669837/150‑1 (MQ=255)
                                                                                                                                          tGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCagagag             >  2:2313155/1‑150 (MQ=255)
                                                                                                                                          tGTTCAGAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCagagag             <  2:2301173/150‑1 (MQ=255)
                                                                                                                                            ttCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGAt           >  1:867994/1‑150 (MQ=255)
                                                                                                                                              cATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTa         <  2:1694700/150‑1 (MQ=255)
                                                                                                                                              cATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTa         <  2:4462349/150‑1 (MQ=255)
                                                                                                                                              cATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGACCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGAGTa         >  2:3275368/1‑150 (MQ=255)
                                                                                                                                                tAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAAc       >  2:2935733/1‑150 (MQ=255)
                                                                                                                                                 aaaCGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACt      >  1:4042009/1‑150 (MQ=255)
                                                                                                                                                 aaaCGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACt      >  1:4931072/1‑150 (MQ=255)
                                                                                                                                                 aaaCGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACt      >  2:2996064/1‑150 (MQ=255)
                                                                                                                                                    cGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGTGCGAGAGGCAGTAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACTTCAGAGAGATTAACTGGt   <  1:4046884/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTC  >  CP009273/305506‑305804

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: