Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,655,680 Δ1 bp coding (679/1284 nt) pepB ← aminopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,655,6770C.89.5% 1079.5 / 129.0 266coding (682/1284 nt)pepBaminopeptidase B
Reads supporting (aligned to +/- strand):  ref base C (22/6);  new base . (118/120);  total (140/126)
Fisher's exact test for biased strand distribution p-value = 4.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 295 total aligned reads displayed.

CGCCCAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTT  >  NC_000913/2655578‑2655776
                                                                                                   |                                                                                                   
cGCCCAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcc                                                                                                    <  2:4785687/101‑1 (MQ=255)
cGCCCAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcc                                                                                                    <  1:4785863/101‑1 (MQ=255)
 gCCCAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGccc                                                                                                   >  1:4786342/1‑101 (MQ=255)
   ccAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcccgg                                                                                                 <  2:4785799/101‑3 (MQ=255)
   ccAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCg                                                                                                 >  1:4786346/1‑101 (MQ=255)
   ccAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCg                                                                                                 >  2:4786347/1‑101 (MQ=255)
    cAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcccggt                                                                                                >  2:4786352/1‑98 (MQ=255)
     aGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcccggta                                                                                               >  1:4786354/1‑97 (MQ=255)
     aGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcccggta                                                                                               <  1:4785836/101‑5 (MQ=255)
       cTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAAc                                                                                            >  1:4786359/1‑101 (MQ=255)
       cTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAAc                                                                                            <  1:4785654/101‑1 (MQ=255)
        ttGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAAcc                                                                                           >  1:4786363/1‑101 (MQ=255)
         tGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCg                                                                                          >  2:4786364/1‑101 (MQ=255)
         tGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCg                                                                                          <  2:4785834/101‑1 (MQ=255)
          gAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGt                                                                                         >  2:4786369/1‑101 (MQ=255)
          gAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGt                                                                                         >  1:4786371/1‑101 (MQ=255)
           aaCGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGtt                                                                                        <  1:4785887/101‑1 (MQ=255)
           aaCGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGtt                                                                                        >  2:4786376/1‑101 (MQ=255)
            aCGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTg                                                                                       >  1:4786379/1‑101 (MQ=255)
            aCGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTg                                                                                       <  2:4785924/101‑1 (MQ=255)
              gCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTgc                                                                                      >  2:4786382/1‑101 (MQ=255)
               cATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccg                                                                                    >  1:4786384/1‑101 (MQ=255)
                  tGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccgcgc                                                                                 >  1:4786388/1‑101 (MQ=255)
                   gCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTgccgcgc                                                                                 >  2:4786391/1‑101 (MQ=255)
                   gCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccgtgcc                                                                                <  1:4785856/101‑1 (MQ=255)
                   gCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccgcgcc                                                                                <  1:4786182/101‑1 (MQ=255)
                    ccGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccgcgccg                                                                               <  2:4786190/101‑1 (MQ=255)
                     cGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccgcgccgc                                                                              >  1:4786395/1‑101 (MQ=255)
                     cGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTgccgcgccgc                                                                              >  2:4786397/1‑101 (MQ=255)
                       cTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGccc                                                                            <  2:4786210/101‑1 (MQ=255)
                         aaTCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCAt                                                                          <  1:4785613/101‑1 (MQ=255)
                         aaTCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCAt                                                                          >  2:4786401/1‑101 (MQ=255)
                          aTCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCAt                                                                          >  2:4786406/1‑101 (MQ=255)
                          aTCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATg                                                                         >  1:4786402/1‑101 (MQ=255)
                             aGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTcc                                                                      >  2:4786411/1‑101 (MQ=255)
                             aGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTcc                                                                      <  2:4785883/101‑1 (MQ=255)
                             aGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTcc                                                                      >  1:4786410/1‑101 (MQ=255)
                             aGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTcc                                                                      <  1:4785930/101‑1 (MQ=255)
                             aGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTcc                                                                      <  2:4785659/101‑1 (MQ=255)
                              ggTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGcccggt                                                                                                >  1:4786414/1‑72 (MQ=255)
                              ggTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTcc                                                                      <  1:4786004/101‑1 (MQ=255)
                                ttATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGAc                                                                   <  1:4786142/101‑1 (MQ=255)
                                 tATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACt                                                                  >  2:4786419/1‑101 (MQ=255)
                                  aTCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACtt                                                                 >  1:4786420/1‑101 (MQ=255)
                                   tCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTc                                                                <  2:4785871/101‑1 (MQ=255)
                                   tCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTc                                                                >  1:4786422/1‑101 (MQ=255)
                                    ccGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCa                                                               <  2:4785564/101‑1 (MQ=255)
                                     cGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCa                                                               >  1:4786427/1‑101 (MQ=255)
                                      gCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATc                                                             >  2:4786429/1‑101 (MQ=255)
                                        acaGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGa                                                           >  1:4786433/1‑101 (MQ=255)
                                          agcagAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGt                                                         <  1:4785972/101‑1 (MQ=255)
                                          agcagAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGt                                                         <  1:4785721/101‑1 (MQ=255)
                                           gcagAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTc                                                        <  2:4785806/101‑1 (MQ=255)
                                            cagAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTcc                                                       <  2:4785765/101‑1 (MQ=255)
                                            cagAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTcc                                                       >  2:4786438/1‑101 (MQ=255)
                                             agagGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCa                                                      >  1:4786440/1‑101 (MQ=255)
                                               agGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATa                                                    <  1:4785626/101‑1 (MQ=255)
                                                 gAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCAt                                                     <  1:4785953/98‑1 (MQ=255)
                                                 gAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATaaa                                                  <  2:4785860/101‑1 (MQ=255)
                                                   aCAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACg                                                >  2:4786443/1‑101 (MQ=255)
                                                     aGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATaaa                                                  >  2:4786446/1‑97 (MQ=255)
                                                       cTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGt                                                                                         >  2:4786450/1‑56 (MQ=255)
                                                       cTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGt                                            <  1:4786108/101‑1 (MQ=255)
                                                        ttCACGCGCTTGTTCAGTCCCCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTc                                           >  1:4786452/1‑101 (MQ=255)
                                                         tCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCt                                          <  1:4785678/101‑1 (MQ=255)
                                                          cACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTg                                         <  1:4785910/101‑1 (MQ=255)
                                                          cACGCGCTTGTTCAGTCCACGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCt                                          >  2:4786454/1‑101 (MQ=255)
                                                            cgcgCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGtt                                       <  2:4786161/101‑1 (MQ=255)
                                                             gcgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGttt                                      >  2:4786458/1‑101 (MQ=255)
                                                              cgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTg                                     >  1:4786462/1‑101 (MQ=255)
                                                                 ttGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATg                                  >  2:4786467/1‑101 (MQ=255)
                                                                 ttGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATg                                  >  1:4786465/1‑101 (MQ=255)
                                                                 ttGTTCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATg                                  >  1:4786469/1‑101 (MQ=255)
                                                                    ttCAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTg                               <  1:4785875/101‑1 (MQ=255)
                                                                      cAGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTa                             <  1:4785899/101‑1 (MQ=255)
                                                                       aGTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAg                            >  2:4786474/1‑101 (MQ=255)
                                                                        gTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAg                            >  2:4786478/1‑101 (MQ=255)
                                                                        gTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAgc                           <  2:4785932/101‑1 (MQ=255)
                                                                        gTCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAgc                           >  1:4786475/1‑101 (MQ=255)
                                                                         tCCGCGCGTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTg                               <  1:4786043/96‑1 (MQ=255)
                                                                          ccGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAgcc                          <  2:4785987/101‑1 (MQ=255)
                                                                            gcgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTc                                                        >  1:4786483/1‑69 (MQ=255)
                                                                               cgTAATGGCAAATGCCCGCG‑CCCGGTAACCGTTgccgcgc                                                                                 >  2:4786485/1‑40 (MQ=37)
                                                                                gTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAg                   <  1:4786076/101‑1 (MQ=255)
                                                                                gTAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAg                   >  1:4786487/1‑101 (MQ=255)
                                                                                 tAATGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAgc                  <  2:4785852/101‑2 (MQ=255)
                                                                                  aaTGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGt                  >  2:4786491/1‑101 (MQ=255)
                                                                                   aTGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCa                <  2:4786021/101‑1 (MQ=255)
                                                                                    tGGCAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCaa               <  2:4786103/101‑1 (MQ=255)
                                                                                       cAAATGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAg            <  1:4785792/101‑1 (MQ=255)
                                                                                        aaaTGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAg            >  1:4786494/1‑101 (MQ=255)
                                                                                        aaaTGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGt           <  2:4786064/101‑1 (MQ=255)
                                                                                        aaaTGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGt           <  1:4786027/101‑1 (MQ=255)
                                                                                        aaaTGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGt           <  1:4785967/101‑1 (MQ=255)
                                                                                        aaaTGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGt           <  1:4785848/101‑1 (MQ=255)
                                                                                        aaaTGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGt           <  1:4785762/101‑1 (MQ=255)
                                                                                         aaTGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCaa               <  2:4786145/96‑1 (MQ=255)
                                                                                          aTGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTg          >  1:4786497/1‑101 (MQ=255)
                                                                                           tGCCAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGAt        <  1:4786036/101‑1 (MQ=255)
                                                                                              cAGCG‑CCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATAcc     >  1:4786501/1‑101 (MQ=255)
                                                                                                   |                                                                                                   
CGCCCAGCTTGAACGCATTGCCGCTAATCAGGTTATCCGCACAGCAGAGGAACAGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTT  >  NC_000913/2655578‑2655776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑