Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,814,302 T→A 35.6% 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,3020TA35.6% 387.1 / inf 328S77C (AGC→TGC) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (99/112);  new base A (53/64);  total (152/176)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 344 total aligned reads displayed.

ATCAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTG  >  CP009273/3814153‑3814451
                                                                                                                                                     |                                                                                                                                                     
aTCAACTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGgctcgct                                                                                                                                                       >  2:3122284/1‑150 (MQ=255)
   aaCTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCAccc                                                                                                                                                    >  2:2896186/1‑150 (MQ=255)
    aCTTTTGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCa                                                                                                                                                   >  2:625287/1‑150 (MQ=255)
      ttttGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTACTGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCaaa                                                                                                                                                 <  1:2022097/150‑1 (MQ=255)
       tttGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAag                                                                                                                                                >  2:687030/1‑150 (MQ=255)
         tGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAagca                                                                                                                                              >  1:1580710/1‑150 (MQ=255)
         tGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAagca                                                                                                                                              >  2:1679730/1‑150 (MQ=255)
         tGGCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAagca                                                                                                                                              <  2:3031836/150‑1 (MQ=255)
           gCTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAagcagc                                                                                                                                            <  2:2746313/150‑1 (MQ=255)
            cTGCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCg                                                                                                                                           <  1:568058/150‑1 (MQ=255)
              gCACCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGtt                                                                                                                                         <  2:2406454/150‑1 (MQ=255)
                aCCCAAAGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGc                                                                                                                                       >  1:1655084/1‑150 (MQ=255)
                    aaaGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACt                                                                                                                                   >  1:2789234/1‑150 (MQ=255)
                    aaaGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGAATCATCACGACGCGGGGCTCGCACCCAAAGCAGCGTTGCCACt                                                                                                                                   >  1:3633783/1‑150 (MQ=255)
                     aaGATCGCTACGTTCTCGCATACACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGAATAACCGCACCATATAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGATGCCACTg                                                                                                                                  >  1:791444/1‑150 (MQ=255)
                      aGATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGc                                                                                                                                 <  1:4744486/150‑1 (MQ=255)
                       gATCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCg                                                                                                                                <  1:2980988/150‑1 (MQ=255)
                        aTCGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGt                                                                                                                               >  2:383893/1‑150 (MQ=255)
                          cGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTaa                                                                                                                             <  1:3799553/150‑1 (MQ=255)
                          cGCTACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTaa                                                                                                                             >  1:3770635/1‑150 (MQ=255)
                             tACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGAc                                                                                                                          <  1:4167916/150‑1 (MQ=255)
                             tACGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGAc                                                                                                                          <  2:1330619/150‑1 (MQ=255)
                              aCGTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACg                                                                                                                         <  2:2286610/150‑1 (MQ=255)
                                gTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACggg                                                                                                                       >  2:2151087/1‑150 (MQ=255)
                                gTTCTCGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACggg                                                                                                                       <  2:1163933/150‑1 (MQ=255)
                                     cGCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCtt                                                                                                                  <  1:1683941/150‑1 (MQ=255)
                                      gCATCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTa                                                                                                                 <  1:4471699/150‑1 (MQ=255)
                                        aTCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAAc                                                                                                               >  2:2255692/1‑150 (MQ=255)
                                        aTCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAAc                                                                                                               >  2:1821776/1‑150 (MQ=255)
                                        aTCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAAc                                                                                                               <  2:4776428/150‑1 (MQ=255)
                                         tCCACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACt                                                                                                              >  1:4784094/1‑150 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             <  2:2382478/150‑1 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             >  2:807132/1‑150 (MQ=255)
                                          ccACAGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTg                                                                                                             >  1:1794160/1‑150 (MQ=255)
                                            acaGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGAt                                                                                                           >  1:1441340/1‑150 (MQ=255)
                                             caGACTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATc                                                                                                          <  2:444584/150‑1 (MQ=255)
                                                aCTTAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTa                                                                                                       <  1:3727342/150‑1 (MQ=255)
                                                   tAATGCGTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAac                                                                                                    <  1:4346503/150‑1 (MQ=255)
                                                         gTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTc                                                                                              >  2:1126563/1‑150 (MQ=255)
                                                         gTTTTTATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTc                                                                                              >  1:4122407/1‑150 (MQ=255)
                                                           ttttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCtt                                                                                            <  1:725842/150‑1 (MQ=255)
                                                             ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTcc                                                                                          <  1:1842347/150‑1 (MQ=255)
                                                             ttATCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTcc                                                                                          <  2:3986807/150‑1 (MQ=255)
                                                               aTCCACCATTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCAc                                                                                        <  1:2669067/150‑1 (MQ=255)
                                                                  caccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTc                                                                                     <  2:2020535/150‑1 (MQ=255)
                                                                   accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                    >  2:2639950/1‑150 (MQ=255)
                                                                   accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                    <  1:400840/150‑1 (MQ=255)
                                                                   accaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCa                                                                                    <  1:4572090/150‑1 (MQ=255)
                                                                    ccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCAc                                                                                   >  1:1061548/1‑150 (MQ=255)
                                                                    ccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCAc                                                                                   >  1:3600399/1‑150 (MQ=255)
                                                                    ccaTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCAc                                                                                   <  2:1025276/150‑1 (MQ=255)
                                                                     caTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACt                                                                                  <  1:687026/150‑1 (MQ=255)
                                                                     caTTTTGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATATAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGAAAACACCGTCTTCCACCTCACt                                                                                  >  2:4170679/1‑150 (MQ=255)
                                                                       ttttGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACttt                                                                                <  2:337455/150‑1 (MQ=255)
                                                                       ttttGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACttt                                                                                <  2:4393546/150‑1 (MQ=255)
                                                                         ttGCGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACttttt                                                                              <  1:1304201/150‑1 (MQ=255)
                                                                            cGTACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttcct                                                                           <  2:412342/150‑1 (MQ=255)
                                                                              tACACCCGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttc                                                                         <  2:1857227/150‑1 (MQ=255)
                                                                                  cccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttc                                                                     <  2:3705583/150‑1 (MQ=255)
                                                                                  cccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttc                                                                     <  1:4332780/150‑1 (MQ=255)
                                                                                  cccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttc                                                                     <  2:656630/150‑1 (MQ=255)
                                                                                   ccGTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTttccttccttcc                                                                    >  1:1605116/1‑150 (MQ=255)
                                                                                     gTATGGGCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCag                                                                  <  1:3365136/150‑1 (MQ=255)
                                                                                          ggCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTc                                                             <  1:514162/150‑1 (MQ=255)
                                                                                          ggCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTc                                                             <  1:2269985/150‑1 (MQ=255)
                                                                                           gCAACCGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCa                                                            <  2:4012597/150‑1 (MQ=255)
                                                                                                cGGATGCATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTc                                                       >  2:4558090/1‑150 (MQ=255)
                                                                                                      cATAACCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTg                                                 <  2:942816/150‑1 (MQ=255)
                                                                                                         aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGttt                                              <  1:3182845/150‑1 (MQ=255)
                                                                                                         aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGttt                                              >  2:1982838/1‑150 (MQ=255)
                                                                                                         aaCCGCACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGttt                                              <  2:3693735/150‑1 (MQ=255)
                                                                                                              cACCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTACAGTATTAATCGTCCCACTGTGTTAttt                                         >  2:4079993/1‑150 (MQ=255)
                                                                                                               aCCATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTg                                        >  2:1489734/1‑150 (MQ=255)
                                                                                                                 cATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCt                                      >  1:3701389/1‑150 (MQ=255)
                                                                                                                 cATTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCt                                      >  1:3003658/1‑150 (MQ=255)
                                                                                                                  aTTTAAAGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCtt                                     <  1:3470891/150‑2 (MQ=255)
                                                                                                                      aaaGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGc                                 >  1:4429125/1‑150 (MQ=255)
                                                                                                                       aaGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTGTATTTGCTGTTGCa                                >  1:1401499/1‑150 (MQ=255)
                                                                                                                       aaGGCGGCATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTGTATTTGCTGTTGCa                                >  2:4316151/1‑150 (MQ=255)
                                                                                                                       aaGGCGGCATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCa                                <  1:3421402/150‑1 (MQ=255)
                                                                                                                         ggcggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAgg                              >  2:1070420/1‑150 (MQ=255)
                                                                                                                           cggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCg                            <  1:902095/150‑1 (MQ=255)
                                                                                                                           cggcATCATCACGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCg                            <  1:1654573/150‑1 (MQ=255)
                                                                                                                            ggcATCATCACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCgg                           <  1:2116657/150‑1 (MQ=255)
                                                                                                                               atcatcACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGaa                        >  2:3527155/1‑150 (MQ=255)
                                                                                                                               atcatcACGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGaa                        >  2:4715872/1‑150 (MQ=255)
                                                                                                                                tcatcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGaaa                       >  2:2877940/1‑150 (MQ=255)
                                                                                                                                 catcaCGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAt                      >  1:1090726/1‑150 (MQ=255)
                                                                                                                                  atcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAtt                     <  1:655145/150‑1 (MQ=255)
                                                                                                                                   tcaCGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAttt                    >  1:2377961/1‑150 (MQ=255)
                                                                                                                                   tcaCGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAAttt                    <  1:2896226/150‑1 (MQ=255)
                                                                                                                                    caCGTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttc                   >  2:2369735/1‑150 (MQ=255)
                                                                                                                                      cGTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttctt                 <  1:4252598/150‑1 (MQ=255)
                                                                                                                                       gTCGCGGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttcttc                >  2:4642127/1‑150 (MQ=255)
                                                                                                                                       gTCGCGGGGCTCGCACCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATttcttc                >  1:2632675/1‑150 (MQ=255)
                                                                                                                                           cgGGGCTCGCTCCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCtgt            >  2:2092858/1‑150 (MQ=255)
                                                                                                                                              ggCTCGCTCCCACAGCAGCGTGTCCACTGCGTAAGTCGGGCGCTTAACTGATCGTAATCCCAGACTTCCACCTCACTTATTCCTTCCATACAGTAGTCATCGTCCCCCTGTTTTATTGTCTGTTGCAGGCGGGAAATTTCTTCCTGTGct         <  2:2310411/150‑1 (MQ=255)
                                                                                                                                                ctcgctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGctct       >  1:2258510/1‑150 (MQ=255)
                                                                                                                                                   gctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctggc    <  1:3122284/150‑1 (MQ=255)
                                                                                                                                                   gctcCCAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTctggc    >  1:4273436/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: