Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,523,317 G→T 100% pseudogene (1159/2037 nt) rhsE → pseudogene, Rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,523,3170GT100.0% 180.7 / NA 66pseudogene (1159/2037 nt)rhsEpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/42);  total (24/42)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTG  >  CP009273/1523284‑1523463
                                 |                                                                                                                                                  
gAACTGGTGCGGCTGCTGGACAGGCTGGAGGAATAACTCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGtt                                <  1:5079775/150‑1 (MQ=14)
gAACTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGtt                                >  2:5090607/1‑150 (MQ=14)
gAACTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGtt                                >  1:4320073/1‑150 (MQ=14)
gAACTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGtt                                >  1:1696438/1‑150 (MQ=14)
gAACTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGtt                                <  1:2943559/150‑1 (MQ=14)
 aaCTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTc                               <  1:4245943/150‑1 (MQ=17)
 aaCTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTc                               <  1:3707021/150‑1 (MQ=17)
 aaCTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTc                               >  2:1765461/1‑150 (MQ=17)
  aCTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCa                              <  2:349665/150‑1 (MQ=17)
  aCTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCa                              >  2:4882584/1‑150 (MQ=17)
   cTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAt                             >  2:1256406/1‑150 (MQ=17)
   cTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAt                             >  2:1256372/1‑150 (MQ=17)
   cTGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAt                             <  1:2517849/150‑1 (MQ=17)
    tGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtt                            <  1:1192819/150‑1 (MQ=17)
    tGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtt                            <  2:1014176/150‑1 (MQ=17)
    tGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtt                            <  2:1014219/150‑1 (MQ=17)
    tGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtt                            <  1:1192835/150‑1 (MQ=17)
    tGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtt                            <  2:1014336/150‑1 (MQ=17)
    tGGTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtt                            <  2:2915079/150‑1 (MQ=17)
     ggTGCGGCTGCTGGACAGGCTGGTGGAATAACTGCGGGCTGTGCGCTTGAGCAGTGAAAGCCGGGCGTGGCTTGCTCAGTTCTGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATTCACACCGGCGCGAAAAGTTCAttt                           <  1:5246877/150‑1 (MQ=11)
     ggTGCGGCTGCTGGACAGGCTGGAGGAATATTTCCGGGCAGACCGCGTGAGCAGTGATAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCTGACAGGTGGAGCCGGAATACTCACCGGCGCGAATAGTTCAttt                           <  2:5173186/150‑1 (MQ=11)
     ggTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAttt                           <  2:286035/150‑1 (MQ=18)
     ggTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAttt                           <  1:4396641/150‑1 (MQ=18)
     ggTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAttt                           <  2:5029297/150‑1 (MQ=18)
      gTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtttt                          <  2:11547/150‑1 (MQ=21)
      gTGCGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCAtttt                          >  1:5304610/1‑150 (MQ=21)
        gCGGCTGCTGGACAGGCTGGAGGTATCAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTa                        <  2:4768963/150‑1 (MQ=255)
         cGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACTCACCGGCGCGAACAGTTCATTTTTAc                       <  2:2368059/150‑1 (MQ=255)
         cGGCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTAc                       <  2:1827539/150‑1 (MQ=25)
          ggCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGACTACACACCGGCGCGAAAAGTTCATTTTTAcc                      >  2:4172554/1‑150 (MQ=25)
          ggCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTAcg                      <  1:1749781/150‑2 (MQ=25)
          ggCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTAcc                      <  1:807045/150‑1 (MQ=25)
          ggCTGCTGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTAcc                      <  2:3222630/150‑1 (MQ=25)
           gctgctGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCa                     <  1:5170616/150‑1 (MQ=37)
           gctgctGGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCa                     <  1:2602811/150‑1 (MQ=37)
            ctgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCAc                    <  2:4709960/150‑1 (MQ=37)
            ctgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCAc                    >  2:2525058/1‑150 (MQ=37)
            ctgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCAc                    >  1:2117078/1‑150 (MQ=37)
            ctgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCAc                    >  1:2716178/1‑150 (MQ=37)
             tgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACt                   <  1:1519832/150‑1 (MQ=37)
             tgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACt                   <  2:4519299/150‑1 (MQ=37)
             tgctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACt                   <  2:656164/150‑1 (MQ=37)
              gctgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTg                  <  2:5307743/150‑1 (MQ=37)
                tgGACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCg                >  1:3734853/1‑150 (MQ=37)
                 ggACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGa               >  1:2698334/1‑150 (MQ=37)
                 ggACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGa               <  1:2669556/150‑1 (MQ=37)
                  gACAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGac              >  2:4509284/1‑150 (MQ=37)
                    cAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGacca            >  2:1646476/1‑150 (MQ=255)
                    cAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGacca            <  2:2768196/150‑1 (MQ=37)
                    cAGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGacca            <  2:100055/150‑1 (MQ=37)
                     aGGCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGaccac           <  2:4410747/150‑1 (MQ=37)
                      ggCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTGCATTTTTACCACGGCGaccacc          >  2:3111252/1‑150 (MQ=255)
                      ggCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGaccacc          <  2:84739/150‑1 (MQ=37)
                      ggCTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGaccacc          <  1:2274291/150‑1 (MQ=37)
                        cTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGACTACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgg        >  1:3362453/1‑150 (MQ=255)
                        cTGGAGGAATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgg        >  2:1099706/1‑150 (MQ=37)
                          ggaggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgggg      <  2:2160218/150‑1 (MQ=37)
                          ggaggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgggg      <  1:4807162/150‑1 (MQ=37)
                          ggaggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgggg      >  2:933357/1‑150 (MQ=37)
                          ggaggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTAGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgggg      >  2:1998548/1‑150 (MQ=255)
                          ggaggaATAAATCCGGGCAGACCGCGTGAGCAGAGAAAGCCGGGCGTAGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCgggg      >  2:1998549/1‑150 (MQ=255)
                            aggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGcc    >  2:3531096/1‑150 (MQ=38)
                            aggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGcc    <  1:1075258/150‑1 (MQ=38)
                             ggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCt   <  1:828590/150‑1 (MQ=255)
                             ggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCt   <  1:4838414/150‑1 (MQ=255)
                             ggaATAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCt   >  1:1648434/1‑150 (MQ=255)
                              gaataaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTg  <  2:3827943/150‑1 (MQ=255)
                              gaataaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTg  <  2:2715147/150‑1 (MQ=255)
                              gaataaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTg  >  1:790898/1‑150 (MQ=255)
                                 |                                                                                                                                                  
GAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTG  >  CP009273/1523284‑1523463

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: ‑