Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,709 T→C 23.5% 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,7090TC23.5% 274.8 / 86.0 156E11G (GAA→GGA) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base T (64/55);  new base C (19/18);  total (83/73)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

GAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCA  >  NC_000913/1413610‑1413808
                                                                                                   |                                                                                                   
gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctt                                                                                                     >  2:3076533/1‑100 (MQ=255)
    cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAccttcctt                                                                                                 >  1:3237928/1‑100 (MQ=255)
     gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCttt                                                                                                >  1:2955814/1‑100 (MQ=255)
     gCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCttt                                                                                                >  1:169841/1‑100 (MQ=255)
      ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCtttt                                                                                               <  2:1549399/100‑1 (MQ=255)
       ccgatctttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTc                                                                                              <  2:3162471/93‑1 (MQ=255)
       cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTc                                                                                              <  1:1305835/100‑1 (MQ=255)
       cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACATTCCTTTTc                                                                                              >  1:1895235/1‑100 (MQ=255)
        tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCa                                                                                             >  2:841456/1‑100 (MQ=255)
        tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCa                                                                                             <  1:3669935/100‑1 (MQ=255)
        tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCa                                                                                             <  1:731383/100‑1 (MQ=255)
         gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCaa                                                                                            <  2:4168345/100‑1 (MQ=255)
         gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCaa                                                                                            <  2:1915163/100‑1 (MQ=255)
         gCTTCTTTGTTACCTGAGAGCATTTTTTCGCATTCGGATTTCGTTAGTTTAGATTTTGACTATCTTGCCCAGTAAGTAGGAGTGCCAccttcc                                                                                                   <  2:106150/93‑1 (MQ=255)
          ttttttttttACTTGAGAGCAGTTTTTCGCATTCGGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCTACCTCCCTTTTCAAt                                                                                           >  2:1650848/5‑100 (MQ=255)
          cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAAt                                                                                           >  1:2243424/1‑100 (MQ=255)
          cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAGTTTGCATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCaa                                                                                            >  1:2911309/1‑99 (MQ=255)
          cttcttTTTTACCTGAGAGCATTTTTTCACATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAAt                                                                                           >  1:3033829/1‑100 (MQ=255)
              ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGt                                                                                        >  2:4090517/1‑99 (MQ=255)
                ttttACCTGATAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGc                                                                                     <  2:3306974/100‑1 (MQ=255)
                ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGc                                                                                     <  1:2248129/100‑1 (MQ=255)
                 tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCg                                                                                    <  2:1372173/100‑1 (MQ=255)
                   tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCgat                                                                                  >  1:2148835/1‑98 (MQ=255)
                     ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTaa                                                                                >  2:2725431/1‑100 (MQ=255)
                     ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTaa                                                                                >  1:154312/1‑100 (MQ=255)
                       tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAtt                                                                              >  1:1895118/1‑100 (MQ=255)
                        gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAttt                                                                             <  2:2955797/100‑1 (MQ=255)
                          gagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTta                                                                           >  1:3635374/1‑100 (MQ=255)
                           agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTtat                                                                          <  2:300233/100‑1 (MQ=255)
                            gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCAACCTTCCTTTTCAATAGTGGCGGTAATTTtata                                                                         >  2:3599835/1‑100 (MQ=255)
                             cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                        <  1:779143/100‑1 (MQ=255)
                              aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACa                                                                       >  2:2976213/1‑100 (MQ=255)
                              aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACa                                                                       <  2:2634696/100‑1 (MQ=255)
                              aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACa                                                                       >  1:3842493/1‑100 (MQ=255)
                                tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATg                                                                     <  2:2784328/100‑1 (MQ=255)
                                 ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAt                                                                      >  2:986650/1‑98 (MQ=255)
                                 ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                    <  2:2013369/100‑1 (MQ=255)
                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                  <  2:1492998/100‑1 (MQ=255)
                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                  <  2:3843633/100‑1 (MQ=255)
                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                  <  1:2855870/100‑1 (MQ=255)
                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                                  <  2:2316924/100‑1 (MQ=255)
                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATGTTATACATGAac                                                                  >  2:3444956/1‑100 (MQ=255)
                                    tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAaca                                                                 <  1:3132051/100‑1 (MQ=255)
                                     cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                                >  2:2055826/1‑100 (MQ=255)
                                      gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAcc                                                               <  2:3831150/100‑1 (MQ=255)
                                       cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCt                                                              >  2:213593/1‑100 (MQ=255)
                                        aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTc                                                             <  2:133355/100‑1 (MQ=255)
                                         ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTcc                                                            >  1:595518/1‑100 (MQ=255)
                                         ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTcc                                                            <  1:442898/100‑1 (MQ=255)
                                             gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCatta                                                        <  1:1937082/100‑1 (MQ=255)
                                             gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACTTCCatta                                                        >  2:2537871/1‑100 (MQ=255)
                                                  cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATTAACACCTCAattattatt                                                   >  2:2174799/1‑100 (MQ=255)
                                                  cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattatt                                                   >  1:3931848/1‑100 (MQ=255)
                                                     tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTcc                                                <  2:1743042/100‑1 (MQ=255)
                                                       gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAg                                              >  1:2665824/1‑100 (MQ=255)
                                                       gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAg                                              >  1:2409564/1‑100 (MQ=255)
                                                       gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCAACTTTACTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAg                                              >  2:3535154/1‑100 (MQ=255)
                                                        tttAGCTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCATTTTCAATAGTGGCGGTAATTTTATCCATGACCACCTCCCTTATTATTTCCAg                                              >  2:3941923/1‑99 (MQ=255)
                                                        tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAACTCCATTATTATTTCCAGt                                             >  2:1201747/1‑100 (MQ=255)
                                                        tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                             <  1:2259875/100‑1 (MQ=255)
                                                         ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTg                                            <  2:1250221/100‑1 (MQ=255)
                                                           aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGt                                          <  2:539282/100‑1 (MQ=255)
                                                             aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTc                                        <  1:3221983/100‑1 (MQ=255)
                                                               tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                      >  1:2595950/1‑100 (MQ=255)
                                                                ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCCTTATTATTTCCAGTGGTTCGtt                                     <  2:2420623/100‑1 (MQ=255)
                                                                ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGtt                                     <  2:2122254/100‑1 (MQ=255)
                                                                   aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                                  >  1:606637/1‑100 (MQ=255)
                                                                    atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                                 <  1:174182/100‑1 (MQ=255)
                                                                    atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                                 >  2:1396147/1‑100 (MQ=255)
                                                                    atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                                 >  2:128437/1‑100 (MQ=255)
                                                                      atCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTcc                               <  1:2395011/100‑1 (MQ=255)
                                                                        cTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt                             >  1:99844/1‑100 (MQ=255)
                                                                           gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCtt                          >  1:3081370/1‑100 (MQ=255)
                                                                            tCCAGTTCGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCGTTATTATTTCCAGTGGTTCGTTTATTCCATCttt                         >  2:2182708/1‑100 (MQ=255)
                                                                            tCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                         >  1:2518755/1‑100 (MQ=255)
                                                                             ccAGTTAGTAGGAGTGCCACCTTCCTTGTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                         >  2:919254/1‑99 (MQ=255)
                                                                              cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg                       >  1:1032180/1‑100 (MQ=255)
                                                                               aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                      <  1:121133/100‑1 (MQ=255)
                                                                                gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                     >  1:2102611/1‑100 (MQ=255)
                                                                                gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                     <  2:606633/100‑1 (MQ=255)
                                                                                 ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                    <  1:2604807/100‑1 (MQ=255)
                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGct                 <  2:1037910/100‑1 (MQ=255)
                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGct                 <  2:215453/100‑1 (MQ=255)
                                                                                      agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttc               <  1:1105834/100‑1 (MQ=255)
                                                                                          gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt           <  2:865580/100‑1 (MQ=255)
                                                                                           tGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttt            <  2:2452100/98‑1 (MQ=255)
                                                                                            gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc         <  1:4009336/100‑1 (MQ=255)
                                                                                            gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc         >  1:760189/1‑100 (MQ=255)
                                                                                            gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc         <  2:3033811/100‑1 (MQ=255)
                                                                                             ccACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc         >  2:20422/1‑99 (MQ=255)
                                                                                             ccACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCa        >  2:484921/1‑100 (MQ=255)
                                                                                              cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc       >  1:1986022/1‑100 (MQ=255)
                                                                                              cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc       <  2:1170772/100‑1 (MQ=255)
                                                                                              cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACACGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc       <  2:4064874/100‑1 (MQ=255)
                                                                                              cACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc       <  1:1595309/100‑1 (MQ=255)
                                                                                               aCCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACt      >  1:307145/1‑100 (MQ=255)
                                                                                                  ttccttTTCAATAGTGTCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcc   <  2:1810112/100‑1 (MQ=255)
                                                                                                  ttccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcc   >  1:1174015/1‑100 (MQ=255)
                                                                                                   tccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcca  <  1:3291023/100‑1 (MQ=255)
                                                                                                   tccttTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTcca  >  1:882842/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCA  >  NC_000913/1413610‑1413808

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: