Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,814,302 T→A 38.3% S77C (AGC→TGC)  ligB ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,814,3020TA38.3% 222.0 / 300.4 235S77C (AGC→TGC) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (72/73);  new base A (49/41);  total (121/114)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 248 total aligned reads displayed.

ATCAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTG  >  CP009273/3814153‑3814451
                                                                                                                                                     |                                                                                                                                                     
aTCAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGgctcgct                                                                                                                                                       >  1:3697591/1‑150 (MQ=255)
  cAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTcc                                                                                                                                                     >  2:1550466/1‑150 (MQ=255)
    aCTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCa                                                                                                                                                   >  1:2239784/1‑150 (MQ=255)
    aCTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCa                                                                                                                                                   >  2:2477013/1‑150 (MQ=255)
          ggCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAagcag                                                                                                                                             <  2:2280230/150‑1 (MQ=255)
           gCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAagcagc                                                                                                                                            <  1:2471019/150‑1 (MQ=255)
               cACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTg                                                                                                                                        <  1:1899824/150‑1 (MQ=255)
                aCCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGc                                                                                                                                       <  2:3980595/150‑1 (MQ=255)
                 cccAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGcc                                                                                                                                      >  1:3231396/1‑150 (MQ=255)
                 cccAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGcc                                                                                                                                      <  1:1876666/150‑1 (MQ=255)
                 cccAAAGATCGCTACGTTCTCGCATCCACAGACTTAAAGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGcc                                                                                                                                      >  2:3651902/1‑150 (MQ=255)
                  ccAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCa                                                                                                                                     >  2:1130343/1‑150 (MQ=255)
                      aGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGc                                                                                                                                 <  1:150574/150‑1 (MQ=255)
                           gCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAg                                                                                                                            >  2:3302662/1‑150 (MQ=255)
                           gCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAg                                                                                                                            >  1:1667071/1‑150 (MQ=255)
                           gCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAg                                                                                                                            >  1:1941617/1‑150 (MQ=255)
                             tACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGAc                                                                                                                          <  2:1752337/150‑1 (MQ=255)
                             tACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCAATTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGACACTGCGTAAGAc                                                                                                                          >  1:482521/1‑150 (MQ=255)
                                gTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACggg                                                                                                                       >  1:2921276/1‑150 (MQ=255)
                                gTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACggg                                                                                                                       >  1:3744936/1‑150 (MQ=255)
                                 ttCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgc                                                                                                                      >  2:3210057/1‑150 (MQ=255)
                                   ctcGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgcgc                                                                                                                    >  1:3003388/1‑150 (MQ=255)
                                    tcGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCt                                                                                                                   <  2:532793/150‑1 (MQ=255)
                                     cGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCtt                                                                                                                  <  2:89955/150‑1 (MQ=255)
                                       cATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGTGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTaa                                                                                                                <  1:1627091/150‑1 (MQ=255)
                                       cATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCACACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTaa                                                                                                                <  1:1229744/150‑1 (MQ=255)
                                         tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACt                                                                                                              >  1:3033013/1‑150 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             >  1:3656155/1‑150 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             >  1:692544/1‑150 (MQ=255)
                                           cacaGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGa                                                                                                            <  1:4126919/150‑1 (MQ=255)
                                             caGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATc                                                                                                          <  1:3368094/150‑1 (MQ=255)
                                              aGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCg                                                                                                         <  1:3067870/150‑1 (MQ=255)
                                               gACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGt                                                                                                        >  2:3662841/1‑150 (MQ=255)
                                                 cTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTaa                                                                                                      >  1:1041132/1‑150 (MQ=255)
                                                 cTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTaa                                                                                                      <  2:1865700/150‑1 (MQ=255)
                                                   tAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAac                                                                                                    <  1:678469/150‑1 (MQ=255)
                                                   tAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAac                                                                                                    <  1:3415788/150‑1 (MQ=255)
                                                       gCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCg                                                                                                <  1:3055929/150‑1 (MQ=255)
                                                           ttttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCtt                                                                                            <  1:3738106/150‑1 (MQ=255)
                                                             ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTcc                                                                                          >  2:3381549/1‑150 (MQ=255)
                                                              tATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCa                                                                                         <  1:3940539/150‑1 (MQ=255)
                                                              tATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCa                                                                                         >  2:2400746/1‑150 (MQ=255)
                                                               aTCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCAc                                                                                        >  2:3942008/1‑150 (MQ=255)
                                                                tCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCAcc                                                                                       <  1:1122118/150‑1 (MQ=255)
                                                                 ccaccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCt                                                                                      >  1:1160871/1‑150 (MQ=255)
                                                                  caccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTc                                                                                     >  2:3933282/1‑150 (MQ=255)
                                                                  caccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTc                                                                                     <  1:2798399/150‑1 (MQ=255)
                                                                      aTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACtt                                                                                 >  1:3546263/1‑150 (MQ=255)
                                                                           gCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcc                                                                            <  2:3050881/150‑1 (MQ=255)
                                                                            cGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcct                                                                           <  2:203513/150‑1 (MQ=255)
                                                                             gTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcctt                                                                          <  2:232658/150‑1 (MQ=255)
                                                                                cacCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttcct                                                                       <  2:2699600/150‑1 (MQ=255)
                                                                                 acCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttcctt                                                                      <  2:1490272/150‑1 (MQ=255)
                                                                                 acCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttcctt                                                                      <  2:2167491/150‑1 (MQ=255)
                                                                                    cGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttcaa                                                                   <  1:3326537/150‑3 (MQ=255)
                                                                                       aTGGACAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCagta                                                                <  2:2882245/150‑1 (MQ=255)
                                                                                         gggCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCagtagt                                                              >  2:1415911/1‑150 (MQ=255)
                                                                                          ggCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTc                                                             >  2:1527690/1‑150 (MQ=255)
                                                                                           gCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCa                                                            <  1:2707313/150‑1 (MQ=255)
                                                                                            cAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCAt                                                           >  2:43811/1‑150 (MQ=255)
                                                                                            cAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCAt                                                           <  2:30930/150‑1 (MQ=255)
                                                                                            cAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCAt                                                           >  1:1771753/1‑150 (MQ=255)
                                                                                             aaCCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATc                                                          <  1:3912749/150‑1 (MQ=255)
                                                                                             aaCCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATc                                                          <  2:1639363/150‑1 (MQ=255)
                                                                                                 ggATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTcc                                                      >  2:2951029/1‑150 (MQ=255)
                                                                                                 ggATGCATAACAGCATCATTTAAAGGCGGCATCATCACGTCGAGGGGCTCGCTCCCAAAGCAGCGTTGACACTGAGTAAGACGGGCGCTTAACTGATCGTAAACACCGTATTCCACCTCACTTTTTCCTTCATTCCAGTAGTCATCGTcc                                                      <  2:1967155/150‑1 (MQ=255)
                                                                                                  gATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTccc                                                     <  1:2048599/150‑1 (MQ=255)
                                                                                                     gCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACt                                                  >  1:2635394/1‑150 (MQ=255)
                                                                                                       aTAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGt                                                >  2:1134679/1‑150 (MQ=255)
                                                                                                        tAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGtt                                               >  1:3403126/1‑150 (MQ=255)
                                                                                                         aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGttt                                              <  1:3470515/150‑1 (MQ=255)
                                                                                                          aCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGtttt                                             <  1:3897717/150‑1 (MQ=255)
                                                                                                             gCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAtt                                          >  1:2348762/1‑150 (MQ=255)
                                                                                                             gCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAtt                                          >  2:2445432/1‑150 (MQ=255)
                                                                                                             gCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAtt                                          >  2:2671800/1‑150 (MQ=255)
                                                                                                                   tttAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGt                                    <  1:873088/150‑1 (MQ=255)
                                                                                                                    ttAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGtt                                   <  1:1930341/150‑1 (MQ=255)
                                                                                                                            ggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCgg                           >  1:4051407/1‑150 (MQ=255)
                                                                                                                            ggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCgg                           >  1:313042/1‑150 (MQ=255)
                                                                                                                            ggcATCATCACGTCGCGGGGCTCACTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCgg                           <  2:1386357/150‑1 (MQ=255)
                                                                                                                                 catcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAt                      <  1:4134294/150‑1 (MQ=255)
                                                                                                                                 catcaCGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAt                      <  1:2418724/150‑1 (MQ=255)
                                                                                                                                 catcaCGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAt                      >  1:1709604/1‑150 (MQ=255)
                                                                                                                                  atcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAtt                     >  2:1907536/1‑150 (MQ=255)
                                                                                                                                     aCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttct                  <  2:1062485/150‑1 (MQ=255)
                                                                                                                                     aCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttct                  <  2:3399198/150‑1 (MQ=255)
                                                                                                                                       gTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttcttc                >  2:69684/1‑150 (MQ=255)
                                                                                                                                       gTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttcttc                >  2:3283741/1‑150 (MQ=255)
                                                                                                                                       gTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttcttc                >  1:758466/1‑150 (MQ=255)
                                                                                                                                       gTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttcttc                <  1:1476488/150‑1 (MQ=255)
                                                                                                                                           cgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgt            >  2:2801748/1‑150 (MQ=255)
                                                                                                                                           cgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgt            <  2:3410580/150‑1 (MQ=255)
                                                                                                                                             gggCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGc          >  1:139187/1‑150 (MQ=255)
                                                                                                                                             gggCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGc          <  1:2477013/150‑1 (MQ=255)
                                                                                                                                                ctcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGctct       >  1:2162338/1‑150 (MQ=255)
                                                                                                                                                 tcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctg      <  2:3839982/150‑1 (MQ=255)
                                                                                                                                                 tcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctg      >  2:3715029/1‑150 (MQ=255)
                                                                                                                                                 tcgcacCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTAAAATGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCAACTTCACACTGTTTTATTTGATTTTGCAGGCGGGAAATTAATTACTTTGCTctg      >  2:4022800/1‑150 (MQ=255)
                                                                                                                                                  cgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctgg     <  1:472300/150‑1 (MQ=255)
                                                                                                                                                   gctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctggc    <  2:1818678/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTG  >  CP009273/3814153‑3814451

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: