Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,129,363 G→A 17.3% A36T (GCG→ACG)  flgF → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,129,3630GA17.3% 869.0 / 154.3 329A36T (GCG→ACG) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base G (155/115);  new base A (32/25);  total (189/140)
Fisher's exact test for biased strand distribution p-value = 8.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

CAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGAT  >  CP009273/1129214‑1129512
                                                                                                                                                     |                                                                                                                                                     
cAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCCATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcg                                                                                                                                                       >  2:1183186/1‑150 (MQ=255)
cAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcg                                                                                                                                                       >  2:2619303/1‑150 (MQ=255)
      gCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGt                                                                                                                                                 >  1:2287845/1‑150 (MQ=255)
       cTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGtt                                                                                                                                                <  1:3473426/150‑1 (MQ=255)
       cTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGtt                                                                                                                                                <  2:3153190/150‑1 (MQ=255)
           ttAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAAt                                                                                                                                            <  2:5468997/150‑1 (MQ=255)
            tAACGTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATg                                                                                                                                           <  1:2678617/150‑1 (MQ=255)
               cTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCtt                                                                                                                                        <  1:583160/150‑1 (MQ=255)
                 tACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTa                                                                                                                                      >  2:4129246/1‑150 (MQ=255)
                     cTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                                                                                  <  2:3423278/150‑1 (MQ=255)
                      tAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgc                                                                                                                                 >  1:4596797/1‑150 (MQ=255)
                      tAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgc                                                                                                                                 >  1:3596082/1‑150 (MQ=255)
                       aaTCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgcg                                                                                                                                >  2:377537/1‑150 (MQ=255)
                         ccGCTGACGCGAGAGCCCACTCGCTCACGCACTCGAGACCGCGAAGGGAGCAGCCAGCCACACACTGAATCAACAGGCGGTAACCGCCAGTAATCAGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGt                                                                                                                              <  2:4426400/149‑1 (MQ=255)
                          cGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTg                                                                                                                             <  1:5267069/150‑1 (MQ=255)
                            cTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGcc                                                                                                                           <  1:3209597/150‑1 (MQ=255)
                                  ggATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGa                                                                                                                     >  1:866368/1‑150 (MQ=255)
                                    aTAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAg                                                                                                                   >  1:281744/1‑150 (MQ=255)
                                    aTAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAg                                                                                                                   >  1:3141913/1‑150 (MQ=255)
                                    aTAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAg                                                                                                                   >  2:218661/1‑150 (MQ=255)
                                    aTAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAg                                                                                                                   >  1:3953666/1‑150 (MQ=255)
                                     tAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAgg                                                                                                                  >  1:5669778/1‑150 (MQ=255)
                                       gCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGc                                                                                                                <  2:389170/150‑1 (MQ=255)
                                       gCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGc                                                                                                                <  2:5566534/150‑1 (MQ=255)
                                          cAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCttt                                                                                                             <  1:3725628/150‑1 (MQ=255)
                                          cAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCttt                                                                                                             <  2:1747688/150‑1 (MQ=255)
                                            aTGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtct                                                                                                           <  2:3091593/150‑1 (MQ=255)
                                                 tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGcc                                                                                                      >  2:549181/1‑150 (MQ=255)
                                                    cGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGCGCTTTCTCTGCCCAc                                                                                                   >  2:3501766/1‑150 (MQ=255)
                                                     gCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcg                                                                                                  <  1:4934096/150‑1 (MQ=255)
                                                      cAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcgc                                                                                                 >  1:1120673/1‑150 (MQ=255)
                                                         tttATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCa                                                                                              >  1:3816619/1‑150 (MQ=255)
                                                           tataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACg                                                                                            >  1:2913481/1‑150 (MQ=255)
                                                            ataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                                                                                           >  1:2534681/1‑150 (MQ=255)
                                                            ataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                                                                                           >  1:2000698/1‑150 (MQ=255)
                                                            ataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                                                                                           >  1:4414828/1‑150 (MQ=255)
                                                            ataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                                                                                           >  1:4650575/1‑150 (MQ=255)
                                                            ataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                                                                                           >  2:664890/1‑150 (MQ=255)
                                                            ataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                                                                                           >  2:23892/1‑150 (MQ=255)
                                                             taCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                                                                                          <  1:5358825/150‑1 (MQ=255)
                                                               ccGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                                                                                        <  1:959529/150‑1 (MQ=255)
                                                                cgcgATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGt                                                                                       >  2:770879/1‑150 (MQ=255)
                                                                 gcgATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTc                                                                                      <  2:19598/150‑1 (MQ=255)
                                                                   gATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCAc                                                                                    >  1:3930712/1‑150 (MQ=255)
                                                                    aTGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACg                                                                                   >  2:3679974/1‑150 (MQ=255)
                                                                    aTGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACg                                                                                   >  2:5667555/1‑150 (MQ=255)
                                                                      gggAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTTGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGc                                                                                 <  2:4822528/150‑1 (MQ=255)
                                                                         agcagcCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTc                                                                              <  2:4788485/150‑1 (MQ=255)
                                                                           cagccagcCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCaa                                                                            >  2:2748559/1‑150 (MQ=255)
                                                                           cagccagcCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCaa                                                                            >  1:466804/1‑150 (MQ=255)
                                                                            agccagccAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAc                                                                           <  2:4249798/150‑1 (MQ=255)
                                                                            agccagccAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAc                                                                           <  2:1912800/150‑1 (MQ=255)
                                                                              ccagccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGc                                                                         >  1:65608/1‑150 (MQ=255)
                                                                               cagccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGcc                                                                        <  1:1683783/150‑1 (MQ=255)
                                                                                agccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCg                                                                       >  1:4831754/1‑150 (MQ=255)
                                                                                agccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCg                                                                       >  2:1456717/1‑150 (MQ=255)
                                                                                   cagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgc                                                                    <  1:1224327/150‑1 (MQ=255)
                                                                                    agACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCGTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcg                                                                   >  2:3827220/1‑150 (MQ=255)
                                                                                    agACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcg                                                                   >  2:5072499/1‑150 (MQ=255)
                                                                                    agACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcg                                                                   >  2:297447/1‑150 (MQ=255)
                                                                                         cTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGat                                                              <  2:20074/150‑1 (MQ=255)
                                                                                           gAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGatat                                                            >  1:4145313/1‑150 (MQ=255)
                                                                                           gAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGGTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGatag                                                            <  2:3809851/150‑2 (MQ=255)
                                                                                              tCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGAc                                                         <  1:1656864/150‑1 (MQ=255)
                                                                                               cAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACg                                                        <  2:1350785/150‑1 (MQ=255)
                                                                                               cAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACg                                                        >  2:3572063/1‑150 (MQ=255)
                                                                                                 accggCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGcc                                                      <  2:5008999/147‑1 (MQ=255)
                                                                                                   aGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCg                                                    <  1:5601989/150‑1 (MQ=255)
                                                                                                    ggCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCgg                                                   <  2:4091304/150‑1 (MQ=255)
                                                                                                            ccGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATgg                                           <  2:4040741/150‑1 (MQ=255)
                                                                                                             cGCCAGTAATCTGGCCAAGGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGa                                          <  2:5447555/150‑1 (MQ=255)
                                                                                                              gCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGAt                                         >  1:5065485/1‑150 (MQ=255)
                                                                                                              gCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGAt                                         <  2:1825306/150‑1 (MQ=255)
                                                                                                              gCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGAt                                         >  1:4786395/1‑150 (MQ=255)
                                                                                                              gCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGAt                                         >  1:1954944/1‑150 (MQ=255)
                                                                                                                cAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTa                                       <  2:4687780/150‑1 (MQ=255)
                                                                                                                cAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTa                                       <  1:3520003/150‑1 (MQ=255)
                                                                                                                    aaTCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACAcc                                   <  1:4705516/150‑1 (MQ=255)
                                                                                                                     aTCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCt                                  >  2:340998/1‑150 (MQ=255)
                                                                                                                     aTCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCt                                  <  2:5221269/150‑1 (MQ=255)
                                                                                                                       cTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTcg                                <  2:277616/150‑1 (MQ=255)
                                                                                                                       cTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTcg                                <  2:2546795/150‑1 (MQ=255)
                                                                                                                          gCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCgcgc                             >  1:4356661/1‑150 (MQ=255)
                                                                                                                             aaTGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCg                          >  1:5525196/1‑150 (MQ=255)
                                                                                                                               tGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCt                        <  1:2619303/150‑1 (MQ=255)
                                                                                                                                 ccTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTgg                      >  1:3303402/1‑150 (MQ=255)
                                                                                                                                  cTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGa                     >  1:1920424/1‑150 (MQ=255)
                                                                                                                                  cTCAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGa                     >  1:1435508/1‑150 (MQ=255)
                                                                                                                                    cAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACg                   <  1:3011208/150‑1 (MQ=255)
                                                                                                                                    cAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACg                   <  2:2079301/150‑1 (MQ=255)
                                                                                                                                    cAACGCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACg                   <  1:3421614/150‑1 (MQ=255)
                                                                                                                                     aaCGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCCGTGGAAGGGCTTTCTCTGCCCCCGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAACTGGATTACACCACGCGCCCGCGGGACgg                  >  2:4561409/1‑149 (MQ=255)
                                                                                                                                     aaCGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGt                  >  1:2220003/1‑150 (MQ=255)
                                                                                                                                     aaCGCCCGGTTTTCGCACGCAGTTGAATGCTGTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGt                  <  2:1134037/150‑1 (MQ=255)
                                                                                                                                       cGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTcg                <  2:1300856/150‑1 (MQ=255)
                                                                                                                                        gCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTcgc               >  1:1898480/1‑150 (MQ=255)
                                                                                                                                        gCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTcgc               <  2:3303707/150‑1 (MQ=255)
                                                                                                                                        gCCCGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTcgc               >  1:5503684/1‑150 (MQ=255)
                                                                                                                                          ccGGTTTTCGCACGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGt             >  1:3046484/1‑150 (MQ=255)
                                                                                                                                                   gcgcgcAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAgg    <  2:2499091/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGCAGATATGACGCCCGGCAAAATGGATTACACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGAT  >  CP009273/1129214‑1129512

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: