Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,814,302 T→A 37.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,3020TA37.4% 373.1 / inf 342S77C (AGC→TGC) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (110/104);  new base A (64/64);  total (174/168)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 359 total aligned reads displayed.

ATCAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTG  >  CP009273/3814153‑3814451
                                                                                                                                                     |                                                                                                                                                     
     cTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCaa                                                                                                                                                  >  2:4347048/1‑150 (MQ=255)
      ttttGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCaaa                                                                                                                                                 >  1:4364373/1‑150 (MQ=255)
       tttGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAag                                                                                                                                                <  2:3483192/150‑1 (MQ=255)
            cTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCg                                                                                                                                           >  1:2667551/1‑150 (MQ=255)
              gCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGtt                                                                                                                                         <  1:224983/150‑1 (MQ=255)
                   cAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCAc                                                                                                                                    <  2:3929815/150‑1 (MQ=255)
                   cAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCAc                                                                                                                                    >  1:2027533/1‑150 (MQ=255)
                      aGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGc                                                                                                                                 >  2:3018646/1‑150 (MQ=255)
                      aGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGc                                                                                                                                 >  1:4161194/1‑150 (MQ=255)
                        aTCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGt                                                                                                                               >  2:1322169/1‑150 (MQ=255)
                        aTCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGt                                                                                                                               >  1:1090747/1‑150 (MQ=255)
                          cGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTaa                                                                                                                             >  1:4480742/1‑150 (MQ=255)
                          cGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTaa                                                                                                                             <  2:430362/150‑1 (MQ=255)
                           gCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAg                                                                                                                            <  1:1825825/150‑1 (MQ=255)
                           gCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAg                                                                                                                            >  2:3741556/1‑150 (MQ=255)
                           gCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAg                                                                                                                            <  2:3714305/150‑1 (MQ=255)
                               cGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACgg                                                                                                                        <  2:2015868/150‑1 (MQ=255)
                                 ttCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgc                                                                                                                      >  2:1063596/1‑150 (MQ=255)
                                 ttCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgc                                                                                                                      <  1:2006091/150‑1 (MQ=255)
                                   ctcGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGgcgc                                                                                                                    >  1:3612043/1‑150 (MQ=255)
                                     cGCATCCACAGACTTAATGCGTTTTTCTCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCACAGCAGCGTTGCCACTGCGTAAGACGGGCGCtt                                                                                                                  <  1:3020139/150‑1 (MQ=255)
                                     cGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCtt                                                                                                                  >  1:3566812/1‑150 (MQ=255)
                                      gCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTa                                                                                                                 <  1:809124/150‑1 (MQ=255)
                                         tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACt                                                                                                              >  2:4069543/1‑150 (MQ=255)
                                         tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACt                                                                                                              <  1:894691/150‑1 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             >  2:4124669/1‑150 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             >  1:279733/1‑150 (MQ=255)
                                           cacaGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGa                                                                                                            <  2:2045098/150‑1 (MQ=255)
                                            acaGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGAt                                                                                                           >  2:2651392/1‑150 (MQ=255)
                                             caGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATc                                                                                                          <  1:3044971/150‑1 (MQ=255)
                                                 cTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTaa                                                                                                      >  1:2999797/1‑150 (MQ=255)
                                                 cTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTaa                                                                                                      >  1:4245351/1‑150 (MQ=255)
                                                        cGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCATCCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGAGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGt                                                                                               <  1:2662529/150‑1 (MQ=255)
                                                            tttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTc                                                                                           <  2:3084155/150‑1 (MQ=255)
                                                            tttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTc                                                                                           <  1:2626792/150‑1 (MQ=255)
                                                             ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTcc                                                                                          >  2:4192296/1‑150 (MQ=255)
                                                             ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTcc                                                                                          >  2:1430310/1‑150 (MQ=255)
                                                              tATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCa                                                                                         <  1:1659900/150‑1 (MQ=255)
                                                                  caccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTc                                                                                     >  2:1865695/1‑150 (MQ=255)
                                                                  caccaTTTTGCGTACACCCGAATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTc                                                                                     >  2:1570025/1‑150 (MQ=255)
                                                                   accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                    >  2:3964443/1‑150 (MQ=255)
                                                                   accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                    >  1:1776350/1‑150 (MQ=255)
                                                                   accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                    >  2:2406775/1‑150 (MQ=255)
                                                                    ccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCAc                                                                                   >  1:2911851/1‑150 (MQ=255)
                                                                       ttttGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACttt                                                                                <  2:618607/150‑1 (MQ=255)
                                                                          tGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttc                                                                             >  1:3954490/1‑150 (MQ=255)
                                                                           gCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcc                                                                            <  2:180346/150‑1 (MQ=255)
                                                                           gCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcc                                                                            <  2:179836/150‑1 (MQ=255)
                                                                             gTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcctt                                                                          >  1:4798485/1‑150 (MQ=255)
                                                                                  cccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttc                                                                     >  2:2227208/1‑150 (MQ=255)
                                                                                   ccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttcc                                                                    <  2:3166384/150‑1 (MQ=255)
                                                                                    cGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCa                                                                   <  1:4528434/150‑1 (MQ=255)
                                                                                    cGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCa                                                                   >  2:2432438/1‑150 (MQ=255)
                                                                                    cGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCa                                                                   <  1:1469807/150‑1 (MQ=255)
                                                                                     gTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCag                                                                  >  2:798879/1‑150 (MQ=255)
                                                                                        tGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCagtag                                                               >  1:4666199/1‑150 (MQ=255)
                                                                                          ggCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTc                                                             <  2:1216759/150‑1 (MQ=255)
                                                                                           gCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCa                                                            <  1:2091541/150‑1 (MQ=255)
                                                                                           gCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCa                                                            <  2:1247096/150‑1 (MQ=255)
                                                                                            cAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCAt                                                           <  1:3723455/150‑1 (MQ=255)
                                                                                             aaCCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATc                                                          >  1:601092/1‑150 (MQ=255)
                                                                                                  gATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTccc                                                     <  2:4567816/150‑1 (MQ=255)
                                                                                                   aTGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCa                                                    <  1:4069543/150‑1 (MQ=255)
                                                                                                   aTGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCa                                                    <  1:3394175/150‑1 (MQ=255)
                                                                                                     gCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACt                                                  <  1:3279033/150‑1 (MQ=255)
                                                                                                       aTAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGt                                                >  2:3336813/1‑150 (MQ=255)
                                                                                                         aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGttt                                              <  2:2718946/150‑1 (MQ=255)
                                                                                                           ccGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTa                                            <  1:1002404/150‑1 (MQ=255)
                                                                                                            cGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAt                                           <  2:1727680/150‑1 (MQ=255)
                                                                                                            cGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAt                                           <  2:1090747/150‑1 (MQ=255)
                                                                                                            cGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAt                                           >  2:3433037/1‑150 (MQ=255)
                                                                                                            cGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAt                                           >  2:2510405/1‑150 (MQ=255)
                                                                                                             gCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAtt                                          <  2:1032624/150‑1 (MQ=255)
                                                                                                              cACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTAttt                                         >  2:465145/1‑150 (MQ=255)
                                                                                                                  aTTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTg                                     <  2:3300577/150‑1 (MQ=255)
                                                                                                                   tttAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCACACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGt                                    >  1:972331/1‑150 (MQ=255)
                                                                                                                     tAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACAACGTCTTCCACCTCACATTTTCCTTCCTTACAGTAGTAATAGTCCCACTGGTTTATTTTCTGTTg                                  >  1:1228782/1‑150 (MQ=255)
                                                                                                                     tAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTg                                  <  2:1601330/150‑1 (MQ=255)
                                                                                                                       aaGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCa                                <  2:2559819/150‑1 (MQ=255)
                                                                                                                       aaGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCa                                >  1:2346901/1‑150 (MQ=255)
                                                                                                                        aGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAg                               <  1:2186817/150‑1 (MQ=255)
                                                                                                                         ggcggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAgg                              >  1:4651837/1‑150 (MQ=255)
                                                                                                                          gcggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGc                             <  2:394051/150‑1 (MQ=255)
                                                                                                                           cggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCg                            <  1:4119234/150‑1 (MQ=255)
                                                                                                                           cggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCg                            >  1:2198397/1‑150 (MQ=255)
                                                                                                                               atcatcACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGaa                        >  1:3834685/1‑150 (MQ=255)
                                                                                                                                  atcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAtt                     <  1:256954/150‑1 (MQ=255)
                                                                                                                                   tcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAttt                    <  2:4425221/150‑1 (MQ=255)
                                                                                                                                   tcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAttt                    <  1:2144087/150‑1 (MQ=255)
                                                                                                                                    caCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttc                   >  1:3170947/1‑150 (MQ=255)
                                                                                                                                     aCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGACATttct                  >  2:4866327/1‑150 (MQ=255)
                                                                                                                                         cgcgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCt              <  1:1017471/150‑1 (MQ=255)
                                                                                                                                         cgcgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCt              >  2:1998872/1‑150 (MQ=255)
                                                                                                                                          gcgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtg             <  1:3418223/150‑1 (MQ=255)
                                                                                                                                          gcgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtg             <  2:1824627/150‑1 (MQ=255)
                                                                                                                                           cgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgt            <  1:3528321/150‑1 (MQ=255)
                                                                                                                                            ggggCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgtg           <  2:2814874/150‑1 (MQ=255)
                                                                                                                                              ggCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGct         <  2:3670883/150‑1 (MQ=255)
                                                                                                                                                ctcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGctct       <  1:1703077/150‑1 (MQ=255)
                                                                                                                                                  cgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTGTGCCTTCCGTCCAGTAGTCAGCGTCCCACTGTGTGATTTGCTGTTGCAGGCGGGCAAGGTCTTCCTGTGCGctgg     >  2:4815525/1‑150 (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: