Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,709 T→C 55.3% E11G (GAA→GGA)  ydaC ← DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7090TC55.3% 56.0 / 188.3 152E11G (GAA→GGA) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base T (27/41);  new base C (43/41);  total (70/82)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 166 total aligned reads displayed.

AACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTC  >  NC_000913/1413612‑1413806
                                                                                                 |                                                                                                 
aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCActttcc                                                                                                 >  1:760259/1‑100 (MQ=255)
aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAccttcc                                                                                                 >  1:4104230/1‑100 (MQ=255)
  cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccctt                                                                                               >  1:3449021/1‑100 (MQ=255)
   gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccctct                                                                                              <  2:2878589/100‑3 (MQ=255)
   gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCttt                                                                                              <  2:4262619/100‑1 (MQ=255)
    ccTGTTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGACTTTGAATATCTTGTCCAGTTAGTAGGACTGCCAccttcc                                                                                                 <  1:1473658/96‑1 (MQ=255)
    ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCtttt                                                                                             <  2:1283894/100‑1 (MQ=255)
    ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCtttt                                                                                             <  2:1997051/100‑1 (MQ=255)
    ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCtttt                                                                                             >  2:1122144/1‑100 (MQ=255)
     cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTATATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTc                                                                                            <  1:2375904/100‑1 (MQ=255)
     cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTc                                                                                            <  2:69895/100‑1 (MQ=255)
        cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGCTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAAt                                                                                         <  1:2653305/100‑1 (MQ=255)
        cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAAt                                                                                         >  2:959440/1‑100 (MQ=255)
         ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATa                                                                                        <  2:2145642/100‑1 (MQ=255)
           cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGt                                                                                      <  1:1219068/100‑1 (MQ=255)
            ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTg                                                                                     >  2:3303559/1‑100 (MQ=255)
            ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTg                                                                                     >  2:3646210/1‑100 (MQ=255)
              ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGc                                                                                   <  1:2747470/100‑1 (MQ=255)
              ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAACAGTGGc                                                                                   >  1:535656/1‑100 (MQ=255)
               tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCg                                                                                  >  2:2341063/1‑100 (MQ=255)
               tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGCAGGAGTGCCACCTCCCTTTTCAATAGTGGCg                                                                                  >  2:95428/1‑100 (MQ=255)
                ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCgg                                                                                 <  1:360367/100‑1 (MQ=255)
                 tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGt                                                                                <  2:2647840/100‑1 (MQ=255)
                  aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTa                                                                               <  2:137988/100‑1 (MQ=255)
                     tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAtt                                                                            >  2:1689889/1‑100 (MQ=255)
                     tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAtt                                                                            >  2:3889423/1‑100 (MQ=255)
                       agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAtttc                                                                          <  2:1497935/100‑2 (MQ=255)
                       agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAtttt                                                                          <  2:3410545/100‑1 (MQ=255)
                        gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTta                                                                         <  1:1095209/100‑1 (MQ=255)
                         agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtat                                                                        <  2:4102595/100‑1 (MQ=255)
                         agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtat                                                                        <  2:1751203/100‑1 (MQ=255)
                          gtaTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGTAGTGCCATATCCCGTTTCAATAGTGCCGGTAAtttt                                                                          <  1:2263512/95‑1 (MQ=255)
                           cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                      <  2:3155894/100‑1 (MQ=255)
                           cATTTTTTCGCATTCTGATTTCGTTAGGTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                      <  1:1048040/100‑1 (MQ=255)
                            aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                      <  1:2065059/99‑1 (MQ=255)
                             ttttttCGCATTCTGATTTCGTTAGTTTAGGTTTTGAATATCTTGTCCTGTTAGTATGAGTGGCACCTTCCTTTTCACTAGTGGCGGTAATTTTATACa                                                                     <  1:1681670/99‑1 (MQ=255)
                                 ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                >  2:2140419/1‑100 (MQ=255)
                                 ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                >  2:3696073/1‑100 (MQ=255)
                                     cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAGATGAACACCt                                                            >  1:2812461/1‑100 (MQ=255)
                                        tccgATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCa                                                         <  1:1774008/97‑1 (MQ=255)
                                        tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTcct                                                         <  1:3696552/100‑2 (MQ=255)
                                            aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattat                                                     <  2:712933/100‑1 (MQ=255)
                                             tttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTATGAGTGCCACCTCCCTTTTCAATAGTTGCGGTAATTTTATACATGAACACCTCCattatt                                                    <  1:131533/100‑1 (MQ=255)
                                              ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatta                                                   <  1:1916345/100‑1 (MQ=255)
                                               tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattat                                                  <  1:849254/100‑1 (MQ=255)
                                                   tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTATTCAATAGTGGCGGTAATTTTATACATGAACACGTCCATTATTATTTcc                                              <  1:3011103/100‑1 (MQ=255)
                                                   tAGTTTAGATTTTGAATATATTGTCCAGTTAGTAGGAGTGCCACATCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcc                                              <  1:1418235/100‑1 (MQ=255)
                                                     gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAAATGAACACCTCCattattattac                                               >  1:3849932/1‑95 (MQ=255)
                                                      tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                           <  2:4104252/100‑1 (MQ=255)
                                                      tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                           >  2:276646/1‑100 (MQ=255)
                                                      tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                           <  1:3160529/100‑1 (MQ=255)
                                                       ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTg                                          <  1:3394209/100‑1 (MQ=255)
                                                         aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGt                                        >  2:227706/1‑100 (MQ=255)
                                                            ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCg                                     <  1:311235/100‑1 (MQ=255)
                                                            ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCg                                     >  2:478367/1‑100 (MQ=255)
                                                            ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTCTACATGAACACTTCCATTATTATTTCCATTGGTTCg                                     >  1:1017171/1‑100 (MQ=255)
                                                             tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                    >  2:1094418/1‑100 (MQ=255)
                                                              ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGtt                                   <  1:1438948/100‑1 (MQ=255)
                                                              ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGtt                                   <  1:2329398/100‑1 (MQ=255)
                                                               tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGttt                                  >  2:4018337/1‑100 (MQ=255)
                                                                  atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                               >  2:1769806/1‑100 (MQ=255)
                                                                   tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTc                              <  2:3407144/100‑1 (MQ=255)
                                                                   tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTc                              >  2:621516/1‑100 (MQ=255)
                                                                   tatCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTc                              <  2:2108260/100‑1 (MQ=255)
                                                                    atCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTcc                             <  2:423989/100‑1 (MQ=255)
                                                                     tCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                            <  2:2715008/100‑1 (MQ=255)
                                                                     tCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATTGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                            <  1:1364106/100‑1 (MQ=255)
                                                                     tCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                            >  2:388895/1‑100 (MQ=255)
                                                                       ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATTAACACCTCCATTATTATTTCCAGCGGTTCGTTTATTCCATc                          <  1:474602/100‑1 (MQ=255)
                                                                           ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                      <  1:1574587/100‑1 (MQ=255)
                                                                            cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                     >  1:401500/1‑100 (MQ=255)
                                                                            cAGTTAGTAGGAGTGCCACCTTCCTTTTCAACAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGTTTCGTTTATTCCATCTTTCg                     <  1:391244/100‑1 (MQ=255)
                                                                            cAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                     <  2:969516/100‑1 (MQ=255)
                                                                             aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                    >  1:3726178/1‑100 (MQ=255)
                                                                             aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                    >  2:2942114/1‑100 (MQ=255)
                                                                             aGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                    >  2:3614597/1‑100 (MQ=255)
                                                                              gTTAGTAGTTGTGGAACCTCGCTATTCAATAGTGGAGGTAATTTTACACATTAACACCTCCATTATTATTTCCTGTGGTTCGTTTATTCCAGCTTTCGa                    <  1:847585/99‑1 (MQ=255)
                                                                              gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                   <  1:728739/100‑1 (MQ=255)
                                                                               ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                  <  1:2140414/100‑1 (MQ=255)
                                                                                tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTg                 <  2:143521/100‑1 (MQ=255)
                                                                                   tagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGctt              >  2:3090545/1‑100 (MQ=255)
                                                                                     ggAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttct            >  1:1342951/1‑100 (MQ=255)
                                                                                     ggAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttct            >  2:3475656/1‑100 (MQ=255)
                                                                                     ggAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttct            <  1:97017/100‑1 (MQ=255)
                                                                                     ggAGTGCCACCTCCCTTTTCAAAAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCCTTTATTCCATCTTTCGAGTGcttc             >  1:2477101/1‑99 (MQ=255)
                                                                                      gAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttc             <  1:1354570/98‑1 (MQ=255)
                                                                                       aGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttt          <  2:2493005/100‑1 (MQ=255)
                                                                                          gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt        >  2:325401/1‑99 (MQ=255)
                                                                                          gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc       >  1:3733256/1‑100 (MQ=255)
                                                                                          gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc       <  2:1764699/100‑1 (MQ=255)
                                                                                          gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc       <  2:3174051/100‑1 (MQ=255)
                                                                                          gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc       <  2:2026936/100‑1 (MQ=255)
                                                                                          gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTAGTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc       >  1:2734777/1‑100 (MQ=255)
                                                                                            cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTAATTCCATCTTTCGAGTGCTTCTTTTTCAc     >  2:1897957/1‑100 (MQ=255)
                                                                                            cACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc     <  2:1607169/100‑1 (MQ=255)
                                                                                            cACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc     <  2:92092/100‑1 (MQ=255)
                                                                                              ccttccttTTCAATAGTGGCGGTAATTTTTGACTTGAACGCCGCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa      <  1:2553456/97‑1 (MQ=255)
                                                                                              cctcccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACtt   >  1:3065864/1‑100 (MQ=255)
                                                                                               cttccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTc  <  2:1478011/100‑1 (MQ=255)
                                                                                               cttccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTGTTTCCATCTTTCGAGTGCTTCCTTTTCa      <  1:864662/96‑1 (MQ=255)
                                                                                                 |                                                                                                 
AACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTC  >  NC_000913/1413612‑1413806

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: