Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,814,302 T→A 5.4% 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,3020TA 5.4% 756.6 / 23.5 242S77C (AGC→TGC) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (122/107);  new base A (9/4);  total (131/111)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 257 total aligned reads displayed.

ACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGC  >  CP009273/3814157‑3814449
                                                                                                                                                 |                                                                                                                                                   
        cTGCACCCAAAGATCGCTACGTTCTAGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAagcag                                                                                                                                           >  1:232008/1‑148 (MQ=255)
          gCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCg                                                                                                                                         <  1:2785839/148‑1 (MQ=255)
             cccAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTg                                                                                                                                      <  1:3266450/148‑1 (MQ=255)
             cccAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTg                                                                                                                                      <  1:1494830/148‑1 (MQ=255)
                aaaGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCa                                                                                                                                   <  1:2440930/148‑1 (MQ=255)
                aaaGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTGTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCa                                                                                                                                   <  1:2399615/148‑1 (MQ=255)
                      cGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGt                                                                                                                             <  1:3447342/148‑1 (MQ=255)
                        cTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTaa                                                                                                                           >  2:1599049/1‑148 (MQ=255)
                          aCGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGa                                                                                                                         <  2:116284/148‑1 (MQ=255)
                          aCGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGa                                                                                                                         >  2:1812678/1‑148 (MQ=255)
                          aCGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGa                                                                                                                         <  2:3144360/148‑1 (MQ=255)
                              tctcGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACggg                                                                                                                     >  1:3494510/1‑148 (MQ=255)
                              tctcGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACggg                                                                                                                     >  1:2634198/1‑148 (MQ=255)
                                tcGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgcg                                                                                                                   <  1:60567/148‑1 (MQ=255)
                                 cGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgcgc                                                                                                                  >  2:1113784/1‑148 (MQ=255)
                                   cATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCtt                                                                                                                >  1:1699233/1‑148 (MQ=255)
                                     tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTaa                                                                                                              <  1:359364/148‑1 (MQ=255)
                                     tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTaa                                                                                                              <  2:1837823/148‑1 (MQ=255)
                                     tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTaa                                                                                                              >  1:1191961/1‑148 (MQ=255)
                                      ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAAc                                                                                                             >  1:1087542/1‑148 (MQ=255)
                                        acaGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                           <  1:747858/148‑1 (MQ=255)
                                        acaGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                           <  2:1544352/148‑1 (MQ=255)
                                         caGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGa                                                                                                          >  2:2979055/1‑148 (MQ=255)
                                         caGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGa                                                                                                          <  1:298797/148‑1 (MQ=255)
                                             cTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGt                                                                                                      >  2:3010352/1‑148 (MQ=255)
                                             cTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGt                                                                                                      >  2:1231890/1‑148 (MQ=255)
                                              ttAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTa                                                                                                     >  1:2480458/1‑148 (MQ=255)
                                              ttAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTa                                                                                                     >  1:2831062/1‑148 (MQ=255)
                                              ttAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTa                                                                                                     >  1:2860211/1‑148 (MQ=255)
                                              ttAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTa                                                                                                     >  2:2009070/1‑148 (MQ=255)
                                               tAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTaa                                                                                                    <  2:209694/148‑1 (MQ=255)
                                                   gCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAacac                                                                                                >  2:2238957/1‑148 (MQ=255)
                                                     gTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCg                                                                                              <  2:1556318/148‑1 (MQ=255)
                                                      tttttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGt                                                                                             <  1:1275492/148‑1 (MQ=255)
                                                        tttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCt                                                                                           <  1:1539211/148‑1 (MQ=255)
                                                         ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCtt                                                                                          <  2:439452/148‑1 (MQ=255)
                                                         ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCtt                                                                                          <  2:437196/148‑1 (MQ=255)
                                                         ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCtt                                                                                          <  2:1883993/148‑1 (MQ=255)
                                                         ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCtt                                                                                          <  2:861910/148‑1 (MQ=255)
                                                          tATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTc                                                                                         <  2:2022724/148‑1 (MQ=255)
                                                             ccaccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTATACACCGTCTTCCAc                                                                                      <  1:3280506/148‑1 (MQ=255)
                                                             ccaccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCAc                                                                                      >  1:472957/1‑148 (MQ=255)
                                                             ccaccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCAc                                                                                      <  1:2029368/148‑1 (MQ=255)
                                                               accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCt                                                                                    >  1:45002/1‑148 (MQ=255)
                                                                 caTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                  <  1:1138387/148‑1 (MQ=255)
                                                                    tttGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACtt                                                                               <  2:3504553/148‑1 (MQ=255)
                                                                     ttGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACttt                                                                              <  1:3227936/148‑1 (MQ=255)
                                                                          tACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcct                                                                         >  2:1739602/1‑148 (MQ=255)
                                                                          tACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcct                                                                         <  2:1655320/148‑1 (MQ=255)
                                                                              cccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttcct                                                                     <  2:24692/148‑1 (MQ=255)
                                                                                  tATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCa                                                                 >  1:1020530/1‑148 (MQ=255)
                                                                                   aTGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCag                                                                >  2:1397491/1‑148 (MQ=255)
                                                                                      ggCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCagtag                                                             >  1:3202831/1‑148 (MQ=255)
                                                                                        cAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTc                                                           >  2:1496733/1‑148 (MQ=255)
                                                                                        cAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTc                                                           >  1:1436349/1‑148 (MQ=255)
                                                                                         aaCCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCa                                                          <  2:3481326/148‑1 (MQ=255)
                                                                                          aCCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACATTTTCCTTCCTTCCAGTAGTCAt                                                         >  2:698333/1‑148 (MQ=255)
                                                                                           ccGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATc                                                        <  1:1599049/148‑1 (MQ=255)
                                                                                           ccGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATc                                                        <  1:1598792/148‑1 (MQ=255)
                                                                                           ccGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATc                                                        <  1:2698864/148‑1 (MQ=255)
                                                                                            cGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCg                                                       <  2:2363191/148‑1 (MQ=255)
                                                                                               aTGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTcc                                                    <  1:3360354/148‑1 (MQ=255)
                                                                                                 gCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTACGACGGGCGCTTAACTGATCGTAAACACCGTCGTCCACCTCACTGTTTCCGTCCGTCCAGTAGTCAGCggggcc                                                  >  1:436739/1‑143 (MQ=255)
                                                                                                 gCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCa                                                  <  2:661548/148‑1 (MQ=255)
                                                                                                  cATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCAc                                                 >  1:2354256/1‑148 (MQ=255)
                                                                                                     aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGt                                              <  2:1565398/148‑1 (MQ=255)
                                                                                                     aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGt                                              >  2:1498365/1‑148 (MQ=255)
                                                                                                     aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGt                                              <  2:1087542/148‑1 (MQ=255)
                                                                                                       ccGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGttt                                            <  2:2907408/148‑1 (MQ=255)
                                                                                                        cGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGtttt                                           <  2:2790207/148‑1 (MQ=255)
                                                                                                             cATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTg                                      <  1:2911189/148‑1 (MQ=255)
                                                                                                              aTTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGc                                     >  1:2919278/1‑148 (MQ=255)
                                                                                                              aTTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGc                                     <  2:1996965/148‑1 (MQ=255)
                                                                                                               tttAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCt                                    >  2:128299/1‑148 (MQ=255)
                                                                                                                ttAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTg                                   >  1:1385590/1‑148 (MQ=255)
                                                                                                                 tAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGt                                  <  1:542584/148‑1 (MQ=255)
                                                                                                                  aaaGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGtt                                 >  1:1499717/1‑148 (MQ=255)
                                                                                                                     ggcggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCa                              >  1:1855953/1‑148 (MQ=255)
                                                                                                                       cggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAgg                            >  1:720997/1‑148 (MQ=255)
                                                                                                                       cggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAgg                            >  1:842409/1‑148 (MQ=255)
                                                                                                                        ggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCGTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGc                           <  2:325713/148‑1 (MQ=255)
                                                                                                                           atcatcACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCggg                        >  2:2673623/1‑148 (MQ=255)
                                                                                                                              atcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGaaa                     <  1:2562614/148‑1 (MQ=255)
                                                                                                                                caCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAtt                   >  2:2587301/1‑148 (MQ=255)
                                                                                                                                caCGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAtt                   >  2:3103930/1‑148 (MQ=255)
                                                                                                                                  cGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttc                 <  2:3494510/148‑1 (MQ=255)
                                                                                                                                  cGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttc                 >  1:3367589/1‑148 (MQ=255)
                                                                                                                                  cGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttc                 <  2:2634198/148‑1 (MQ=255)
                                                                                                                                    tCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttctt               >  1:980931/1‑148 (MQ=255)
                                                                                                                                    tCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttctt               <  1:3255754/148‑1 (MQ=255)
                                                                                                                                     cgcgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGACATttcttc              >  2:1933166/1‑148 (MQ=255)
                                                                                                                                      gcgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTcc             <  2:2831062/148‑1 (MQ=255)
                                                                                                                                          ggCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCCGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgtg         <  1:3158609/148‑1 (MQ=255)
                                                                                                                                          ggCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgtg         >  2:51164/1‑148 (MQ=255)
                                                                                                                                          ggCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTACTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgtg         >  1:3095324/1‑148 (MQ=255)
                                                                                                                                           gctcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGc        >  2:3351680/1‑148 (MQ=255)
                                                                                                                                            ctcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGct       <  1:122923/148‑1 (MQ=255)
                                                                                                                                                ctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctgg   <  2:1762843/148‑1 (MQ=255)
                                                                                                                                                ctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctgg   >  1:385716/1‑148 (MQ=255)
                                                                                                                                                ctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctgg   >  2:548614/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGC  >  CP009273/3814157‑3814449

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: