Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,262,356 A→G K230E (AAG→GAG)  dusA → tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,3560AG98.6% 250.1 / ‑4.8 71K230E (AAG→GAG) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (33/37);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGA  >  NC_000913/4262222‑4262489
                                                                                                                                      |                                                                                                                                     
tcCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCg                                                                                                                                  >  1:3395/1‑140 (MQ=255)
 ccACGCAGGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGc                                                                                                                                 <  1:420503/140‑1 (MQ=255)
 ccACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGt                                                                                                                                    >  1:437735/1‑137 (MQ=255)
 ccACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGt                                                                                                                                    <  2:437735/137‑1 (MQ=255)
   acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTg                                                                                                                               >  1:94440/1‑140 (MQ=255)
   acgcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTg                                                                                                                               >  1:414255/1‑140 (MQ=255)
     gcacgTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGa                                                                                                                             >  2:353709/1‑140 (MQ=255)
          tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAgg                                                                                                                        >  1:98261/1‑140 (MQ=255)
                 tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGaa                                                                                                                            >  1:537617/1‑129 (MQ=255)
                 tGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGaa                                                                                                                            <  2:537617/129‑1 (MQ=255)
                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGcac                                                                                                               >  2:711907/1‑140 (MQ=255)
                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGcac                                                                                                               >  2:124444/1‑140 (MQ=255)
                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGca                                                                                                                <  1:542634/139‑1 (MQ=255)
                   gCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGca                                                                                                                >  2:542634/1‑139 (MQ=255)
                         gggtTTAACCCCGAAAGAAAACCGTAAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                         <  1:244760/140‑1 (MQ=255)
                         ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACTTGc                                                                                                         <  1:554549/140‑1 (MQ=255)
                         ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                         <  1:124444/140‑1 (MQ=255)
                         ggggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGc                                                                                                         <  2:212149/140‑1 (MQ=255)
                             ttAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACa                                                                                                     <  2:220048/140‑1 (MQ=255)
                                 gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                                                                                 <  2:3395/140‑1 (MQ=255)
                                   ccGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga                                                                                               <  1:458462/140‑1 (MQ=255)
                                   ccGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga                                                                                               <  1:286218/140‑1 (MQ=255)
                                    cGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggat                                                                                              <  2:414255/140‑1 (MQ=255)
                                     gaaagaaaACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggatg                                                                                             >  2:297306/1‑140 (MQ=255)
                                     gaaagaaaACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggatg                                                                                             >  1:616346/1‑140 (MQ=255)
                                      aaagaaaACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatggatgg                                                                                            >  2:628494/1‑140 (MQ=255)
                                            aaCCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                                                                                 >  2:864803/1‑129 (MQ=255)
                                            aaCCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatg                                                                                                 <  1:864803/129‑1 (MQ=255)
                                                 tGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCg                                                                                 >  2:30916/1‑140 (MQ=255)
                                                 tGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCg                                                                                 >  1:556498/1‑140 (MQ=255)
                                                 tGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCg                                                                                 >  2:628960/1‑140 (MQ=255)
                                                 tGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCg                                                                                 >  2:177716/1‑140 (MQ=255)
                                                    aaTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGgc                                                                              >  2:488575/1‑140 (MQ=255)
                                                     aTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGgcg                                                                             <  2:98261/140‑1 (MQ=255)
                                                     aTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGgcg                                                                             <  2:924529/140‑1 (MQ=255)
                                                     aTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGgcg                                                                             <  1:297306/140‑1 (MQ=255)
                                                     aTCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGgcg                                                                             <  1:51695/140‑1 (MQ=255)
                                                       cccGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGcgcg                                                                           <  1:305723/140‑1 (MQ=255)
                                                       cccGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGcgcg                                                                           <  1:488575/140‑1 (MQ=255)
                                                       cccGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGcgcg                                                                           <  1:521146/140‑1 (MQ=255)
                                                            cgcTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCg                                                                      <  1:627307/140‑1 (MQ=255)
                                                               tCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTAt                                                                   <  2:556498/140‑1 (MQ=255)
                                                                    tATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGaa                                                              <  2:282071/140‑1 (MQ=255)
                                                                              tgtATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTAtt                                                    >  1:531904/1‑140 (MQ=255)
                                                                                tATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCt                                                  >  1:652104/1‑140 (MQ=255)
                                                                                tATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCt                                                  >  1:112413/1‑140 (MQ=255)
                                                                                tATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCt                                                  >  1:872702/1‑140 (MQ=255)
                                                                                     aCTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTggcgg                                             >  1:855487/1‑140 (MQ=255)
                                                                                        gAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGcggcgg                                          <  1:628494/140‑1 (MQ=255)
                                                                                        gAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGcggcgg                                          <  1:615999/140‑1 (MQ=255)
                                                                                        gAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGcggcgg                                          >  2:314158/1‑140 (MQ=255)
                                                                                         aaGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGt                                         >  2:48196/1‑140 (MQ=255)
                                                                                          aGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTgg                                                <  1:498776/132‑1 (MQ=255)
                                                                                          aGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTgg                                                >  2:498776/1‑132 (MQ=255)
                                                                                              tGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGAcc                                    >  1:332562/1‑140 (MQ=255)
                                                                                                  tttCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGga                                <  2:94440/140‑1 (MQ=255)
                                                                                                  tttCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGga                                <  1:30916/140‑1 (MQ=255)
                                                                                                  tttCCGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGga                                <  1:177716/140‑1 (MQ=255)
                                                                                                     ccGCATCTGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGAt                             >  1:899923/1‑140 (MQ=255)
                                                                                                            tGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGt                      <  1:705299/140‑1 (MQ=255)
                                                                                                            tGACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGt                      <  1:314158/140‑1 (MQ=255)
                                                                                                             gACAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGtt                     <  1:628960/140‑1 (MQ=255)
                                                                                                               cAATGTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTcc                   <  1:353709/140‑1 (MQ=255)
                                                                                                                   gTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGt                      <  1:420750/133‑1 (MQ=255)
                                                                                                                   gTCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGt                      >  2:420750/1‑133 (MQ=255)
                                                                                                                    tCGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGa              <  2:652104/140‑1 (MQ=255)
                                                                                                                     cGATTAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGAt             >  2:901570/1‑140 (MQ=255)
                                                                                                                        ttAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGt                      >  2:695100/1‑128 (MQ=255)
                                                                                                                        ttAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGt                      <  1:695100/128‑1 (MQ=255)
                                                                                                                        ttAACGGTGGTATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATAcc          <  2:332562/140‑1 (MQ=255)
                                                                                                                                ggtATCGAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGa  >  1:293900/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTGATGGTCGGGCGCGAGGCGTATCAGAATCCGGGTATTCTGGCGGCGGTAGACCGGGAGATCTTTGGTTCCTCGGATACCGATGCCGA  >  NC_000913/4262222‑4262489

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