Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,297 (C)8→9 noncoding (79/87 nt) leuT → tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,2891.C89.1% 429.8 / 49.0 156noncoding (71/87 nt)leuTtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (4/13);  new base C (42/97);  total (46/110)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAC  >  NC_000913/3980197‑3980376
                                                                                               |                                                                                           
cATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGt                                                                                              <  2:4393461/93‑1 (MQ=255)
cATTATTAAAAGTTGTGACAATACGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGt                                                                                              <  2:2882033/93‑1 (MQ=255)
 attattAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑c                                                                                           >  1:1134376/1‑93 (MQ=255)
 attattAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑c                                                                                           >  1:3640039/1‑93 (MQ=255)
  ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGT‑‑cc                                                                                         >  2:4306035/1‑93 (MQ=255)
   tattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccc                                                                                        >  2:3054562/1‑93 (MQ=255)
     ttaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccc                                                                                     >  2:2052006/1‑93 (MQ=255)
     ttaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccc                                                                                     >  1:2215973/1‑93 (MQ=255)
       aGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccc                                                                                   >  2:1040558/1‑93 (MQ=255)
        gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGC‑‑cccccccc                                                                                  >  2:3663035/1‑93 (MQ=255)
         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cgccccccc                                                                                 >  1:188325/1‑92 (MQ=255)
         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccc                                                                                 <  1:633616/93‑2 (MQ=255)
          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAAT‑‑ccccccccct                                                                                >  1:3391415/1‑91 (MQ=255)
           gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccctc                                                                               >  1:4158779/1‑90 (MQ=255)
           gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTTCAGT‑‑cgccccccccc                                                                               >  2:866720/1‑90 (MQ=14)
           attGTGACAATGCGAAGGTGCCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccctc                                                                               <  1:3920314/92‑4 (MQ=255)
                  cAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCCcaccac                                                                         >  2:1849711/1‑93 (MQ=255)
                  cAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCCTCGcaccac                                                                         >  1:4023163/1‑93 (MQ=255)
                    aTGCGAAGGTGGCGGAATTGGTAGACGCTCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGa                                                                       <  1:442183/93‑1 (MQ=255)
                    aTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCCCACCACgc                                                                       >  1:1542513/1‑92 (MQ=11)
                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGAc                                                                      <  1:1585515/93‑1 (MQ=255)
                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGAc                                                                      <  2:3412818/93‑1 (MQ=255)
                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGAc                                                                      <  2:1783512/93‑1 (MQ=255)
                         aaGGTGGCGGAATGGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTa                                                                  <  2:4043615/93‑1 (MQ=255)
                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaa                                                                 >  2:3618974/1‑93 (MQ=255)
                           ggcggCGGAATGGGTAGGCGCGCCAGGTTCCGGTGTTTGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaaa                                                                <  2:78021/90‑1 (MQ=14)
                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaaa                                                                <  2:4270468/93‑1 (MQ=255)
                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaaa                                                                <  1:2155763/93‑1 (MQ=255)
                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaaa                                                                <  2:2557725/93‑1 (MQ=255)
                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaaa                                                                <  1:4199348/93‑1 (MQ=255)
                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGaa                                                            <  1:4114913/93‑1 (MQ=255)
                               gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAt                                                           >  2:3520810/1‑93 (MQ=255)
                                cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAtt                                                          <  1:224311/93‑1 (MQ=255)
                                cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAtt                                                          <  1:3663882/93‑1 (MQ=255)
                                cGGAATTGGTAGACGCGCTAGCTTCAGGAGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAtt                                                          <  1:4510495/93‑1 (MQ=255)
                                  gAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGa                                                        <  2:621531/93‑1 (MQ=255)
                                   aaTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGaa                                                       >  2:3853407/1‑93 (MQ=255)
                                   aaTTGGTAGACGCGCTAGCTTCAGGCGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTCaa                                                       >  1:4398586/1‑93 (MQ=17)
                                    aTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAAc                                                      >  2:2737507/1‑93 (MQ=255)
                                     ttGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACt                                                     <  2:2592402/93‑1 (MQ=255)
                                      tGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTa                                                    <  1:2890702/93‑1 (MQ=255)
                                       tgTAGACGCGCTAGCTTCATGTGTTAGTGTCCTGACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaa                                                   <  1:510797/92‑1 (MQ=18)
                                       ggTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaa                                                   <  1:2788687/93‑1 (MQ=255)
                                       ggTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaa                                                   <  2:3940360/93‑1 (MQ=255)
                                         tAGACGCGCTAGTTTCAGGTGTTAGTGTCCTTATGGAAGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaaa                                                 <  1:709464/93‑1 (MQ=18)
                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaa                                                  <  1:2364274/93‑1 (MQ=18)
                                            aCGCGCTTGGCTCCGGTGTTTGTGGCCTTTCCGGCCGTGGGGTTCCAGGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAt                                               <  1:1616677/93‑1 (MQ=255)
                                                gcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                          >  2:711837/1‑93 (MQ=255)
                                                 cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaa                                         <  2:3487762/93‑1 (MQ=255)
                                                 cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaa                                         >  2:3771187/1‑93 (MQ=255)
                                                 cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaa                                         <  2:647225/93‑1 (MQ=255)
                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGCATTTAAAGAATTGAACTAAAAATTCaaaaa                                        >  1:1938878/1‑93 (MQ=18)
                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                        <  1:618812/93‑1 (MQ=255)
                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                        <  2:3448405/93‑1 (MQ=255)
                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                        <  2:3703904/93‑1 (MQ=255)
                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                        <  2:2141392/93‑1 (MQ=255)
                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGAC‑‑CCCCCCCCTCGCCCCACGACCTTAAAAAAATGAAAGAAAAAATaaaaaa                                        >  2:1190695/1‑93 (MQ=9)
                                                    gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGc                                      <  2:2623251/93‑1 (MQ=255)
                                                    gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGc                                      <  1:1749733/93‑1 (MQ=255)
                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                     <  1:2740465/93‑1 (MQ=255)
                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                     <  1:1006900/93‑1 (MQ=255)
                                                      ttCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGTT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    >  2:1579911/1‑93 (MQ=255)
                                                       tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGt                                   <  1:1992194/93‑1 (MQ=255)
                                                       tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGt                                   >  2:2874284/1‑93 (MQ=255)
                                                       tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCCCACCACGACCTTAAAGAAATTAAACAAAAAATTaaaaaacagt                                    >  2:3355572/1‑89 (MQ=17)
                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                  <  2:4176875/93‑1 (MQ=255)
                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                  <  1:2949359/93‑1 (MQ=255)
                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                  <  1:4633058/93‑1 (MQ=255)
                                                         aGGTGTTAGTGTCCTTACTGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAt                                 <  2:243608/93‑1 (MQ=255)
                                                          ggTGTTAGTTTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  2:3289878/93‑1 (MQ=255)
                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                >  1:4335298/1‑93 (MQ=255)
                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  1:837898/93‑1 (MQ=255)
                                                           gtgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAttt                               <  2:250238/93‑1 (MQ=255)
                                                           gtgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCCCGCACCCCGACTTTAAAAAATTGAACTAAAAATTCAAAAAACAGTAttt                               >  2:4099013/1‑93 (MQ=21)
                                                            tgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              <  1:3188882/93‑1 (MQ=255)
                                                            tgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              <  2:2028932/93‑1 (MQ=255)
                                                            tgtTAGGGTCCTTACGGACGTGGGGGTTCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              <  2:4451718/93‑1 (MQ=37)
                                                             gtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                             <  1:3493806/93‑1 (MQ=255)
                                                              ttAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCgg                            >  1:2101058/1‑93 (MQ=255)
                                                                    tCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTa                      <  2:2226140/93‑1 (MQ=255)
                                                                    tCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTa                      >  2:1265822/1‑93 (MQ=255)
                                                                      cTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgc                    <  2:4241163/93‑1 (MQ=255)
                                                                      cTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgc                    <  2:496947/93‑1 (MQ=255)
                                                                         tACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcgc                  >  1:4542933/1‑93 (MQ=255)
                                                                           cGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAg                <  2:1521722/93‑1 (MQ=255)
                                                                            ggACGTAGGGGTAAAAGT‑‑CCCCTCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGc               <  1:2167368/93‑1 (MQ=255)
                                                                            ggACGGGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGc               <  1:349720/93‑1 (MQ=255)
                                                                             gACGTGGGGGTTCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              <  1:2471170/93‑1 (MQ=255)
                                                                             gACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              >  2:1999677/1‑93 (MQ=255)
                                                                             gACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              >  2:2313599/1‑93 (MQ=255)
                                                                             gACGTGGGGAGTCAAGC‑‑CACCCCCCCTCGCACCACGACTATAAAGAATTGAACTAAAAATTCAAAAAGCACTATTTCGGCGAGTAGCGCAGCt              <  2:1974481/93‑1 (MQ=255)
                                                                              aCGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCtt             >  1:4663774/1‑93 (MQ=255)
                                                                                gTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTgg           <  1:3042561/93‑1 (MQ=255)
                                                                                    gggTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAgc       <  2:270553/93‑1 (MQ=255)
                                                                                      gTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAgcgc     <  2:1736983/93‑1 (MQ=255)
                                                                                       ttCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCaa   >  2:3167783/1‑93 (MQ=255)
                                                                                       ttCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGGAGCGCaa   >  1:4265354/1‑93 (MQ=255)
                                                                                         cAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc  <  2:1460657/93‑1 (MQ=255)
                                                                                         cAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc  <  1:2987813/93‑1 (MQ=255)
                                                                                         cAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc  <  1:2592241/93‑1 (MQ=255)
                                                                                               |                                                                                           
CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAC  >  NC_000913/3980197‑3980376

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