Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 292,302 C→G 22.3% 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,3020CG22.3% 470.9 / 142.7 219intergenic (+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 (85/84);  new base G (17/32);  total (103/116)
Fisher's exact test for biased strand distribution p-value = 7.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

Only 100 of 227 total aligned reads displayed.

GTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCA  >  NC_000913/292203‑292401
                                                                                                   |                                                                                                   
gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                     <  2:56717/100‑1 (MQ=255)
gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                     <  2:2464089/100‑1 (MQ=255)
gtagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                     <  1:513462/100‑1 (MQ=255)
 tagCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCCTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAc                                                                                                     <  1:2180384/99‑1 (MQ=255)
  agCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTg                                                                                                   <  2:924598/100‑1 (MQ=255)
     ggACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTgaga                                                                                                <  2:19340/100‑1 (MQ=255)
       aCTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                              <  2:1892125/100‑1 (MQ=255)
        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                             <  2:3144142/100‑1 (MQ=255)
        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                             <  2:403063/100‑1 (MQ=255)
        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAgg                                                                                             <  2:4051900/100‑1 (MQ=255)
         ttGCTGATGACTATCAGTTATCTATAGAAACTTGCGTAACGAAACTGCATTTTAGTAATTTGTTGAAACAATACAAATAAGGCCTTCTCAGGGAGAAgg                                                                                             <  1:1492876/99‑1 (MQ=255)
         ttGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                            <  2:249651/100‑1 (MQ=255)
         ttGCTGATGACTATCAGGAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAggn                                                                                            <  2:1656032/100‑2 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                           <  2:2897781/100‑1 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                           <  1:232133/100‑1 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                           >  2:3392462/1‑100 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGcc                                                                                           >  2:3682134/1‑100 (MQ=255)
          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAAGTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                           <  1:2092905/100‑1 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAATATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                          <  1:3828872/100‑2 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                          <  2:2791333/100‑2 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                          <  1:247422/100‑2 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCa                                                                                          >  1:1415447/1‑100 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCa                                                                                          <  1:1303467/100‑1 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCa                                                                                          >  1:4271393/1‑100 (MQ=255)
           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCa                                                                                          >  1:4035118/1‑100 (MQ=255)
           gCTGATAACTATCAGTAATCTATATCAACTTGCGCAACGCAACTGTATTTTAGTAATTGTTTGAAACAATACACATAAGGCCTTCTCAGTGAGAACGcct                                                                                          <  1:3486156/100‑2 (MQ=255)
            cTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCAt                                                                                         >  1:1561273/1‑100 (MQ=255)
            cTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCAt                                                                                         >  2:4236919/1‑100 (MQ=255)
            cTGATGACACTCACTTATATATTGTAAATCGCCTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAATTCCTTCTCAGTTAGAAGGcc                                                                                           <  1:3668686/98‑1 (MQ=255)
             tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcctta                                                                                        <  1:2958145/100‑4 (MQ=255)
             tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcctta                                                                                        <  2:800089/100‑4 (MQ=255)
             tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcctta                                                                                        <  2:2057454/100‑4 (MQ=255)
             tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCAtt                                                                                        <  2:3792470/100‑1 (MQ=255)
             tgatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCAtt                                                                                        >  1:934857/1‑100 (MQ=255)
              gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTa                                                                                       <  2:634251/100‑1 (MQ=255)
                tgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTAcc                                                                                     <  2:2664777/100‑1 (MQ=255)
                tgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTAcc                                                                                     <  1:2020070/100‑1 (MQ=255)
                 gaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCg                                                                                    >  2:3214480/1‑100 (MQ=255)
                 gaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCg                                                                                    <  2:1282353/100‑1 (MQ=255)
                   cTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGAc                                                                                  <  2:3049648/100‑1 (MQ=255)
                   cTATCAGTAATCTATAGAAAATTGCGTAAAGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGGCTTCTCACTGAGAAGGCCATTACCAAc                                                                                  >  1:1118209/1‑100 (MQ=255)
                    tgtCAGTAATCTAGAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACt                                                                                 >  2:2979680/3‑100 (MQ=255)
                    tATCAATAATCTATAGAAAATTGCATAACCAAACTGCATTTTAGTAATTTTTTTAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTAcc                                                                                     <  1:3852465/96‑1 (MQ=255)
                       cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTTAATCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAc                                                                              >  1:2777228/1‑100 (MQ=255)
                       cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAc                                                                              >  2:1070466/1‑100 (MQ=255)
                         gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAg                                                                            >  1:821671/1‑100 (MQ=255)
                           aaTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGtt                                                                          >  2:48112/1‑100 (MQ=255)
                            aTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTc                                                                         <  1:1963422/100‑1 (MQ=255)
                            aTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTc                                                                         <  1:1757135/100‑1 (MQ=255)
                            aTCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTc                                                                         <  2:1540868/100‑1 (MQ=255)
                              cTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGa                                                                       <  1:3309468/100‑1 (MQ=255)
                                ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACACTTCGAtt                                                                     >  1:3581310/1‑100 (MQ=255)
                                 taGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGCAGGCCATTACCGACTTACAGTTCGAt                                                                      <  1:3704700/98‑1 (MQ=255)
                                  aGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTg                                                                   >  1:393348/1‑100 (MQ=255)
                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                  >  2:2163425/1‑100 (MQ=255)
                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                  >  1:3055549/1‑100 (MQ=255)
                                     aaaTTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCt                                                                >  2:167753/1‑100 (MQ=255)
                                        ttGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACa                                                             <  2:1603637/100‑1 (MQ=255)
                                           cGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGcc                                                          >  1:1609126/1‑100 (MQ=255)
                                              aaCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAAAGTTCGATTTGGCTACAGCCGaa                                                       >  1:1975925/1‑100 (MQ=255)
                                               aCGAAACTGCTTTTTAGTAATTTTTTGATATAATACAAATAATGCCTTCTTAATGATAAGGCCATTATCGACCTACAGTTGCATTAGGCTACAGCCGAAg                                                      >  1:57189/1‑100 (MQ=255)
                                               aCGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAg                                                      <  2:4315484/100‑1 (MQ=255)
                                                cGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGt                                                     >  2:1442280/1‑100 (MQ=255)
                                                 gAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGtt                                                    >  2:1303433/1‑100 (MQ=255)
                                                  aaaCTTCATTTTAGTAATTTTTTGAAACATTCCAAATAAGCCCCTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAg                                                      <  1:672779/97‑1 (MQ=255)
                                                     cTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACgg                                                <  1:645747/100‑1 (MQ=255)
                                                     cTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACg                                                 >  1:3529626/1‑99 (MQ=255)
                                                       gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTATGGtc                                              <  2:203201/100‑1 (MQ=255)
                                                       gCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtc                                              >  1:1254558/1‑100 (MQ=255)
                                                         aTTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtctc                                            >  1:313298/1‑100 (MQ=255)
                                                          ttttAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTcc                                           >  2:732677/1‑100 (MQ=255)
                                                                   ttttttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCaaa                                  <  2:3492108/100‑1 (MQ=255)
                                                                      tttGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAacac                               >  2:858237/1‑100 (MQ=255)
                                                                        tGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACtt                             >  2:1865555/1‑100 (MQ=255)
                                                                        tGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACtt                             >  1:3216121/1‑100 (MQ=255)
                                                                         gAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACttt                            <  2:398058/100‑1 (MQ=255)
                                                                         gAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACttt                            >  2:3233514/1‑100 (MQ=255)
                                                                          aaaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTg                           >  2:4109024/1‑100 (MQ=255)
                                                                           aaCAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGc                          >  1:4061739/1‑100 (MQ=255)
                                                                               aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTc                      <  2:1239401/100‑1 (MQ=255)
                                                                                  cAAATAAGGCCTTCTGACCGAGAAGGCCATTACCGACGTACAGGTCGGTTCAGCTAAAGGCGAAGTTACGGTCTCCTTCCACAAACACTTTGTTTACa                     >  1:3364628/1‑98 (MQ=255)
                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCaa                    >  1:3335622/1‑99 (MQ=255)
                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                   >  1:2623085/1‑100 (MQ=255)
                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                   >  2:2713542/1‑100 (MQ=255)
                                                                                  cAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                   >  2:2141984/1‑100 (MQ=255)
                                                                                   ggtaatggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACt                  >  2:3326356/7‑100 (MQ=255)
                                                                                      tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg               <  1:1158107/100‑1 (MQ=255)
                                                                                      tAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg               >  1:3871821/1‑100 (MQ=255)
                                                                                       atggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg              >  2:4404076/3‑100 (MQ=255)
                                                                                       atggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg               >  2:198163/3‑99 (MQ=255)
                                                                                        tggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGt             >  1:2712357/2‑100 (MQ=255)
                                                                                         ggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTc            >  2:3770548/1‑100 (MQ=255)
                                                                                         ggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTc            >  2:3551165/1‑100 (MQ=255)
                                                                                           ccTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAg          >  2:2703860/1‑100 (MQ=255)
                                                                                             ttCTCAGTGAGAAGGCCTTTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGt        <  1:3183210/100‑1 (MQ=255)
                                                                                               ctcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAg      <  1:4123933/100‑1 (MQ=255)
                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGt     >  2:292576/1‑100 (MQ=255)
                                                                                                tcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGt     >  2:91132/1‑100 (MQ=255)
                                                                                                   cTGAGAAGGCCATTACCGACTTTCAGTTCGATTTGGCTACAGCCTAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAgg         <  1:4332743/93‑1 (MQ=255)
                                                                                                   cTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCa  >  2:178885/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCA  >  NC_000913/292203‑292401

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