Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,4090CA41.3% 166.7 / 39.3 93S13R (AGC→AGArusADLP12 prophage; endonuclease RUS
Reads supporting (aligned to +/- strand):  ref base C (48/6);  new base A (1/37);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 2.91e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCAC  >  NC_000913/573259‑573556
                                                                                                                                                           |                                                                                                                                                     
acgaAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGc                                                                                                                                                       <  4:1097798/151‑1 (MQ=255)
 tgaAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGc                                                                                                                                                       >  1:26319/2‑150 (MQ=255)
    aaTCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATAACTACAGCATCACATTACCCTGGCCTCCGAGCaat                                                                                                                                                    >  2:1508955/1‑150 (MQ=255)
     aTCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCaata                                                                                                                                                   >  1:202576/1‑150 (MQ=255)
      tCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCaataat                                                                                                                                                 >  3:423974/1‑151 (MQ=255)
        gACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATc                                                                                                                                                >  1:415469/1‑150 (MQ=255)
        gACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCg                                                                                                                                               >  2:50590/1‑151 (MQ=255)
        gACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCg                                                                                                                                               >  1:249255/1‑151 (MQ=255)
          ccgccgCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCt                                                                                                                                             >  2:1231377/1‑151 (MQ=255)
           cgccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCt                                                                                                                                             >  2:890964/1‑150 (MQ=255)
            gccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTAt                                                                                                                                          >  1:222306/1‑151 (MQ=255)
             ccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTAtt                                                                                                                                         <  1:144887/151‑1 (MQ=255)
             ccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTAtt                                                                                                                                         >  4:708301/1‑151 (MQ=255)
                 acgcaTTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCg                                                                                                                                     >  4:1585237/1‑151 (MQ=255)
                    caTTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGcc                                                                                                                                   >  2:1224058/1‑150 (MQ=255)
                         tgttgACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAAt                                                                                                                              >  2:1392741/1‑150 (MQ=255)
                           ttgACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATcg                                                                                                                            >  3:44141/1‑150 (MQ=255)
                             gACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCgg                                                                                                                         >  2:1552012/1‑151 (MQ=255)
                              aCGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCggg                                                                                                                        >  2:840080/1‑151 (MQ=255)
                                gtggaatataaaaaaGAATGCGCGCTGGAAGGTAGGGAAAAAACACAGGTTTTCGGGCTGAAAGGGGGGTTTATTAGGGCGTGAATACCTACAGCACCACATTCCCCGGGCCCCGGAGAAATAATCTCTTTTCCCGCCATAATCGCgg                                                                                                                         <  2:1163713/137‑1 (MQ=255)
                                gCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGgcg                                                                                                                      >  3:164430/1‑151 (MQ=255)
                                   ggATATAAAAAAGAAGGCGGGTGGGAAGTTAGGGGGAAAACACAGGTTTTTGGGTTAAAGGGGGGGTTTATTAGGGCGTAAATCCCTACAGCATCACATTCCCCGGCCCCCCGAGAAATAACCGCTATTACCGCCATAATCGCGGgcgc                                                                                                                     <  4:758437/149‑1 (MQ=255)
                                   gattaaaaaaaaGAAGGCGCGGGAAAAGGTAGGGCGAAAACACGGGTTATCGGGAAGAAAGAGGGGTTTATAAAGGCGTAAATCCCACCACCCCCCATTCCCCCGGCCCCCCGACCAATAACCCCTATTACCGCCATAATCGCgg                                                                                                                         <  2:107098/137‑1 (MQ=255)
                                    gATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGcgcac                                                                                                                  >  3:269080/1‑150 (MQ=255)
                                        taaaaaaGAAGGCGCGCGGGAGGGTATGGGGAAAACACGGTTTATCTGGCTAAAAGGGGGGTTTATTAAGGGGTGAATCCCTACAGCATCACATTCCCCTGGCCCCCGAGAAATAATCGTTTTTACCGCCATAATCGCGGGcgcacgc                                                                                                                <  3:222321/145‑1 (MQ=255)
                                        tGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGcacgcacg                                                                                                             >  1:415152/1‑151 (MQ=255)
                                        tGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGcacgcacg                                                                                                             >  4:310755/1‑151 (MQ=255)
                                         gCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGt                                                                                                            >  4:1437146/1‑151 (MQ=255)
                                          aaaaaaaAGGCGCGCGGGAAGGTGGGGTGAAAACACTGTTTTTCGGGTGGAAAGGGGGGTTTTTTAAGGAGGAAATACTTCCAGCACCACATTCCCCTGTCCCCAAAGAAATACACGTTATTCCCGCCATAATCGCGGGcgcacgc                                                                                                                <  2:1310678/145‑1 (MQ=255)
                                          aaaaaGAAGGCGGGCGGGAAGGTAGGGGGAGAACCCAGTTTATCTGGCGAAAAGGGGGGTTTATTAGGGCGTAAATCCCTACAGCATCACATCCCCCTGCCCCCGGAGAAATATCCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTc                                                                                                           <  3:97821/150‑1 (MQ=255)
                                           aaaaGAATGCGCGTTGGAAGGTATGGCGAGAACACGGTTTATGTGGCTGAAGGAGGGGGTTATTAAGGCGTGAATACCTACACCATCACATTACCCTGGCCTCCGACCAATAATCGCTATTACCGCCATAATCGCGGGcgcacg                                                                                                                 <  2:1573165/144‑1 (MQ=255)
                                           aaaaGAATGCGCGGGGAAAGGTAGGGGGAAAACACAGTTTATCTGGCTGAAGGGGGGGTTTATAAGGGCGTAAATCCCTCCACCACCACATCCCCCTGCCCTCCGAGCAAAAATCTTTATTACCGCCATAATCGCGGGCGCACGCACGTCa                                                                                                          <  1:442836/151‑1 (MQ=255)
                                           aaaaGAAGGCGGGCGGAAAGGAAGGGAAAAAACCCAGTTTATCTGGGAAAAGGGGGGGTTTTTGAAGGGGTAAATCCCCACAGCATCACATTCCCCGGCCCACGGGAAAATAATCGATATTGCCGCCATAATCGCGGGcgcacgca                                                                                                               <  2:840117/146‑1 (MQ=255)
                                            aaaGAATGCGCGCTGGAAGGTAGGGCGAGAACCGAGGTTGTCGGGAAGAAAGGGGGGTTTATTAGGGCGTTAATCCCTTCAGATTATCCTCCCCCCGCCGCCAGAGAAAAAATGGCTTTCACCGCCATAATCGCGGGCGCACGCACGTCAg                                                                                                         <  2:1018410/151‑1 (MQ=255)
                                                aaTGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAgcgc                                                                                                      >  2:1447652/1‑150 (MQ=255)
                                                        cgggAAGGTAGGGCGAAAACCGAGGTTATGGGGAAAAAAGGGGGGTTTATAAGGGAGAAAATTCTTACAGCAACTCTTCCCCTTGCCCTCCGAGAAATAATCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGggc                                                                                              <  3:202590/148‑1 (MQ=255)
                                                           gAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGggcaggc                                                                                          >  3:138716/1‑151 (MQ=255)
                                                             aGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCAt                                                                                        >  4:1647265/1‑151 (MQ=255)
                                                              ggTATGGCGAAAACACGGTTTATCTGGCTGAAAGGGGGGGTTATAAAGGCGTAAATCCCTCCACCATCACATTCCCCTGCCCTCAGACAAAAAACCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGggcaggc                                                                                          <  4:1593113/148‑1 (MQ=255)
                                                                tATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGTGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATAcc                                                                                     >  1:469820/1‑151 (MQ=255)
                                                                tATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATAcc                                                                                     >  1:272382/1‑151 (MQ=255)
                                                                   ggggagaACACAGTTTATGTGGCTGAAAGGGGGGTTTATTAAGGGGGAAATCCCTACAGCACCACTTCCCCCGGCCCTCCGACCAATAATCGTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATa                                                                                       <  2:1512090/143‑1 (MQ=255)
                                                                   ggagaaaccacaGTTTATCGGGTGGAAAGGGGGGTTTATTAAGGCGTAAATCCCTACAGCACCACATTCCCCGGCCCCCCGAAAAAAAATCCTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATAc                                                                                      <  4:630366/139‑1 (MQ=255)
                                                                     cGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGAt                                                                                >  2:720328/1‑151 (MQ=255)
                                                                          acacAGGTTATCGGGCAAAAGGGGGGGGTTATAAGGGCGTAAACCCCCACACCACCACATCCCCCGGCCCCCCGAGAAAAAACCTCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAAc                                                                             <  2:150406/149‑1 (MQ=255)
                                                                            acaGGTTATCTGTATGAAAGGGGGTTTTATAAAGGCGTGAATCCCCACAACATCCCTTTCCCCTGCCCCCGGAGAAATACTCGCTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATa                                                                               <  2:1399356/145‑1 (MQ=255)
                                                                             caGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGc                                                                        >  2:1162594/1‑151 (MQ=255)
                                                                                 ttctgcggCTGAAAGGGGGTGTTTATAAGGTGTGAATCCCGAATCTATCACTTCCCCGTCTCCACAAAGCAATAAAAGGTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATCCCGCGATAACGTc                                                                          <  1:144578/139‑1 (MQ=255)
                                                                                 ttatttggggaaaaGGGGGGGTTAATAGGCGGTAAATCCCCACACCACCACTTCCCCGTGGCCCCCGACAAATATGCCCTATCGCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGtt                                                                          <  2:1595049/134‑2 (MQ=255)
                                                                                 ttATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGccc                                                                      >  2:1620642/1‑149 (MQ=255)
                                                                                 ttATCTGGCTAAAAGAGGGGGTTATTAGGGCGTAAATCCCTACAGCATCACATTCCCCTGGCCCCGGAGAAATAATCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGccc                                                                      <  3:26321/149‑1 (MQ=255)
                                                                                  tATCTGGCTGATAGCGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGTGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGa                                                                    >  2:377838/1‑150 (MQ=255)
                                                                                   aTCGGGCTAAAAGGGGGGTTTTTAAGGGCGAGAAACCCCACAGCACCACTTCCCCCGGGCCCCGGACAAATAGCCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAAt                                                                  <  2:950433/151‑1 (MQ=255)
                                                                                       gtgtGAAAGAGGGGTTTATAAGGGCGTAATTCCCCACAGCAACCCATCCCCCTGCCCCCGGAGAAAAAACATTTATTTCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATc                                                                 <  4:1280587/145‑1 (MQ=255)
                                                                                       ggCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCAtt                                                              >  2:945237/1‑151 (MQ=255)
                                                                                         cTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTa                                                             >  4:1015029/1‑150 (MQ=255)
                                                                                          tGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTaa                                                            >  2:768196/1‑150 (MQ=255)
                                                                                            aaaGGGGGGTTTATTAGGGCGTAAATCCCTACAACATCACATTCCCCTGCCCCCCGACAAATAATCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTa                                                             <  3:272398/147‑1 (MQ=255)
                                                                                            aaaGGGGGGGTTATTAAGCCGTAAAACCCTACAACACCCCATTCCCCTGGCCTCCGAGAAAAAATCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTaa                                                            <  3:379741/148‑1 (MQ=255)
                                                                                               ggggggtTTATAAAGGCGTAAATCCCTACAGCATCCCATTCCCCTGCCCCCCGAGAAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                       <  2:1009545/148‑1 (MQ=255)
                                                                                               gagGGGTTTATTAAGGCGTAAATCCCAACACCATCACATTCCCCTGCCCCCCGACAAATAACCCTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                       <  4:524104/150‑1 (MQ=255)
                                                                                                 ggggtTTATTAGGGCGTAAATCCCAACAGCATCACATTCCCCGGGCCCCCGAAAAAAAACCGCTATTCCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                     <  1:228356/150‑1 (MQ=255)
                                                                                                 ggggtTTATTAAGGCGAAAACCCCTACACCACCACTTCCCCCTGCCCCCCGAGAAATAACCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCaa                                                    <  2:1107810/151‑1 (MQ=255)
                                                                                                  gggtTTATAAGGGGGTAAATCCCTACAGAACCACATCCCCCGGCCCCCCGAGAAAAAAGCCTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAAt                                                   <  2:707351/151‑1 (MQ=255)
                                                                                                  gggtTTATAAAGGGGTAAACCCCACCAGCATCACTTCCCCCGGGCCTCGGAGAAATAACCGTTATTCCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAAt                                                   <  4:1128560/151‑1 (MQ=255)
                                                                                                       ttatgaAGGGGTGATGCCCAACACCATCACATCCCCCTGCCCCCCGAATAATAGTCGCTCTTGCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATTCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATg                                                  <  2:1448877/148‑1 (MQ=255)
                                                                                                       ttataaGGGCGTGAATCCCTACACCATCACATTCCCCTGCCCCCCGAGAAATAATCGTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                <  3:279176/150‑1 (MQ=255)
                                                                                                         attaGGGGGTAAATCCCTACACCACCACTTCCCCCGGCCCCCGGAAAAATAATCGTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              <  2:35249/150‑1 (MQ=255)
                                                                                                         attaAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              >  1:17159/1‑150 (MQ=255)
                                                                                                           taAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGat                                            >  2:1526238/1‑150 (MQ=255)
                                                                                                             aaGGCGTAAATACCTACACCATCACATCCCCCTGCCCCCCGAGAAATAATCGTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              <  3:146158/147‑1 (MQ=255)
                                                                                                              aggggtaaaaTCCCTAAAGCATCACATCCCCCGGCCCCCGGAAAAATACTCCTTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              <  2:1583283/139‑1 (MQ=255)
                                                                                                                 gCGTAAATACCTACAGCACCACTTTCCCCTGGCCTCGGAGAAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCgg                                       <  3:245809/151‑1 (MQ=255)
                                                                                                                  cGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                      >  3:434015/1‑151 (MQ=255)
                                                                                                                   gaaattcccTACAGCACCACTTTCCCCGGCCCCCCGAGAAATAACCGTTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGcc                                     <  4:1372509/144‑1 (MQ=255)
                                                                                                                   gTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGcc                                     >  2:441700/1‑151 (MQ=255)
                                                                                                                    tGAATCCATACAGCATCACATTCCCCGGCCCCCCGAGAAATAATCCCTATACCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                      <  4:18038/149‑1 (MQ=255)
                                                                                                                    tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCt                                    >  4:973257/1‑151 (MQ=255)
                                                                                                                      aacccccacacCACCACACTCCCCTCCCCCCCGACAAATAATCACTATACCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCg                                        <  2:1165313/138‑1 (MQ=255)
                                                                                                                      aaTCCCTACAGCACCACATCCCCCTGCCCCCCGAGAAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTgg                                  <  4:1520112/151‑1 (MQ=255)
                                                                                                                       aTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGc                                 >  3:42074/1‑151 (MQ=255)
                                                                                                                        tACCTACAGCATCACATTACCCTGCCTCCGGAGAATTAATGCTATTTACCGCATTATTCGGGGCGGACCGCACGTCAGCGCGGGGGGGCAGGAAACCCGCGTTAAGGCCGCCCGAATCATTAAAAAGGCAAGGCTGGAAATCGGCTTg                                   >  2:1114021/1‑148 (MQ=255)
                                                                                                                          ccaacacCATCACCTCCCCCGCGCCCCCGAGAAATAATCGTTTTTCCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTAt                              <  2:1144449/148‑1 (MQ=255)
                                                                                                                              cAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCt                          >  2:87910/1‑151 (MQ=255)
                                                                                                                               aGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtg                         >  2:172004/1‑151 (MQ=255)
                                                                                                                                ccATCACATTCCCCGGCCCTCCAAGAAATAATCGCTTTTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgt                        <  2:1435348/150‑1 (MQ=255)
                                                                                                                                   tCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGa                      >  2:1508844/1‑150 (MQ=255)
                                                                                                                                      caTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaaaa                   >  3:246373/1‑150 (MQ=255)
                                                                                                                                       aTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATc                 >  2:520144/1‑151 (MQ=255)
                                                                                                                                        tcccccTGGCCCCCGAGAAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTATGaaaa                   <  4:122143/145‑1 (MQ=255)
                                                                                                                                           ccccggccccccGAGAAATAATCTTTATAACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGc              <  4:1449462/141‑1 (MQ=255)
                                                                                                                                            ccTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCAt            >  4:78209/1‑151 (MQ=255)
                                                                                                                                                  ccTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGt       >  4:1126309/1‑151 (MQ=255)
                                                                                                                                                   cTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGt       >  2:1023288/1‑150 (MQ=255)
                                                                                                                                                         aGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcac  >  1:141310/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                                     
ACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCAC  >  NC_000913/573259‑573556

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: