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.C96.0% 476.5 / 11.8 149noncoding (71/87 nt)leuTtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base C (45/98);  total (46/103)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 190 total aligned reads displayed.

CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAC  >  NC_000913/3980197‑3980376
                                                                                              |                                                                                          
 attattAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑c                                                                                          >  1:1991391/1‑93 (MQ=255)
  ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cc                                                                                         <  2:2343851/93‑1 (MQ=255)
   tattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccc                                                                                        <  1:4721836/93‑1 (MQ=255)
   tattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccc                                                                                        <  1:2709921/93‑1 (MQ=255)
      taGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGT‑‑cccccc                                                                                    >  1:3671440/1‑93 (MQ=255)
        gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑cccccccc                                                                                  >  2:2044030/1‑93 (MQ=255)
        gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGT‑‑cccccccc                                                                                  >  2:3496330/1‑93 (MQ=255)
         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccc                                                                                 <  1:4587606/93‑2 (MQ=255)
         aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGG‑‑ccccccccc                                                                                 <  2:1792406/93‑2 (MQ=255)
          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccct                                                                                <  2:3596532/93‑3 (MQ=255)
          aGTTGTGACAATGCGAAGGTGGCGAAATTGGTAGCCGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccct                                                                                <  2:768275/93‑3 (MQ=255)
           gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccctc                                                                               >  2:6040593/1‑90 (MQ=255)
           gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccctc                                                                               <  2:369150/93‑4 (MQ=255)
           gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑ccccccccctc                                                                               >  2:4948481/1‑90 (MQ=255)
                gACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGcacc                                                                           <  1:103559/93‑1 (MQ=255)
                 aCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGcacca                                                                          <  2:4023954/93‑1 (MQ=16)
                  cAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGcaccac                                                                         <  1:4835918/93‑1 (MQ=16)
                    aTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGa                                                                       <  2:3280606/93‑1 (MQ=16)
                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGAc                                                                     <  2:2717349/93‑1 (MQ=255)
                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGAc                                                                     <  2:127722/93‑1 (MQ=255)
                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGAc                                                                     <  1:1455288/93‑1 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGCCGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    <  2:5295196/93‑1 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    <  1:4895585/93‑1 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    >  2:5648186/1‑93 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    <  2:541183/93‑1 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    <  1:3103018/93‑1 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    <  1:2835946/93‑1 (MQ=255)
                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACt                                                                    >  2:5110536/1‑93 (MQ=255)
                        gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACttt                                                                  >  1:305370/1‑93 (MQ=255)
                        gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACttt                                                                  <  1:292160/93‑1 (MQ=255)
                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTa                                                                 >  2:6304388/1‑93 (MQ=255)
                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTa                                                                 <  1:3881153/93‑1 (MQ=255)
                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCCTCGCACCACGACTTTa                                                                 >  1:5742205/1‑93 (MQ=255)
                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTaa                                                                >  1:2771589/1‑93 (MQ=255)
                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCCCACCACGACTTTaa                                                                >  1:1034729/1‑93 (MQ=255)
                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTT‑‑CCCCCCCCTCTCACCACCACCTTaa                                                                >  1:387983/1‑93 (MQ=11)
                            gTGGCGGAATTGGTAGACGCGCTAGCTTCATGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAg                                                              <  1:2497738/93‑1 (MQ=255)
                             tGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGa                                                             <  2:3825348/93‑1 (MQ=255)
                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGT‑‑CCCCCCCCCCGCACCACCACCTTTaaaaaa                                                            >  1:3392239/1‑93 (MQ=11)
                               gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAt                                                           <  1:3663035/93‑1 (MQ=255)
                               gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAt                                                           <  1:503035/93‑1 (MQ=255)
                                cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAtt                                                          >  2:2692846/1‑93 (MQ=255)
                                  gAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGa                                                        >  2:6253541/1‑93 (MQ=255)
                                  gAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGa                                                        >  1:1075538/1‑93 (MQ=255)
                                    aTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAAc                                                      <  2:5033965/93‑1 (MQ=255)
                                    aTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAAc                                                      >  1:865947/1‑93 (MQ=255)
                                      tGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTa                                                    >  2:5739697/1‑93 (MQ=255)
                                      tGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTa                                                    <  1:2082713/93‑1 (MQ=255)
                                       ggTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTCAACAATTGAACtca                                                   >  1:5204959/1‑91 (MQ=17)
                                        gTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaa                                                  <  2:5882648/93‑1 (MQ=255)
                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTTAAGTCCCCCCCCCTCGCCCCACGACCTTTAAGAATTGAGCTaaaa                                                 >  1:3568776/1‑93 (MQ=17)
                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCGCCCCCTCGCACCACGACTTTAAAGAATTGAACTACAAAt                                               >  2:3017755/1‑93 (MQ=255)
                                             cgcgCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTc                                             <  2:2941601/93‑1 (MQ=255)
                                             cgcgCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTc                                             <  1:4362292/93‑1 (MQ=255)
                                              gcgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                            <  2:5141182/93‑1 (MQ=255)
                                                ccTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                          <  1:3236910/92‑1 (MQ=255)
                                                 cTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaa                                         >  2:6156614/1‑93 (MQ=255)
                                                   aGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCCCCCCCACCTTTAAAAAATTAAATAAAAAATCaaaaaa                                       >  1:3693390/1‑92 (MQ=9)
                                                   aGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                       <  2:2953444/93‑1 (MQ=255)
                                                   aGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                       >  2:3912215/1‑93 (MQ=255)
                                                    gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGc                                      <  1:2629419/93‑1 (MQ=255)
                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                     <  1:2019963/93‑1 (MQ=255)
                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                     <  2:3988646/93‑1 (MQ=255)
                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  2:1084166/93‑1 (MQ=255)
                                                      ttCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  2:1250894/93‑1 (MQ=255)
                                                      ttCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  1:73192/93‑1 (MQ=255)
                                                      tacAGATGTTAGTGGCCTTACGGGCTTGGCGGTTCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAg                                    <  1:5576633/91‑1 (MQ=25)
                                                       tCAGGTGTTAGTGTCCTTACGGGCGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAAGTGAACTAAAAATTCAAAAAGCAGt                                   <  1:127557/93‑1 (MQ=25)
                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                  <  1:6238664/93‑1 (MQ=255)
                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAACAGTAt                                 >  2:5809280/1‑93 (MQ=255)
                                                          ggTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                <  2:4149701/93‑1 (MQ=255)
                                                           gtgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAttt                               <  1:2054885/93‑1 (MQ=255)
                                                           gtgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAttt                               <  1:3305974/93‑1 (MQ=255)
                                                            tgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              <  1:495912/93‑1 (MQ=255)
                                                            tgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              <  2:4233212/93‑1 (MQ=255)
                                                            tgtTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTc                              >  1:687528/1‑93 (MQ=255)
                                                              ttAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCCCCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCgg                            <  2:883120/93‑1 (MQ=255)
                                                               tAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGc                           <  1:1761526/93‑1 (MQ=255)
                                                               tAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGc                           <  2:1464933/93‑1 (MQ=255)
                                                                   gtCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGt                       <  2:3699461/93‑1 (MQ=255)
                                                                   gtCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGt                       <  2:4481121/93‑1 (MQ=255)
                                                                    tCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTa                      <  1:2285971/93‑1 (MQ=255)
                                                                    tCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAAATAAAAATTCAAAAAGCAGTATTTCGGCGAGTa                      >  2:2681252/1‑93 (MQ=255)
                                                                      cTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgc                    <  2:222621/93‑1 (MQ=255)
                                                                      cTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgc                    <  1:3749114/93‑1 (MQ=255)
                                                                       ttACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcg                   <  2:3845142/93‑1 (MQ=255)
                                                                        tccGGACGTGGGGTTTCAAGCC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcgc                  <  1:1117065/91‑1 (MQ=38)
                                                                        tACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcgc                  <  1:5389982/93‑1 (MQ=255)
                                                                            gACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt              >  2:3637314/1‑93 (MQ=255)
                                                                             aCGTGGGGATTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCtt             <  1:4276425/93‑1 (MQ=255)
                                                                               gTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTgg           <  1:1488979/93‑1 (MQ=255)
                                                                                tGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGt          <  2:4752649/93‑1 (MQ=255)
                                                                                tGGGGGTTCAAGT‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTa         >  1:4037319/1‑93 (MQ=255)
                                                                                 gggggTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTa         >  2:3054350/1‑93 (MQ=255)
                                                                                 gggggTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTa         <  2:3092180/93‑1 (MQ=255)
                                                                                  ggggTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg        >  2:3707930/1‑93 (MQ=255)
                                                                                  ggggTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGAGCAGCTTCGTAg        >  1:2377960/1‑93 (MQ=255)
                                                                                  ggggTTCAAGCCCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTa         <  1:3087345/93‑1 (MQ=255)
                                                                                   gggTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAgc       >  2:2519025/1‑93 (MQ=255)
                                                                                        cAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAAc  <  1:3591012/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 < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: