Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,982 G→C 17.9% intergenic (+28/+32) yhcG → / ← yhcH DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,9820GC17.9% 440.6 / 91.6 179intergenic (+28/+32)yhcG/yhcHDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily
Reads supporting (aligned to +/- strand):  ref base G (83/64);  new base C (22/10);  total (105/74)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 187 total aligned reads displayed.

GATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACG  >  NC_000913/3368883‑3369080
                                                                                                   |                                                                                                  
gatgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTTCAGACAg                                                                                                    >  2:283482/1‑100 (MQ=255)
  tgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtg                                                                                                  >  1:2714249/1‑100 (MQ=255)
  tgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAAAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtg                                                                                                  >  1:2909785/1‑100 (MQ=255)
   gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                 >  2:2474831/1‑100 (MQ=255)
     gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgaga                                                                                               >  2:3347489/1‑97 (MQ=255)
     gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCt                                                                                               >  1:1684942/1‑100 (MQ=255)
      aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTg                                                                                              <  1:4054019/100‑1 (MQ=255)
       aaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGc                                                                                             >  1:4004225/1‑100 (MQ=255)
       aaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGc                                                                                             >  1:1900751/1‑100 (MQ=255)
        aCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCa                                                                                            <  1:1574840/100‑1 (MQ=255)
            cTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                        >  1:1473738/1‑100 (MQ=255)
            cTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                        <  1:95869/100‑1 (MQ=255)
            cTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                         <  1:3132307/99‑1 (MQ=255)
             tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATc                                                                                       >  2:3285266/1‑100 (MQ=255)
              cGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCt                                                                                      <  1:1535978/100‑1 (MQ=255)
              cGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCt                                                                                      <  1:3335789/100‑1 (MQ=255)
               gCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCGGACACTGTCTGCAGAAttcg                                                                                     <  2:2561135/100‑4 (MQ=255)
                cAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCttt                                                                                    >  1:2106205/1‑100 (MQ=255)
                 aGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTa                                                                                   <  1:135796/100‑1 (MQ=255)
                  gATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaa                                                                                  >  1:2785436/1‑100 (MQ=255)
                  gATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGAAttcgaaa                                                                                  <  2:2459616/100‑7 (MQ=255)
                   atctaCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaaa                                                                                 <  1:1996942/96‑1 (MQ=255)
                      aaCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTGAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACg                                                                              <  1:1625141/100‑1 (MQ=255)
                       aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaaaagat                                                                             >  2:1996939/1‑97 (MQ=255)
                       aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGa                                                                             >  2:3749032/1‑100 (MQ=255)
                        cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                            >  1:108193/1‑100 (MQ=255)
                         tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAAc                                                                           <  1:3486560/100‑1 (MQ=255)
                         tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAAc                                                                           >  2:2626747/1‑100 (MQ=255)
                             tCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAg                                                                       <  1:2676660/100‑1 (MQ=255)
                               gAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGCACTAAGt                                                                      >  1:2269737/1‑99 (MQ=255)
                                aaCACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGttt                                                                    >  2:125120/1‑100 (MQ=255)
                                 acagAGGCGGTGCTGGAGGAAGGCTATGGAGGGCGTTAATTTCAGAGGTTTGGAATTCTGAAGACAGTGTCTGCAGAATCTTGAAAACGAACTAAGttt                                                                    <  2:1165505/99‑1 (MQ=255)
                                   acaGGCGGTGCTAGAGGAAAGCTATAGACGCCGCTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                 >  2:2740669/1‑100 (MQ=255)
                                    caGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                >  2:4047191/1‑100 (MQ=255)
                                    caGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                <  1:1527284/100‑1 (MQ=255)
                                    caGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                >  2:302023/1‑100 (MQ=255)
                                      ggCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGcc                                                              >  2:352587/1‑100 (MQ=255)
                                        cGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAt                                                            <  1:4107366/100‑1 (MQ=255)
                                         ggTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAtt                                                           <  1:632476/100‑1 (MQ=255)
                                           tGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTaa                                                         >  2:2646659/1‑100 (MQ=255)
                                           tGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTaa                                                         >  2:1111838/1‑100 (MQ=255)
                                           tGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTaa                                                         <  1:3945627/100‑1 (MQ=255)
                                            gCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAACCTTTAAAACGAACTAAGTTTTAAGCCATTaaa                                                        >  1:3810172/1‑100 (MQ=255)
                                             cTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAt                                                       >  1:3694579/1‑100 (MQ=255)
                                              tAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACCAAGTTTTAAGCCATTAAATc                                                      <  2:2566515/100‑1 (MQ=255)
                                               agatGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCATAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCa                                                     <  2:3885041/100‑1 (MQ=255)
                                               agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCa                                                     <  1:1692638/100‑1 (MQ=255)
                                               agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCa                                                     >  1:699047/1‑100 (MQ=255)
                                                gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGACTCTTTAAAACGAACTAAGTTTTAAGCCCTTAAATc                                                      >  2:3878516/1‑98 (MQ=255)
                                                gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAg                                                    <  2:2895909/100‑1 (MQ=255)
                                                  ggAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGcc                                                  >  1:2404926/1‑100 (MQ=255)
                                                  ggAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGcc                                                  <  2:1900748/100‑1 (MQ=255)
                                                  ggAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGcc                                                  <  1:2071828/100‑1 (MQ=255)
                                                  ggAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGcc                                                  >  1:538532/1‑100 (MQ=255)
                                                        caataGACGCCGTCAATTTCAGCGGTTTCGAGTTCTGCAGACGGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                            >  1:2831721/3‑100 (MQ=255)
                                                        cTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                            >  1:939444/1‑100 (MQ=255)
                                                         tataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGTAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCt                                           >  1:501655/1‑100 (MQ=255)
                                                           taGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGAAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAACTTc                                         >  1:1982591/1‑100 (MQ=255)
                                                            aGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTca                                        <  1:2361254/100‑1 (MQ=255)
                                                             gACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcac                                       >  1:4197968/1‑100 (MQ=255)
                                                               cGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                      <  2:1670715/99‑1 (MQ=255)
                                                                   gTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAAc                                 <  1:3269420/100‑1 (MQ=255)
                                                                   gTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAAc                                 >  1:2996249/1‑100 (MQ=255)
                                                                     tAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCt                               >  2:3661007/1‑100 (MQ=255)
                                                                      aaTTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtt                              >  1:1234517/1‑100 (MQ=255)
                                                                        tttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                            >  1:785315/1‑100 (MQ=255)
                                                                        tttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                            <  1:2022002/100‑1 (MQ=255)
                                                                        tttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                            >  2:1429592/1‑100 (MQ=255)
                                                                        tttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                            <  1:3491009/100‑1 (MQ=255)
                                                                         ttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                           <  1:3655186/100‑1 (MQ=255)
                                                                          tCAGCGGTTTCGGATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                          <  1:850096/100‑1 (MQ=255)
                                                                          tCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCGCCACAACCtttttn                          <  2:1826188/100‑2 (MQ=255)
                                                                             gCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtc                       <  2:329252/100‑1 (MQ=255)
                                                                               ggTATCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctc                     <  2:3467346/100‑1 (MQ=255)
                                                                                 tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctctc                   <  1:1809904/100‑1 (MQ=255)
                                                                                 tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctctc                   <  2:538531/100‑1 (MQ=255)
                                                                                 tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCAATTAAATCAGCATTAACGTTCACCACAACCTTTTTAAtctctc                   >  1:2012523/1‑100 (MQ=255)
                                                                                 tttCGAATTCTGCAGACAGCGTTTGCAGAATCTTTAAAACGAACTAAGTTTTAATCCATAAAATCATCCTTAACCATCACCACAACCTTTTTAAtc                       >  2:2918328/1‑96 (MQ=255)
                                                                                  ttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTcc                  >  2:80453/1‑100 (MQ=255)
                                                                                    cGAATTCTGCAGACTGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAGTCTCTCCa                 <  2:3162774/99‑1 (MQ=255)
                                                                                    cGAATTCTGCAGACAGTGTCTGCAGAATCTTTACAACGAACTAAGTTTTACGACATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAg                <  2:129173/100‑1 (MQ=255)
                                                                                    aagaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACAACAACCTTTTTAATCTCTCCAg                >  2:3909312/4‑100 (MQ=255)
                                                                                     agaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCa                 >  2:2601803/3‑98 (MQ=255)
                                                                                      gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACAACAACCTTTTTAATCTCTTCAGGc              >  1:323337/2‑100 (MQ=255)
                                                                                      aaTTCTGCAGACAGTGTCTGCAGAATATTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGACTTAACATTCACAAAAACATTTTTAAtatatcaa                 >  2:2221447/1‑90 (MQ=255)
                                                                                       aTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTACCCACAACCTTTTTAATCTCTCCAGGc              >  1:3787961/1‑99 (MQ=255)
                                                                                        ttCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTc            >  1:1420141/1‑100 (MQ=255)
                                                                                        ttCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTc            >  2:1139436/1‑100 (MQ=255)
                                                                                          cTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGGTTTAAGCTATGAAATCTGCTTTAACCTTCACCAAACCCTTTTTAATCGCTCCAGGCTccc          >  2:2356028/1‑98 (MQ=255)
                                                                                          cTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc          >  2:4187666/1‑100 (MQ=255)
                                                                                          cTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc          >  1:2529433/1‑100 (MQ=255)
                                                                                          cTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACAACAACCTTTTTAATCTCTCCAGGCTCg           >  2:2129849/1‑99 (MQ=255)
                                                                                           tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGcc         >  1:3042040/1‑100 (MQ=255)
                                                                                           tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc          >  2:655019/1‑99 (MQ=255)
                                                                                           tGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGcc         >  1:2154472/1‑100 (MQ=255)
                                                                                            gCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCg        >  1:1134411/1‑100 (MQ=255)
                                                                                            gCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCg        >  1:3067678/1‑100 (MQ=255)
                                                                                              aGACAGTGTCTCCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAACTCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCg        <  2:3806171/98‑1 (MQ=255)
                                                                                               gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACa     >  1:3879188/1‑100 (MQ=255)
                                                                                                 cAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc   <  1:3039778/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACG  >  NC_000913/3368883‑3369080

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: