Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,608,595 2 bp→GC intergenic (+23/‑369) prfC → / → osmY peptide chain release factor RF‑3/OsmY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,608,5950CG100.0% 713.3 / NA 225intergenic (+23/‑370)prfC/osmYpeptide chain release factor RF‑3/OsmY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (99/126);  total (99/126)
*NC_0009134,608,5960GC99.6% 721.5 / ‑4.8 225intergenic (+24/‑369)prfC/osmYpeptide chain release factor RF‑3/OsmY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (99/125);  minor base T (0/1);  total (99/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

Only 100 of 229 total aligned reads displayed.

CAACCTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTACGCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCTGATA  >  NC_000913/4608509‑4608682
                                                                                      ||                                                                                      
 aaCCTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGCCCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTgc                                                                                  <  2:3743443/93‑1 (MQ=255)
 aaCCTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTgc                                                                                  >  1:6292047/1‑93 (MQ=255)
 aaCCTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTgc                                                                                  <  2:376242/93‑1 (MQ=255)
    cTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGcc                                                                               <  1:3737489/93‑1 (MQ=255)
    cTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGcc                                                                               >  2:3043316/1‑93 (MQ=255)
     tGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGccc                                                                              <  1:5631429/93‑1 (MQ=255)
      gcgcCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCt                                                                             <  1:5591567/93‑1 (MQ=255)
         ccTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTc                                                                          <  1:4635205/93‑1 (MQ=255)
          cTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCa                                                                         >  2:182910/1‑93 (MQ=255)
          cTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCa                                                                         >  1:1939515/1‑93 (MQ=255)
          cTGGCAAAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCa                                                                         >  2:1979067/1‑93 (MQ=255)
           tGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCaa                                                                        >  1:2104783/1‑93 (MQ=255)
             gCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAAtt                                                                      >  2:3491011/1‑93 (MQ=255)
              cacaGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAAttt                                                                     <  1:2930885/93‑1 (MQ=255)
               acaGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTc                                                                    >  1:5088619/1‑93 (MQ=255)
               acaGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTc                                                                    <  2:5231496/93‑1 (MQ=255)
                 aGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCCATTTcgc                                                                  >  2:3777960/1‑91 (MQ=255)
                  ggAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCcct                                                                 >  1:1493413/1‑93 (MQ=255)
                  ggAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCcct                                                                 >  2:4923402/1‑93 (MQ=255)
                     aCGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCccttcc                                                              >  2:274767/1‑93 (MQ=255)
                      cGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCccttcct                                                             <  1:223657/93‑1 (MQ=255)
                       gTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCccttcctt                                                            >  1:56791/1‑93 (MQ=255)
                       gTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCccttcctt                                                            >  2:5830022/1‑93 (MQ=255)
                        ttATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCtta                                                           >  1:468931/1‑93 (MQ=255)
                        ttATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCtta                                                           <  1:4250418/93‑1 (MQ=255)
                         tATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCttat                                                          <  1:6034010/93‑1 (MQ=255)
                          aTCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCttatt                                                         >  1:4791200/1‑93 (MQ=255)
                          aTCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCttatt                                                         >  1:608571/1‑93 (MQ=255)
                          aTCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCttatt                                                         <  2:4318260/93‑1 (MQ=255)
                          aTCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCttatt                                                         >  1:1353694/1‑93 (MQ=255)
                           tCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCACTTCCttatta                                                        >  2:5615936/1‑93 (MQ=255)
                            ccGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAg                                                       <  2:5082978/93‑1 (MQ=255)
                              ggACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGcc                                                     <  1:4550987/93‑1 (MQ=255)
                              ggACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGcc                                                     >  2:5689558/1‑93 (MQ=255)
                               gACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCg                                                    >  2:4579558/1‑93 (MQ=255)
                                 cGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCt                                                  <  2:6123033/93‑1 (MQ=255)
                                  gTTCAGTTCCACCAGCCCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCtt                                                 <  2:4052002/93‑1 (MQ=255)
                                  gTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCtt                                                 <  1:5814673/93‑1 (MQ=255)
                                  gTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCtt                                                 >  2:725987/1‑93 (MQ=255)
                                   ttCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTa                                                <  2:5376246/93‑1 (MQ=255)
                                    tCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTAc                                               >  2:3137330/1‑93 (MQ=255)
                                    tCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTAc                                               >  1:371469/1‑93 (MQ=255)
                                    tCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTAc                                               >  1:5970760/1‑93 (MQ=255)
                                    tCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTAc                                               <  1:2123497/93‑1 (MQ=255)
                                    tCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTAc                                               <  2:3624624/93‑1 (MQ=255)
                                     taGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACg                                              <  1:4149434/92‑1 (MQ=255)
                                      aGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACgg                                             <  1:496705/93‑1 (MQ=255)
                                       gTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGa                                            >  2:100371/1‑93 (MQ=255)
                                          ccaccaGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATg                                         <  2:4110335/93‑1 (MQ=255)
                                          ccaccaGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATg                                         >  2:2202208/1‑93 (MQ=255)
                                           caccaGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGt                                        <  2:4707291/93‑1 (MQ=255)
                                            cccaGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGtt                                       <  2:1444616/92‑1 (MQ=255)
                                             ccaGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTc                                      <  1:3148724/93‑1 (MQ=255)
                                             ccaGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTc                                      <  2:809548/93‑1 (MQ=255)
                                              caGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCt                                     <  1:5317487/93‑1 (MQ=255)
                                              caGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCt                                     <  1:2641305/93‑1 (MQ=255)
                                               aGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCtt                                    <  2:853243/93‑1 (MQ=255)
                                                gACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTa                                   >  2:4172010/1‑93 (MQ=255)
                                                   ccGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTaaaa                                <  2:5643814/93‑1 (MQ=255)
                                                   ccGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTaaaa                                >  1:2787268/1‑93 (MQ=255)
                                                    cgcgAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAAc                               <  2:634749/93‑1 (MQ=255)
                                                     gcgAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACa                              >  2:4273747/1‑93 (MQ=255)
                                                     gcgAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACa                              <  2:2417514/93‑1 (MQ=255)
                                                      cgAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACAt                             <  1:6235146/93‑1 (MQ=255)
                                                      cgAGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACAt                             <  2:3537155/93‑1 (MQ=255)
                                                        aGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTc                           >  1:2888407/1‑93 (MQ=255)
                                                        aGCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTc                           <  1:5554584/93‑1 (MQ=255)
                                                         gCATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCa                          <  1:1580128/93‑1 (MQ=255)
                                                         acATTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTACCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCa                          <  2:4949099/92‑1 (MQ=255)
                                                          cATTAATTCTTGCCGCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCAc                         <  1:5057425/93‑1 (MQ=255)
                                                           aTTAATTCTTGCCTCCAGGGCGCGTTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACt                        <  1:2592518/93‑1 (MQ=255)
                                                           aTTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACt                        >  1:3599036/1‑93 (MQ=255)
                                                           aTTAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACt                        >  1:4354992/1‑93 (MQ=255)
                                                            ttAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACtt                       >  2:4982118/1‑93 (MQ=255)
                                                            ttAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACtt                       >  1:3775356/1‑93 (MQ=255)
                                                             tAATTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACttt                      <  1:3966743/93‑1 (MQ=255)
                                                              aaTTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACtttt                     >  1:1692845/1‑93 (MQ=255)
                                                               aTTCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTg                    <  1:6052654/93‑1 (MQ=255)
                                                                 tCTTGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCt                  <  1:104550/93‑1 (MQ=255)
                                                                   ttGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTa                <  1:790065/93‑1 (MQ=255)
                                                                   ttGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTa                <  2:1979788/93‑1 (MQ=255)
                                                                    tGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTAt               <  1:52614/93‑1 (MQ=255)
                                                                    tGCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTAt               >  1:3307036/1‑93 (MQ=255)
                                                                     gCCTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATg              <  1:2159508/93‑1 (MQ=255)
                                                                      ccTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGt             <  2:1264707/93‑1 (MQ=255)
                                                                      ccTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGt             <  2:1057992/93‑1 (MQ=255)
                                                                       cTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAGAACATTCACTTTTGCTTATGtt            >  1:6112383/1‑93 (MQ=255)
                                                                       cTCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGtt            <  2:3232701/93‑1 (MQ=255)
                                                                        tCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGttt           >  2:1595514/1‑93 (MQ=255)
                                                                        tCCAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGttt           >  2:5456040/1‑93 (MQ=255)
                                                                          cAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCCCTTTTGCTTATGTTTTc         >  2:1221850/1‑93 (MQ=255)
                                                                          cAGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTc         <  2:2533751/93‑1 (MQ=255)
                                                                           aGGGCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCg        >  2:4073716/1‑93 (MQ=255)
                                                                            gggCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGc       >  1:122833/1‑93 (MQ=255)
                                                                             ggCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCt      >  2:192120/1‑93 (MQ=255)
                                                                             ggCGCGGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCt      <  1:1146570/93‑1 (MQ=255)
                                                                                gcgGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCTGat   >  1:5167969/1‑93 (MQ=255)
                                                                                 cgGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCTGata  <  2:3331548/93‑1 (MQ=255)
                                                                                 cgGTAGCCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCTGata  <  1:5467325/93‑1 (MQ=255)
                                                                                 cgGTAGCCGCTGCGCCCTGTCAATTGCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCTGata  <  1:3012865/93‑1 (MQ=255)
                                                                                      ||                                                                                      
CAACCTGCGCCTGGCACAGGAACGTTATCCGGACGTTCAGTTCCACCAGACCCGCGAGCATTAATTCTTGCCTCCAGGGCGCGGTACGCGCTGCGCCCTGTCAATTTCCCTTCCTTATTAGCCGCTTACGGAATGTTCTTAAAACATTCACTTTTGCTTATGTTTTCGCTGATA  >  NC_000913/4608509‑4608682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: