Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,581 G→A 22.8% R135W (CGG→TGG)  yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,5810GA22.8% 390.0 / 111.1 207R135W (CGG→TGG) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (69/89);  new base A (21/26);  total (92/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 212 total aligned reads displayed.

CGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTC  >  NC_000913/281482‑281680
                                                                                                   |                                                                                                   
 gtcCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCgg                                                                                                    >  2:2944110/1‑100 (MQ=255)
  tcCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGc                                                                                                   <  2:1217181/100‑1 (MQ=255)
    cagcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGccc                                                                                                 >  1:2272142/1‑100 (MQ=255)
    cagcagGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGccc                                                                                                 >  2:2179065/1‑100 (MQ=255)
     agcagGGTGAGCGGATGGCAGCGTCCACCGCCGAGCGGAAAGTGGCGCTGAAAATCCGTCTGCCAGAGGCGGCTCGGCGCTTCGTGGTCGAACCGGCCCg                                                                                                >  2:2464006/1‑100 (MQ=255)
           gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCgg                                                                                          >  1:376335/1‑100 (MQ=255)
           gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCgg                                                                                          <  1:3105459/100‑1 (MQ=255)
           gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCgg                                                                                          <  1:2043497/100‑1 (MQ=255)
            tGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCggg                                                                                         <  1:3535680/100‑1 (MQ=255)
            tGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCggg                                                                                         <  2:3220433/100‑1 (MQ=255)
            tGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCggg                                                                                         <  1:2878853/100‑1 (MQ=255)
             gAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGa                                                                                        <  1:649423/100‑1 (MQ=255)
               gCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAAt                                                                                      <  2:1907117/100‑1 (MQ=255)
                 ggATGGCAGCGTCCACCGCCAAAGGGGAGGGGGCCCTTAAAATCCAGCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATg                                                                                     <  2:3090325/99‑1 (MQ=255)
                    tGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCg                                                                                 >  1:3945523/1‑100 (MQ=255)
                    tGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCg                                                                                 <  1:1527853/100‑1 (MQ=255)
                    tGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCg                                                                                 >  2:3584856/1‑100 (MQ=255)
                      gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCgg                                                                                >  2:1543432/1‑99 (MQ=255)
                      gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGt                                                                               <  2:1406716/100‑1 (MQ=255)
                      gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGt                                                                               >  1:2739970/1‑100 (MQ=255)
                          cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAg                                                                           <  1:3910510/100‑1 (MQ=255)
                          cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAg                                                                           <  1:1113747/100‑1 (MQ=255)
                           gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAgc                                                                          <  2:901438/100‑1 (MQ=255)
                            tCCACCGCCAAAGGGAAAGTGGCCCTTCAAATCCCTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAgcg                                                                         <  2:1777832/100‑1 (MQ=255)
                             ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAgcgc                                                                        >  1:2241138/1‑100 (MQ=255)
                              cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAgcgc                                                                        >  2:2426727/1‑99 (MQ=255)
                               aCCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGccc                                                                      <  1:3340177/100‑1 (MQ=255)
                                  gacAAAAGGAAAGTGGCCCTTAAACTCCATCTGCCAGGGGCTGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgc                                                                   <  2:2507196/98‑1 (MQ=255)
                                    cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcag                                                                 <  1:1713944/100‑1 (MQ=255)
                                     aaaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcagc                                                                <  2:1313821/100‑1 (MQ=255)
                                      aaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcagca                                                               >  2:569995/1‑100 (MQ=255)
                                       aaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGcagcag                                                              <  1:3614000/100‑1 (MQ=255)
                                       aaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGcagcag                                                              <  2:1188586/100‑1 (MQ=255)
                                        aGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAgg                                                             <  1:3236284/100‑1 (MQ=255)
                                              gTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATgg                                                       <  1:1130701/100‑1 (MQ=255)
                                              gTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATgg                                                       >  1:3130146/1‑100 (MQ=255)
                                              gTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATgg                                                       >  2:232423/1‑100 (MQ=255)
                                                ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCg                                                     >  2:1567882/1‑100 (MQ=255)
                                                 gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCgg                                                    >  1:2626421/1‑100 (MQ=255)
                                                  cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                   >  1:2573044/1‑100 (MQ=255)
                                                  cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                   >  1:3970600/1‑100 (MQ=255)
                                                   ctttAAAGGCCATCTGCCAGAGGCGGTTCGGGGTGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCTCGGCAGCAGGCCATGGCggg                                                   <  2:265413/97‑1 (MQ=255)
                                                   ccTTAAAATCCATCTGCCAGAGGCGGTTCTGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  <  2:3686837/100‑1 (MQ=255)
                                                   ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  <  2:2573029/100‑1 (MQ=255)
                                                   ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  <  2:2390585/100‑1 (MQ=255)
                                                   ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                  >  2:4078474/1‑100 (MQ=255)
                                                    cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                                 <  2:291930/100‑1 (MQ=255)
                                                    cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                                 <  2:1255853/100‑1 (MQ=255)
                                                    cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCCGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                                 <  1:2799349/100‑1 (MQ=255)
                                                      tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCAt                                               >  1:2126060/1‑100 (MQ=255)
                                                       aaaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATc                                              <  1:1005160/100‑1 (MQ=255)
                                                         aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCan                                            <  1:453236/100‑2 (MQ=255)
                                                         aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAg                                            <  1:3864739/100‑1 (MQ=255)
                                                         aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAg                                            >  2:3007119/1‑100 (MQ=255)
                                                         aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAg                                            <  1:2680048/100‑1 (MQ=255)
                                                          aTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                           <  1:4127990/100‑1 (MQ=255)
                                                             cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                        >  1:25223/1‑100 (MQ=255)
                                                             cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                        >  1:2972479/1‑100 (MQ=255)
                                                             cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                        <  1:3816065/100‑1 (MQ=255)
                                                              aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTCGCGGGTGTGCACGGTGGAGGGCCCGGCAGCAGGCCAAGGCGGGCCATCAGGt                                          >  2:43710/1‑97 (MQ=255)
                                                                cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATg                                      <  2:3152403/99‑1 (MQ=255)
                                                                cTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGGTATgg                                     >  2:2042107/1‑100 (MQ=255)
                                                                 tGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                                    >  1:2340587/1‑100 (MQ=255)
                                                                 tGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                                    >  1:1227093/1‑100 (MQ=255)
                                                                  gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCTCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                   >  1:2911785/1‑100 (MQ=255)
                                                                  gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGaa                                   >  1:917007/1‑99 (MQ=255)
                                                                  gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                   >  1:1154174/1‑100 (MQ=255)
                                                                  gCCAGAGGCGGTTCGGCGCGTCGGGTTCGAACCGGCCCGGGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                   <  2:3845278/100‑1 (MQ=255)
                                                                   ccAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACg                                  >  2:2121715/1‑100 (MQ=255)
                                                                     agagGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGGGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGt                                <  1:588237/100‑1 (MQ=255)
                                                                      gagGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGAAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTg                               <  2:350756/100‑1 (MQ=255)
                                                                        ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCTATCAGGTTATGGACGGTGCt                             >  2:3317409/1‑100 (MQ=255)
                                                                        ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                             >  1:780531/1‑100 (MQ=255)
                                                                        ggCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                             >  1:1183117/1‑100 (MQ=255)
                                                                         gCGGTTCGGCGCGTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTg                            <  2:705198/100‑1 (MQ=255)
                                                                            gTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAg                         <  2:464317/100‑1 (MQ=255)
                                                                               cggcgcggcgTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAg                         <  2:1940641/97‑1 (MQ=255)
                                                                                ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTCTGGCCGGTGCTGAAGGCgg                     >  1:764163/1‑100 (MQ=255)
                                                                                ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                     >  2:337940/1‑100 (MQ=255)
                                                                                  cgcgTGGTGTCCGACCCAGCCCGTGGCCGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggg                    <  1:4134885/99‑1 (MQ=255)
                                                                                  cgcgTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGc                   >  2:3729020/1‑100 (MQ=255)
                                                                                   gcgTCGTGTTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCa                  <  1:3488589/100‑1 (MQ=255)
                                                                                       cgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt              >  2:3797899/1‑100 (MQ=255)
                                                                                       cgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt              <  2:3937070/100‑1 (MQ=255)
                                                                                        gtgtTCGCACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg             <  2:1589896/100‑1 (MQ=255)
                                                                                        gtgtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg             <  1:1805052/100‑1 (MQ=255)
                                                                                        gtgtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg             <  1:322547/100‑1 (MQ=255)
                                                                                         tgtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtgt            >  1:141361/1‑100 (MQ=255)
                                                                                          gtTCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg           >  1:1912578/1‑100 (MQ=255)
                                                                                            tCGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc         >  2:312034/1‑100 (MQ=255)
                                                                                             cGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccc        >  1:893290/1‑100 (MQ=255)
                                                                                             cGAACCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc         >  1:1058583/1‑99 (MQ=255)
                                                                                              gAACGGGCCCGTGGCGGGAATGCGCGGTGGGGCGCGCGGCAGCAGGCCATGGGGGGGCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccgc       <  2:954363/100‑3 (MQ=255)
                                                                                              gAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc       >  2:86050/1‑100 (MQ=255)
                                                                                              gAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcccc       <  2:894685/100‑1 (MQ=255)
                                                                                               aaCCAGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGAAGCAGGCCAGGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCt      <  2:1017582/100‑1 (MQ=255)
                                                                                                aCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTg     >  1:1760848/1‑100 (MQ=255)
                                                                                                 ccGGCCCGTGGCGGGAATGCCCGGTGAAGCGCTCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg    <  1:3882096/100‑1 (MQ=255)
                                                                                                 ccGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGTTGTGCCCCTgg    >  1:516299/1‑100 (MQ=255)
                                                                                                 ccGGCCCGGGGCGGGACTGCCCGCTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCAAGGTGTGcccc       <  2:1346649/97‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTC  >  NC_000913/281482‑281680

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: