Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,788,981 A→T 5.5% S143R (AGT→AGA btuC ← vitamin B12 transporter subunit: membrane component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,788,9810AT 5.5% 853.2 / 27.8 256S143R (AGT→AGAbtuCvitamin B12 transporter subunit: membrane component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (134/108);  new base T (7/7);  total (141/115)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 257 total aligned reads displayed.

AGCCAGCTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTAG  >  CP009273/1788832‑1789129
                                                                                                                                                     |                                                                                                                                                    
agccagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGa                                                                                                                                                      >  1:1251043/1‑150 (MQ=255)
  ccagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACt                                                                                                                                                    >  1:1111964/1‑150 (MQ=255)
  ccagcTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACt                                                                                                                                                    >  2:4260558/1‑150 (MQ=255)
       tttGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATATTCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCg                                                                                                                                               <  1:3698680/150‑1 (MQ=255)
        ttGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGa                                                                                                                                              >  1:4437330/1‑150 (MQ=255)
           ccgccAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAg                                                                                                                                           >  2:1250501/1‑150 (MQ=255)
               cAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATga                                                                                                                                       <  2:4147307/150‑1 (MQ=255)
                aGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATgac                                                                                                                                      <  2:22126/150‑1 (MQ=255)
                 gTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATgacg                                                                                                                                     <  2:2407187/150‑1 (MQ=255)
                   cTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATgacgac                                                                                                                                   >  1:973082/1‑150 (MQ=255)
                    tACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGacgacg                                                                                                                                  >  2:4075963/1‑150 (MQ=255)
                    tACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGacgacg                                                                                                                                  >  1:3638325/1‑150 (MQ=255)
                    tACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGacgacg                                                                                                                                  >  1:2540546/1‑150 (MQ=255)
                    tACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGacgacg                                                                                                                                  >  1:2539453/1‑150 (MQ=255)
                    tACGCCGCCAAAACCGCCCATCATCCAGGACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGacgacg                                                                                                                                  <  2:2324355/150‑1 (MQ=255)
                     aCGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACgg                                                                                                                                 >  2:3651790/1‑150 (MQ=255)
                      cgccgcCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCTACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGATACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACggg                                                                                                                                <  1:351986/150‑1 (MQ=255)
                       gccgccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGc                                                                                                                               <  2:3367093/150‑1 (MQ=255)
                         cgccAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGa                                                                                                                             >  2:954729/1‑150 (MQ=255)
                            cAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAAc                                                                                                                          <  1:2646224/150‑1 (MQ=255)
                             aaaaCCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGTAAAAGTTGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACg                                                                                                                         <  1:1178779/150‑1 (MQ=255)
                             aaaaCCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCATCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACg                                                                                                                         <  2:2025108/150‑1 (MQ=255)
                               aaCCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTa                                                                                                                       <  1:2847129/150‑1 (MQ=255)
                                 ccGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAag                                                                                                                     <  2:398000/150‑1 (MQ=255)
                                 ccGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAag                                                                                                                     <  2:397906/150‑1 (MQ=255)
                                    cccATCATCCAGTACATCCGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGt                                                                                                                  <  1:1788446/150‑1 (MQ=255)
                                    cccATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGt                                                                                                                  <  1:1347079/150‑1 (MQ=255)
                                     ccATCATCCAGTACATCAGCTGACGCACAGCACCTGTGGTGGAACAGTAGATAGCACACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTa                                                                                                                 <  2:4409972/150‑1 (MQ=255)
                                       atcatcCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTAtt                                                                                                               <  1:2921775/150‑1 (MQ=255)
                                          atcCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTaaa                                                                                                            <  1:3039063/150‑1 (MQ=255)
                                          atcCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTaaa                                                                                                            >  2:590045/1‑150 (MQ=255)
                                          atcCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTaaa                                                                                                            >  2:2608641/1‑150 (MQ=255)
                                            ccAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGt                                                                                                          >  1:1352559/1‑150 (MQ=255)
                                            ccAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGt                                                                                                          >  1:2104441/1‑150 (MQ=255)
                                            ccAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGt                                                                                                          >  1:2104405/1‑150 (MQ=255)
                                                  cATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTgatgat                                                                                                    <  2:1370438/150‑1 (MQ=255)
                                                   aTCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATa                                                                                                   >  1:4176626/1‑150 (MQ=255)
                                                   aTCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATa                                                                                                   <  1:2745847/150‑1 (MQ=255)
                                                   aTCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATa                                                                                                   >  1:3342140/1‑150 (MQ=255)
                                                     cAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAg                                                                                                 >  2:1969104/1‑150 (MQ=255)
                                                     cAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAg                                                                                                 >  2:2640761/1‑150 (MQ=255)
                                                       gCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgcg                                                                                               >  1:61728/1‑150 (MQ=255)
                                                        cTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgcgc                                                                                              >  1:3461235/1‑150 (MQ=255)
                                                            cGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCa                                                                                          <  1:4513764/150‑1 (MQ=255)
                                                              cAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGc                                                                                        >  2:226924/1‑150 (MQ=255)
                                                              cAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGAAAAGCGCGCCAGc                                                                                        >  1:4091519/1‑150 (MQ=255)
                                                                aaTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                                                                                      >  2:1514499/1‑150 (MQ=255)
                                                                 aTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAAt                                                                                     >  2:982987/1‑150 (MQ=255)
                                                                 aTCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAAt                                                                                     <  1:797441/150‑1 (MQ=255)
                                                                    aaCTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGc                                                                                  >  1:1492720/1‑150 (MQ=255)
                                                                         aGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGcacaca                                                                             >  1:2398402/1‑150 (MQ=255)
                                                                         aGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCAAGCAATCGcacaca                                                                             >  1:93246/1‑150 (MQ=255)
                                                                          ggTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAg                                                                            >  1:3097997/1‑150 (MQ=255)
                                                                          ggTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAg                                                                            >  1:3939203/1‑150 (MQ=255)
                                                                            tGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGcc                                                                          >  2:2339572/1‑150 (MQ=255)
                                                                            tGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGcc                                                                          >  2:2339467/1‑150 (MQ=255)
                                                                              gAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCt                                                                        >  1:2347578/1‑150 (MQ=255)
                                                                                 aaGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgc                                                                     >  2:119573/1‑150 (MQ=255)
                                                                                 aaGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgc                                                                     <  2:3501244/150‑1 (MQ=255)
                                                                                  aGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgcg                                                                    >  2:9213/1‑150 (MQ=255)
                                                                                  aGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgcg                                                                    >  1:3902183/1‑150 (MQ=255)
                                                                                       aTAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAg                                                               >  2:1590370/1‑150 (MQ=255)
                                                                                        tAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGt                                                              >  2:2413259/1‑150 (MQ=255)
                                                                                        tAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGt                                                              >  2:2413128/1‑150 (MQ=255)
                                                                                              aCGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCACTCGCACACAGCCCTAGCGCCCAGTTggggc                                                        >  2:3251840/1‑149 (MQ=255)
                                                                                                  cATTAGTGCGCTACAGATAATCCCTACTGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttg                                                    <  2:1479640/150‑1 (MQ=255)
                                                                                                  cATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttg                                                    <  2:1480003/150‑1 (MQ=255)
                                                                                                      aGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttgccct                                                >  1:18372/1‑150 (MQ=255)
                                                                                                        tGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttgcccttg                                              <  2:1977394/150‑1 (MQ=255)
                                                                                                             tACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCaa                                         >  2:3915108/1‑150 (MQ=255)
                                                                                                                  aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGttgccctt                                               >  1:95313/1‑139 (MQ=255)
                                                                                                                  aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAAt                                    <  1:1486077/150‑1 (MQ=255)
                                                                                                                  aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAAt                                    <  2:1119585/150‑1 (MQ=255)
                                                                                                                    aaTCCCTAATGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATAc                                  <  1:2685349/150‑1 (MQ=255)
                                                                                                                    aaTCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGAGCGCCAGCAATCGCACACAGCCATAGCGCCCAGTTTGGGAGTTTCCCTTGCCCAAGCAATAc                                  >  1:2987176/1‑150 (MQ=255)
                                                                                                                      tCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcg                                >  1:1016336/1‑150 (MQ=255)
                                                                                                                      tCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcg                                >  1:3293221/1‑150 (MQ=255)
                                                                                                                       cccTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcgc                               >  1:2956337/1‑150 (MQ=255)
                                                                                                                        ccTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcgcg                              >  2:3767437/1‑150 (MQ=255)
                                                                                                                        ccTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACcgcg                              >  1:1581288/1‑150 (MQ=255)
                                                                                                                          tAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCgcggc                            <  1:2640761/150‑1 (MQ=255)
                                                                                                                          tAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCgcggc                            <  2:768057/150‑1 (MQ=255)
                                                                                                                           aaTGCAACGCCAGCCAGCAATAACCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCgcggcg                           <  1:4075959/150‑1 (MQ=255)
                                                                                                                             tGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGAt                         >  1:1315196/1‑150 (MQ=255)
                                                                                                                                    ccagccagCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGcccc                  <  2:3422073/150‑1 (MQ=255)
                                                                                                                                     cagccagcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCa                 >  1:3887667/1‑150 (MQ=255)
                                                                                                                                      agccagcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCAc                <  2:2813885/150‑1 (MQ=255)
                                                                                                                                        ccagcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGc              >  1:3514598/1‑150 (MQ=255)
                                                                                                                                          agcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCt            >  1:243113/1‑150 (MQ=255)
                                                                                                                                           gcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTg           <  1:4254946/150‑1 (MQ=255)
                                                                                                                                            cAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTgc          <  1:2036066/150‑1 (MQ=255)
                                                                                                                                             aaTAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTgcg         <  2:4176627/150‑1 (MQ=255)
                                                                                                                                             aaTAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTgcg         >  2:1663198/1‑150 (MQ=255)
                                                                                                                                              aTAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTgcgc        <  2:3075114/150‑1 (MQ=255)
                                                                                                                                                aaCCGTCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGTGATAAGCCCCACGCCTGCGCCg      <  2:2127503/150‑1 (MQ=255)
                                                                                                                                                  ccGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGtt    <  1:2702728/150‑1 (MQ=255)
                                                                                                                                                   cGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTa   <  2:583756/150‑1 (MQ=255)
                                                                                                                                                   cGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTa   >  1:3149288/1‑150 (MQ=255)
                                                                                                                                                   cGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTa   <  2:2850394/150‑1 (MQ=255)
                                                                                                                                                    gACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTag  <  2:672413/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCCAGCTTTGCCGCCAGTCTACGCCGCCAAAACCGCCCATCATCCAGTACATCAGCTGACGCAAATCAACTGAGGTGGAAAAGTAGATAGCCCACGTCATTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTGGGGAGTTGCCCTTGCCCAAGCAATACCGCGGCGATAAGCCCCACGCCTGCGCCGTTAG  >  CP009273/1788832‑1789129

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: