Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 292,302 C→G 16.2% intergenic (+71/+20) yagJ → / ← yagK CP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,3020CG16.2% 403.2 / 70.6 164intergenic (+71/+20)yagJ/yagKCP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base C (61/74);  new base G (12/14);  total (74/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-02

Only 100 of 166 total aligned reads displayed.

GTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTC  >  NC_000913/292203‑292399
                                                                                                   |                                                                                                 
gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                   <  2:2735958/100‑1 (MQ=255)
  agCGGACTTGCTGATTACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTg                                                                                                 <  1:2821577/100‑1 (MQ=255)
  agCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAccg                                                                                                 <  2:1572693/100‑3 (MQ=255)
  agCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTg                                                                                                 >  2:187682/1‑100 (MQ=255)
    cGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTgag                                                                                               <  1:1314127/100‑1 (MQ=255)
     ggACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTgaga                                                                                              <  2:998510/100‑1 (MQ=255)
        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                           >  2:849607/1‑100 (MQ=255)
        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                           >  2:1841566/1‑100 (MQ=255)
         ttGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACTAATCTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                          <  1:3255296/100‑1 (MQ=255)
         ttGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                          <  1:2035665/100‑1 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                         <  1:2493443/100‑1 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                         >  2:390604/1‑100 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                         <  1:1758070/100‑1 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                         >  1:3816831/1‑100 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGACCTTCTCACTGAGAAGGc                                                                                          >  1:890302/1‑99 (MQ=255)
           gCTGATTATTATCAGTAATCTATATAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                        <  1:1718455/100‑2 (MQ=255)
           gCTGATGACTATCAGTAATCTATATAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                        <  1:486784/100‑2 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                         >  2:1281048/1‑99 (MQ=255)
            cTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTTGATTTTAGTAATTTTTTGAAACAATACAAATAAGAAATTCACACAAAGAAGGCAAt                                                                                       >  1:2692844/1‑100 (MQ=255)
             tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcctta                                                                                      <  1:2189009/100‑4 (MQ=255)
              gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAACTAAGGCCTTCTCACTGAGAAGGCCATTa                                                                                     >  2:3264004/1‑100 (MQ=255)
              gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTa                                                                                     >  2:739703/1‑100 (MQ=255)
                ttaCTATCAGTAATCTATAGAATATTGCGTAATGAAACTGCATTTTAGTAATTTTATGAAACAATACAAATATGGCCTTCTCAGTGAGAAGGccttat                                                                                     <  1:918497/96‑5 (MQ=37)
                 gaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGccttatttg                                                                                  <  2:1850689/100‑8 (MQ=255)
                 gaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCg                                                                                  >  2:2854060/1‑100 (MQ=255)
                  aCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGa                                                                                 >  2:193098/1‑100 (MQ=255)
                   cTATCAGTAATCTATAGAAAATTGCGTAAAGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCTTTCTCACTAAGAAGGCCATTAACGAc                                                                                >  2:2426949/1‑100 (MQ=255)
                   cTATCAGTAATCTATAGAAAATTGCGGAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGccttatttgta                                                                                <  2:1534501/100‑10 (MQ=38)
                   cTATCAGCAATTTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAACTAAGGCCTTCTCACTGAGAAGGCCATTAcc                                                                                   <  2:2575604/97‑1 (MQ=255)
                    tATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACt                                                                               <  2:2616897/100‑1 (MQ=255)
                     aTCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                              <  2:3220002/100‑1 (MQ=255)
                     aTCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                              <  2:3219014/100‑1 (MQ=255)
                     aTCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                              <  1:3106237/100‑1 (MQ=255)
                       cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAc                                                                            >  1:2863829/1‑100 (MQ=255)
                       cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAc                                                                            <  1:2363294/100‑1 (MQ=255)
                        aGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACa                                                                           <  1:2663629/100‑1 (MQ=255)
                         gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAg                                                                          <  2:1254955/100‑1 (MQ=255)
                         gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAg                                                                          <  1:2095575/100‑1 (MQ=255)
                             tCTATAGAAAATTGCGTAACTAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGTCCTTCTCTCTGAGAAGGCCATTACCGACTTACAGtgc                                                                       <  2:1804418/99‑3 (MQ=255)
                                ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGAtt                                                                   >  1:2584786/1‑100 (MQ=255)
                                ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGAtt                                                                   >  1:666560/1‑100 (MQ=255)
                                  aGAAAATTGCGTAACTAAACTGAATTTTAGTAATTTTTTGAAACAATAAAAATAAGGCCTTATCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTg                                                                 >  2:3806883/1‑100 (MQ=255)
                                  aGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTg                                                                 <  2:511041/100‑1 (MQ=255)
                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAATTCGATTTgg                                                                >  1:2735213/1‑100 (MQ=255)
                                    acaaTTGCGTAACGAAACCGAATTATAGTAAATCTTTGAAACAATACAAATAAGTCCTCCTCACCAACAAAGCCATTACCGACTTACCGATCGATTTgac                                                               >  2:3496802/3‑98 (MQ=255)
                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGc                                                               <  1:3653619/100‑1 (MQ=255)
                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGc                                                               <  1:2155690/100‑1 (MQ=255)
                                           cGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGcc                                                        >  2:485198/1‑100 (MQ=255)
                                             tAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCAATTACCGACTTACAATTCGATTTGGCTACAGCCGa                                                      >  1:3274584/1‑100 (MQ=255)
                                                    aCTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAAAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACg                                               >  1:1408854/1‑100 (MQ=255)
                                                    aCTGCATTTTAGTAATTTATTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTAACAGTTCGATTTGGCTACAGCCGAAATTACg                                               >  2:3710960/1‑100 (MQ=255)
                                                      tGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGt                                             >  2:262375/1‑100 (MQ=255)
                                                      tGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGt                                             >  1:2131467/1‑100 (MQ=255)
                                                        cATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtct                                           >  2:2610961/1‑100 (MQ=255)
                                                        cATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtct                                           >  2:3380470/1‑100 (MQ=255)
                                                           tttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCt                                        <  2:2250963/100‑1 (MQ=255)
                                                           tttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCGCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCt                                        <  2:2632422/100‑1 (MQ=255)
                                                            ttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCtt                                       >  2:166081/1‑100 (MQ=255)
                                                              aGTAATTTTTTGAATTAATATAAATAAGGTCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                     <  1:1801176/100‑1 (MQ=255)
                                                              aGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCTATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                     <  2:577632/100‑1 (MQ=255)
                                                               gTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAc                                    >  2:3650196/1‑100 (MQ=255)
                                                               gTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAc                                    <  2:266101/100‑1 (MQ=255)
                                                                tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGGAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAcc                                   <  2:2011374/100‑1 (MQ=255)
                                                                tAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAcc                                   <  2:540188/100‑1 (MQ=255)
                                                                     ttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAaca                              <  1:1412212/100‑1 (MQ=255)
                                                                      tttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGTCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAacac                             >  2:3292606/1‑100 (MQ=255)
                                                                        tGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGCCTCCTTCACCAAACACtt                           >  2:2329407/1‑100 (MQ=255)
                                                                           aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGc                        >  2:3202461/1‑100 (MQ=255)
                                                                           aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGc                        <  1:2415317/100‑1 (MQ=255)
                                                                             cAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCt                       <  2:2714692/99‑1 (MQ=255)
                                                                              aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACTGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCttt                     <  1:1012618/100‑1 (MQ=255)
                                                                              aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAACCACTTTGCatt                     >  2:1207006/1‑100 (MQ=255)
                                                                              aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCttt                     <  1:391861/100‑1 (MQ=255)
                                                                               aTACATATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTGACGGTCTCCTTCACCAAACACTTTGCTTTc                    <  1:122134/100‑1 (MQ=255)
                                                                               aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTc                    <  1:464324/100‑1 (MQ=255)
                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCa                   <  2:2817184/100‑1 (MQ=255)
                                                                                tACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCa                   <  2:299724/100‑1 (MQ=255)
                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCCTTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                 >  1:1227344/1‑100 (MQ=255)
                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                 >  2:1228476/1‑100 (MQ=255)
                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACATACGAAGTTACAATCTCCTTCACCAAACACTTTGCTTTCAACtt               >  1:3527838/1‑100 (MQ=255)
                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACtt               >  2:2192211/1‑100 (MQ=255)
                                                                                    aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCAAAAATACGATCTCCTTCACCAAACACTTTGCTTTCAACtt               >  1:1237293/1‑100 (MQ=255)
                                                                                     aTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACttt              <  1:660496/100‑1 (MQ=255)
                                                                                      tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg             <  1:2314842/100‑1 (MQ=255)
                                                                                      tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg             <  1:397547/100‑1 (MQ=255)
                                                                                      tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg             <  1:1952898/100‑1 (MQ=255)
                                                                                       atggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCAAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg            >  1:354640/3‑100 (MQ=255)
                                                                                       atggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTAGGCAAAAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg            >  1:2947135/3‑100 (MQ=255)
                                                                                       aaGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg            <  1:187687/100‑1 (MQ=255)
                                                                                        tggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGt           >  1:826656/2‑100 (MQ=255)
                                                                                         ggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTc          >  1:3012404/1‑100 (MQ=255)
                                                                                         ggCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTc          <  1:2071173/100‑1 (MQ=255)
                                                                                          gCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCa         >  1:2580900/1‑100 (MQ=255)
                                                                                           ccTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAg        <  2:1322118/100‑1 (MQ=255)
                                                                                            cTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAgg       <  1:3032026/100‑1 (MQ=255)
                                                                                             ttCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGt      >  1:598363/1‑100 (MQ=255)
                                                                                               ctcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAg    <  1:1562068/100‑1 (MQ=255)
                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGt   >  2:3634670/1‑100 (MQ=255)
                                                                                                 cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTc  >  1:2517552/1‑100 (MQ=255)
                                                                                                 cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTc  >  1:3256024/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTC  >  NC_000913/292203‑292399

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