Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 50,798 A→G V142A (GTG→GCG)  apaH ← diadenosine tetraphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091350,7980AG99.6% 823.3 / ‑4.7 235V142A (GTG→GCG) apaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (128/105);  total (129/106)
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.17e-01

Only 100 of 239 total aligned reads displayed.

TTCCCAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCACCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCT  >  NC_000913/50709‑50887
                                                                                         |                                                                                         
ttCCCAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCg                                                                                        >  1:2567333/1‑93 (MQ=255)
ttCCCAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCg                                                                                        >  2:2783854/1‑93 (MQ=255)
   ccAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCtt                                                                                     >  1:1411016/1‑93 (MQ=255)
    cAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTc                                                                                    <  1:1171276/93‑1 (MQ=255)
    cAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTc                                                                                    <  1:4191686/93‑1 (MQ=255)
    cAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTc                                                                                    >  2:3772747/1‑93 (MQ=255)
     aGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCt                                                                                   >  2:4762140/1‑93 (MQ=255)
      gCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTa                                                                                  <  2:3749950/93‑1 (MQ=255)
       cccccGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTAc                                                                                 <  1:5562056/93‑1 (MQ=255)
          ccGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATc                                                                              <  1:226191/93‑1 (MQ=255)
           cGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATcg                                                                             >  2:2383500/1‑93 (MQ=255)
            gCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATcgc                                                                            <  2:2846660/93‑1 (MQ=255)
             cAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATcgcg                                                                           >  2:2558987/1‑93 (MQ=255)
              aaTTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGt                                                                          <  2:3382576/93‑1 (MQ=255)
              aaTTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGt                                                                          >  2:5644214/1‑93 (MQ=255)
              aaTTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGt                                                                          >  2:854053/1‑93 (MQ=255)
               aTTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTg                                                                         >  2:4113768/1‑93 (MQ=255)
               aTTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTg                                                                         >  1:4009399/1‑93 (MQ=255)
                ttCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTgc                                                                        >  2:1225231/1‑93 (MQ=255)
                   cGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCa                                                                     <  2:1472684/93‑1 (MQ=255)
                   cGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCa                                                                     <  1:1885247/93‑1 (MQ=255)
                    ggTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCAc                                                                    >  1:4818023/1‑93 (MQ=255)
                     gTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCAct                                                                   <  2:5345085/93‑1 (MQ=255)
                       gACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCActct                                                                 <  1:3737045/93‑1 (MQ=255)
                         ccAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCttt                                                               <  2:813980/93‑1 (MQ=255)
                          cAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTg                                                              >  2:1996063/1‑93 (MQ=255)
                            gTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGc                                                            <  1:2118159/93‑1 (MQ=255)
                             ttATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGCTAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCg                                                           >  2:4995399/1‑93 (MQ=255)
                              tATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCgg                                                          >  1:3662041/1‑93 (MQ=255)
                              tATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCgg                                                          <  2:4205666/93‑1 (MQ=255)
                              tATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCgg                                                          >  1:2909729/1‑93 (MQ=255)
                                tttGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTc                                                        >  1:2765898/1‑93 (MQ=255)
                                tttGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCAGCTTCTACATCGCGCGCGCACTCTTTGGCGGTc                                                        >  2:1688054/1‑93 (MQ=255)
                                 ttGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCt                                                       >  2:3451848/1‑93 (MQ=255)
                                   ggCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGc                                                     <  1:1344899/93‑1 (MQ=255)
                                     cATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAg                                                   <  2:2879406/93‑1 (MQ=255)
                                      atatCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGa                                                  <  1:2991265/93‑1 (MQ=255)
                                      atatCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGa                                                  >  2:4010024/1‑93 (MQ=255)
                                      atatCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGa                                                  >  2:2550058/1‑93 (MQ=255)
                                       tatCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGAt                                                 <  1:3965363/93‑1 (MQ=255)
                                        atCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATc                                                <  1:4364228/93‑1 (MQ=255)
                                           gCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCa                                             >  1:4868743/1‑93 (MQ=255)
                                            ccGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCAc                                            <  1:1411162/93‑1 (MQ=255)
                                            ccGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCAc                                            <  1:1587666/93‑1 (MQ=255)
                                             cGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGACCCCACt                                           >  2:365054/1‑93 (MQ=255)
                                              gTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTg                                          >  2:6076400/1‑93 (MQ=255)
                                              gTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTg                                          >  1:2753478/1‑93 (MQ=255)
                                               tACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTgc                                         >  1:596952/1‑93 (MQ=255)
                                                aCATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTgcg                                        >  2:5246826/1‑93 (MQ=255)
                                                 cATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTgcgg                                       >  1:1685543/1‑93 (MQ=255)
                                                 cATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTgcgg                                       <  2:3187404/93‑1 (MQ=255)
                                                 cATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTgcgg                                       >  1:4620059/1‑93 (MQ=255)
                                                  aTGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTgcggc                                      <  1:2661724/93‑1 (MQ=255)
                                                    ggCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGt                                    <  2:5846398/93‑1 (MQ=255)
                                                      cATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGa                                  >  2:5548488/1‑93 (MQ=255)
                                                          aagaaagaAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              >  1:3247493/1‑93 (MQ=255)
                                                          aagaaagaAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              <  2:5598514/93‑1 (MQ=255)
                                                          aagaaagaAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              >  2:92912/1‑93 (MQ=255)
                                                          aagaaagaAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              >  2:5779595/1‑93 (MQ=255)
                                                          aagaaagaAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              <  1:1013581/93‑1 (MQ=255)
                                                          aagaaagaAGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              >  2:1874732/1‑93 (MQ=255)
                                                          aagaaagaAGGGATAGGAGGCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATccc                              <  2:3478687/93‑1 (MQ=255)
                                                           agaaagaaGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCt                             <  2:4779036/93‑1 (MQ=255)
                                                           agaaagaaGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCt                             <  1:450005/93‑1 (MQ=255)
                                                           agaaagaaGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCt                             >  1:3532204/1‑93 (MQ=255)
                                                           agaaagaaGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCt                             <  2:1939128/93‑1 (MQ=255)
                                                            gaaagaaGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTg                            <  2:4220894/93‑1 (MQ=255)
                                                             aaagaaGGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGc                           >  1:6021699/1‑93 (MQ=255)
                                                              aagaagGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCg                          >  1:5102287/1‑93 (MQ=255)
                                                               agaatGGATAGGAGTCTCTCGATAGCGCCTCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGt                         <  1:5782897/93‑1 (MQ=255)
                                                                gaagGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTg                        >  1:1504087/1‑93 (MQ=255)
                                                                 aagGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTgg                       <  2:3011011/93‑1 (MQ=255)
                                                                 aagGGATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTgg                       <  2:3894626/93‑1 (MQ=255)
                                                                 aagGGATAGGAGTCACTCGATAGCGCCGCATCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTgg                       >  1:3888439/1‑93 (MQ=255)
                                                                   gggATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGc                     <  2:2048153/93‑1 (MQ=255)
                                                                    ggATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGcc                    <  2:5708852/93‑1 (MQ=255)
                                                                     gATAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCa                   <  2:3603287/93‑1 (MQ=255)
                                                                      aTAGGAGTCACTCGATAGCACCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCAt                  >  1:2170187/1‑93 (MQ=255)
                                                                       tAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATc                 <  2:3649872/93‑1 (MQ=255)
                                                                       tAGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATc                 >  2:149875/1‑93 (MQ=255)
                                                                       aaGGAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATc                 <  2:1432215/92‑1 (MQ=255)
                                                                         ggAGTCACTCGATAGCGCCGCTTCTACATCGCGGTCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCAc               <  2:5936532/93‑1 (MQ=255)
                                                                          gAGTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCAcc              >  1:5874753/1‑93 (MQ=255)
                                                                            gTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAg            >  2:2290464/1‑93 (MQ=255)
                                                                            gTCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAg            >  1:1219988/1‑93 (MQ=255)
                                                                             tCACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGc           <  1:5193220/93‑1 (MQ=255)
                                                                              cACTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCt          >  2:2574001/1‑93 (MQ=255)
                                                                               aCTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCtt         >  1:2616793/1‑93 (MQ=255)
                                                                               aCTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCtt         >  2:4827707/1‑93 (MQ=255)
                                                                                cTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCttt        <  2:5215975/93‑1 (MQ=255)
                                                                                cTCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCttt        <  1:617835/93‑1 (MQ=255)
                                                                                cTCGATAGCGCCGCTACTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCATGATCCCTGCGTGGGCCATCACCAGCttt        <  2:5496867/93‑1 (MQ=255)
                                                                                 tCGATAGCTCCTCTGCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCtttt       <  1:899585/93‑1 (MQ=255)
                                                                                 tCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCtttt       <  1:701661/93‑1 (MQ=255)
                                                                                 tCGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCtttt       >  1:44103/1‑93 (MQ=255)
                                                                                  cGATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGGCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCttttt      <  2:5056950/93‑1 (MQ=255)
                                                                                   gATAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTtc     >  1:3396466/1‑93 (MQ=255)
                                                                                     tAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTtctc   >  1:658608/1‑93 (MQ=255)
                                                                                     tAGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTtctc   >  1:3365616/1‑93 (MQ=255)
                                                                                      aGCGCCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTTctct  <  1:4965682/93‑1 (MQ=255)
                                                                                         |                                                                                         
TTCCCAGCCCCCGCAATTCCGGTGACCAGTTATTTGGCATATCGCCGTACATGGCATCAAGAAAGAAGGGATAGGAGTCACTCGATAGCACCGCTTCTACATCGCGTGCGCACTCTTTGGCGGTCTGCAGATCCCACTGCGGCGTGATCCCTGCGTGGGCCATCACCAGCTTTTTCTCT  >  NC_000913/50709‑50887

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