Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,173,376 A→G noncoding (44/76 nt) thrT → tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,173,3760AG99.7% 1235.3 / ‑5.5 387noncoding (44/76 nt)thrTtRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (209/176);  total (209/178)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 395 total aligned reads displayed.

CAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAT  >  NC_000913/4173287‑4173465
                                                                                         |                                                                                         
 aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTc                                                                                       <  2:2281233/93‑1 (MQ=255)
 aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTc                                                                                       >  1:3954634/1‑93 (MQ=255)
 aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTc                                                                                       >  2:795389/1‑93 (MQ=255)
 aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTc                                                                                       <  2:3371187/93‑1 (MQ=255)
 aaGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTc                                                                                       >  2:2822404/1‑93 (MQ=255)
     tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCa                                                                                   <  1:3070127/93‑1 (MQ=255)
     tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCa                                                                                   >  2:7563458/1‑93 (MQ=255)
     tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCa                                                                                   <  2:4319144/93‑1 (MQ=255)
      gatgatGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGACATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGAGGGGTCGGCAg                                                                                  >  1:3793786/1‑93 (MQ=255)
       atgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGt                                                                                 <  2:860958/93‑1 (MQ=255)
        tgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGtt                                                                                >  1:6074534/1‑93 (MQ=255)
         gatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTc                                                                               >  1:4771717/1‑93 (MQ=255)
          atgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCg                                                                              >  1:179794/1‑93 (MQ=255)
           tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGa                                                                             >  1:1231126/1‑93 (MQ=255)
           tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGa                                                                             >  1:5447185/1‑93 (MQ=255)
            gCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGaa                                                                            <  1:5065006/93‑1 (MQ=255)
            gCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGaa                                                                            <  1:5732437/93‑1 (MQ=255)
             cGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAAt                                                                           <  1:3441379/93‑1 (MQ=255)
              gggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATc                                                                          <  2:1710721/93‑1 (MQ=255)
              gggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATc                                                                          <  1:5873780/93‑1 (MQ=255)
               ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCt                                                                         >  2:1547150/1‑93 (MQ=255)
               ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCt                                                                         >  1:7231500/1‑93 (MQ=255)
                gTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTg                                                                        <  1:9639855/93‑1 (MQ=255)
                 ttCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGc                                                                       >  2:8867685/1‑93 (MQ=255)
                  tCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGcc                                                                      <  1:1169292/93‑1 (MQ=255)
                  tCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGcc                                                                      <  2:2430332/93‑1 (MQ=255)
                    gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTa                                                                    <  1:5301742/93‑1 (MQ=255)
                    gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTa                                                                    >  2:5989903/1‑93 (MQ=255)
                    gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTa                                                                    <  2:3722307/93‑1 (MQ=255)
                       tCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCa                                                                 <  1:3571158/93‑1 (MQ=255)
                       tCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCa                                                                 <  2:7883285/93‑1 (MQ=255)
                          cGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGca                                                              <  1:8709073/93‑1 (MQ=255)
                           gCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcac                                                             <  2:4776696/93‑1 (MQ=255)
                            cTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcacc                                                            <  1:9987236/93‑1 (MQ=255)
                            cTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGcacc                                                            >  2:5838507/1‑93 (MQ=255)
                               cccGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCAct                                                         >  2:9210007/1‑93 (MQ=255)
                                 cGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtc                                                       >  1:5511980/1‑93 (MQ=255)
                                  gCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtct                                                      >  2:9486164/1‑93 (MQ=255)
                                  gCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtct                                                      >  1:8865778/1‑93 (MQ=255)
                                   cTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtctt                                                     >  1:3316877/1‑93 (MQ=255)
                                   cTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCActtctt                                                     >  2:3464273/1‑93 (MQ=255)
                                    tCCTAGATGTGCTGATATAGCTCAGTTGGTAGAGCTCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCttt                                                    <  2:2384881/93‑1 (MQ=255)
                                    tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCttt                                                    >  2:351396/1‑93 (MQ=255)
                                    tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCttt                                                    >  1:7656711/1‑93 (MQ=255)
                                     ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCtttt                                                   >  1:2422432/1‑93 (MQ=255)
                                     ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCtttt                                                   <  2:7734320/93‑1 (MQ=255)
                                     ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCtttt                                                   >  2:1869145/1‑93 (MQ=255)
                                      cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTtc                                                  <  1:9812935/93‑1 (MQ=255)
                                      cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTtc                                                  <  2:7319708/93‑1 (MQ=255)
                                        aGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTtctc                                                >  2:394192/1‑93 (MQ=255)
                                           tgtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTctcctc                                             >  1:6912211/1‑93 (MQ=255)
                                            gtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCtcctcc                                            >  1:1620860/1‑93 (MQ=255)
                                            gtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCtcctcc                                            >  1:10010212/1‑93 (MQ=255)
                                               cTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTATTTTCTCCTCCCTg                                         <  2:3872967/93‑1 (MQ=255)
                                                 gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtt                                       >  1:8419816/1‑93 (MQ=255)
                                                 gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtt                                       >  1:2598547/1‑93 (MQ=255)
                                                   tataGCTCAGTTGTAAGAGCTCACCCTAGGTAAGGGTGGGGTCGGCAATTCGAATCTACCTATCAGCAACACTTCTTTTCTCCTCCCTGtttt                                     <  2:6030135/93‑1 (MQ=255)
                                                   tataGCTCAGTTGGTAGAGCTCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttt                                     <  2:9820255/93‑1 (MQ=255)
                                                   tataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttt                                     >  2:3125437/1‑93 (MQ=255)
                                                   tataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttt                                     >  2:6233727/1‑93 (MQ=255)
                                                   tataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttct                                     >  1:1453748/1‑91 (MQ=255)
                                                       gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTcc                                 >  2:522441/1‑93 (MQ=255)
                                                          cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  1:6619378/93‑1 (MQ=255)
                                                          cAGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              >  1:2303957/1‑93 (MQ=255)
                                                          cAGTTGGTAGAGCGCACCCTTGGNAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                              <  1:4233721/93‑1 (MQ=255)
                                                           aGTTGGTAGAGCGCTCCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCt                             <  1:8316774/93‑1 (MQ=255)
                                                           aGTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCt                             <  2:9319644/93‑1 (MQ=255)
                                                            gTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                            <  2:9038795/93‑1 (MQ=255)
                                                            gTTGGTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                            <  1:8250216/93‑1 (MQ=255)
                                                               ggTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGttt                         >  2:2942604/1‑93 (MQ=255)
                                                               ggTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGttt                         <  2:7195245/93‑1 (MQ=255)
                                                                gTAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTa                        <  2:872293/93‑1 (MQ=255)
                                                                 tAGAGCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAt                       >  2:7020733/1‑93 (MQ=255)
                                                                  agagCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAtt                      >  1:2897963/1‑93 (MQ=255)
                                                                  agagCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAtt                      <  1:9206782/93‑1 (MQ=255)
                                                                   gagCGCACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTg                     >  1:1347169/1‑93 (MQ=255)
                                                                     gcgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCa                   >  1:3171500/1‑93 (MQ=255)
                                                                     gcgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCa                   <  2:2710057/93‑1 (MQ=255)
                                                                     gcgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCa                   <  1:7365518/93‑1 (MQ=255)
                                                                     gcgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCa                   <  1:9081248/93‑1 (MQ=255)
                                                                     gcgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCa                   <  1:7940278/93‑1 (MQ=255)
                                                                      cgcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAt                  <  2:4615341/93‑1 (MQ=255)
                                                                       gcACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAtt                 <  2:6343728/93‑1 (MQ=255)
                                                                        cACCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATTAGGAGCACTTCTTTTCTCCTCTCTGTTTTTTCTTTCTGTTTATTGGAttt                >  1:2084326/1‑92 (MQ=37)
                                                                         aCCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca               >  2:1361932/1‑93 (MQ=255)
                                                                         aCCCTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca               >  2:5064797/1‑93 (MQ=255)
                                                                          cccTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaa              >  2:6834720/1‑93 (MQ=255)
                                                                          cccTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaa              >  2:5544543/1‑93 (MQ=255)
                                                                           ccTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaac             <  2:210219/93‑1 (MQ=255)
                                                                            cTTGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaaca            <  2:5768704/93‑1 (MQ=255)
                                                                             ttGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaacaa           >  1:5212372/1‑93 (MQ=255)
                                                                              tGGTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAg          <  2:4045991/93‑1 (MQ=255)
                                                                               ggTGAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGt         <  2:109465/93‑1 (MQ=255)
                                                                                gTAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTc        >  1:7754860/1‑93 (MQ=255)
                                                                                 tAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCg       >  2:8272287/1‑93 (MQ=255)
                                                                                 tAAGGGTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCg       <  1:1974233/93‑1 (MQ=255)
                                                                                    gggTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGc    <  2:9113113/93‑1 (MQ=255)
                                                                                    gggTGGGGTCGGCAGTTCGAATCTGCCTAACAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGc    <  1:6398216/93‑1 (MQ=255)
                                                                                     ggTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCa   <  1:8096606/93‑1 (MQ=255)
                                                                                     ggTGGGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCa   <  1:4391810/93‑1 (MQ=255)
                                                                                         |                                                                                         
CAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAT  >  NC_000913/4173287‑4173465

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