Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,413,709 T→C 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,7090TC97.1% 345.4 / ‑1.7 103E11G (GAA→GGA) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (47/53);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 109 total aligned reads displayed.

AACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTT  >  NC_000913/1413612‑1413805
                                                                                                 |                                                                                                
aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccc                                                                                                <  2:2058933/100‑1 (MQ=255)
aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccc                                                                                                >  2:832408/1‑100 (MQ=255)
aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccc                                                                                                <  2:1722185/100‑1 (MQ=255)
aaCGCCTGCTCCTTCTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccc                                                                                                >  2:2100772/1‑100 (MQ=255)
   gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCttt                                                                                             >  1:1957360/1‑100 (MQ=255)
    ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCtttt                                                                                            >  1:2240085/1‑100 (MQ=255)
    ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCtttt                                                                                            >  2:569798/1‑100 (MQ=255)
    ccTGATTCTTTTTTGCCTGAAAGCATTTTTTCGCACTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGGAGCAGTGCCACCTCCCtttt                                                                                            <  2:1938523/100‑1 (MQ=255)
     tctcttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTc                                                                                           <  1:372408/97‑1 (MQ=255)
     cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTc                                                                                           <  1:486289/100‑1 (MQ=255)
      tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCa                                                                                          >  1:2324522/1‑100 (MQ=255)
      tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCa                                                                                          <  2:1758539/100‑1 (MQ=255)
      tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCa                                                                                          <  1:723769/100‑1 (MQ=255)
        cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGCTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAAt                                                                                        >  2:1393469/1‑100 (MQ=255)
        cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAAt                                                                                        >  2:2140242/1‑100 (MQ=255)
         ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATa                                                                                       >  1:616485/1‑100 (MQ=255)
           cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGt                                                                                     >  1:713452/1‑100 (MQ=255)
           cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGt                                                                                     <  2:1374022/100‑1 (MQ=255)
              ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGc                                                                                  <  1:1294176/100‑1 (MQ=255)
               tttATCTGACAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCg                                                                                 <  2:706833/100‑1 (MQ=255)
                 tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGt                                                                               >  1:288954/1‑100 (MQ=255)
                  aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGCTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTa                                                                              <  2:1440176/100‑1 (MQ=255)
                  aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGATAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTa                                                                              <  2:1440074/100‑1 (MQ=255)
                   ccTGAGAGCATTTTTTTTTTTTTTTATTTCGTTAGTTTAGATTTTGAATATCTTGTTCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTaa                                                                             >  2:1112216/1‑100 (MQ=255)
                    cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAt                                                                            >  1:2259620/1‑100 (MQ=255)
                     tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAtt                                                                           <  2:1609531/100‑1 (MQ=255)
                     tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAtt                                                                           >  2:269846/1‑100 (MQ=255)
                      gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTggagatcggaa                                                                          >  2:897721/1‑91 (MQ=255)
                       agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAtttt                                                                         <  2:223952/100‑1 (MQ=255)
                        gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTta                                                                        <  1:1864968/100‑1 (MQ=255)
                         agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTtat                                                                       >  2:2144264/1‑100 (MQ=255)
                          gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTtata                                                                      >  1:1841959/1‑100 (MQ=255)
                           cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                     <  2:1541232/100‑1 (MQ=255)
                              tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATg                                                                  <  1:2144271/100‑1 (MQ=255)
                              tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATg                                                                  >  1:1545376/1‑100 (MQ=255)
                              tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGCCGGTAATTTTATACATg                                                                  <  2:1371311/100‑1 (MQ=255)
                                  tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAaca                                                              >  2:2127399/1‑100 (MQ=255)
                                   cGCATTCTGATTTCGTTAGTTTAGATTTTGAATGTCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                             >  2:1411023/1‑100 (MQ=255)
                                   cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTTCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTTGCGGTAATTTTATACATGAaca                                                              <  2:2104332/99‑1 (MQ=255)
                                   cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                             >  1:2064949/1‑100 (MQ=255)
                                   cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                             <  1:1685833/100‑1 (MQ=255)
                                    gCATTCTGATTTCGTTAGATTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAcc                                                            >  2:773835/1‑100 (MQ=255)
                                     cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCt                                                           >  2:1337707/1‑100 (MQ=255)
                                     cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCACTTTTCAATAGTGGCGGTAATTTTATACATGAACACCt                                                           >  2:2066658/1‑100 (MQ=255)
                                         cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCat                                                       <  2:153095/100‑1 (MQ=255)
                                           gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGCACACCTCCatta                                                     <  1:1536618/100‑1 (MQ=255)
                                                 gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCCACATTATTAttt                                               >  2:717803/1‑100 (MQ=255)
                                                  ttAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTc                                              >  1:1191942/1‑100 (MQ=255)
                                                   tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCGCCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTAttt                                               >  2:404645/1‑98 (MQ=255)
                                                   tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcc                                             <  2:2259612/100‑1 (MQ=255)
                                                   tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcc                                             <  2:1545372/100‑1 (MQ=255)
                                                   cgatctAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTc                                              <  2:2099079/94‑1 (MQ=255)
                                                     gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTAttt                                               >  2:2269459/1‑96 (MQ=255)
                                                      tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                          <  2:1841954/100‑1 (MQ=255)
                                                       ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTg                                         >  2:871057/1‑100 (MQ=255)
                                                        tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTgg                                        <  1:1021246/100‑1 (MQ=255)
                                                        tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTgg                                        >  2:207001/1‑100 (MQ=255)
                                                        tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTgg                                        <  1:1621837/100‑1 (MQ=255)
                                                         aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcagatcgg                                       >  1:2099086/1‑93 (MQ=255)
                                                          gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                                      >  1:1504724/1‑100 (MQ=255)
                                                             tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                   <  1:569799/100‑1 (MQ=255)
                                                             tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                   >  1:1566784/1‑100 (MQ=255)
                                                             tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                   <  1:1287045/100‑1 (MQ=255)
                                                              ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGtt                                  <  1:892303/100‑1 (MQ=255)
                                                               tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGttt                                 >  2:833520/1‑100 (MQ=255)
                                                                 aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                               <  1:1215018/100‑1 (MQ=255)
                                                                 aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                               <  1:1206751/100‑1 (MQ=255)
                                                                  atatCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                              <  1:2115433/100‑1 (MQ=255)
                                                                   tatCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTc                             <  1:1669116/100‑1 (MQ=255)
                                                                     tCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                           <  2:1995939/100‑1 (MQ=255)
                                                                      cTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                          >  2:1657825/1‑100 (MQ=255)
                                                                       ttGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTCTACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATc                         >  1:682811/1‑100 (MQ=255)
                                                                       ttGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGGGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATc                         <  1:760127/100‑1 (MQ=255)
                                                                        tctgcAGTTAGTAGGAGAGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACTCCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                          <  1:2059774/94‑1 (MQ=255)
                                                                        tGTCCAGTTAGTAGTAGTGCCACCTCCCTTTTCAATATTGGCAGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTTGTTCGTTTATTCCATc                         <  2:1346553/99‑1 (MQ=255)
                                                                        tGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCt                        <  2:288952/100‑1 (MQ=255)
                                                                          tCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                      >  1:1730775/1‑100 (MQ=255)
                                                                           ccAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAtagatc                     >  2:2059767/1‑95 (MQ=255)
                                                                           ccAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                     >  1:281489/1‑100 (MQ=255)
                                                                           ccAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                     >  1:282546/1‑100 (MQ=255)
                                                                           ccAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                     <  2:1730770/100‑1 (MQ=255)
                                                                             aGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                   <  1:773836/100‑1 (MQ=255)
                                                                              gTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                  <  1:459996/100‑1 (MQ=255)
                                                                              gTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCT‑CATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                 <  1:1386422/100‑1 (MQ=255)
                                                                                 agtagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGc               <  1:832410/100‑1 (MQ=255)
                                                                                  gtagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGct              >  2:2048682/1‑100 (MQ=255)
                                                                                   tagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGctt             <  1:833522/100‑1 (MQ=255)
                                                                                    agGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttc            >  2:1675101/1‑100 (MQ=255)
                                                                                         tGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt       <  2:713451/100‑1 (MQ=255)
                                                                                          tccACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc      <  1:263006/99‑1 (MQ=255)
                                                                                          gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc      >  2:258096/1‑100 (MQ=255)
                                                                                          gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc      <  1:2053683/100‑1 (MQ=255)
                                                                                           ccACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGCACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt       >  2:2221053/1‑98 (MQ=255)
                                                                                           ccACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa     >  2:263004/1‑100 (MQ=255)
                                                                                           ccACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa     <  2:1112015/100‑1 (MQ=255)
                                                                                            cACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt       >  2:1802356/1‑97 (MQ=255)
                                                                                            cACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc    <  1:873895/100‑1 (MQ=255)
                                                                                             aCCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACt   >  2:112255/1‑100 (MQ=255)
                                                                                              ccttccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACt   <  2:938480/99‑1 (MQ=255)
                                                                                              cctcccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACtt  >  2:1520009/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTT  >  NC_000913/1413612‑1413805

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