Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,710,555 (T)8→7 coding (9/699 nt) b3548 ← putative lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,710,5480T.99.6% 1256.2 / ‑3.4 272coding (16/699 nt)b3548putative lipase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (115/156);  total (115/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 317 total aligned reads displayed.

GATATATTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCATAGTAGC  >  NC_000913/3710456‑3710636
                                                                                            |                                                                                        
  tatatTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCttt                                                                                        >  1:44305/1‑93 (MQ=255)
  tatatTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCttt                                                                                        >  2:5305698/1‑93 (MQ=255)
   atatTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCtttt                                                                                       <  1:1747048/93‑1 (MQ=255)
    tatTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCttttt                                                                                      >  2:2802094/1‑93 (MQ=255)
    tatTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCttttt                                                                                      <  1:1832046/93‑1 (MQ=255)
      ttCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCttttttt                                                                                    >  2:3154159/1‑93 (MQ=255)
      ttCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCttttttt                                                                                    >  2:5563247/1‑93 (MQ=255)
        cTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCtttttttat                                                                                  <  1:2350415/93‑3 (MQ=255)
         ttgttgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCtttttttatg                                                                                 >  1:4539032/1‑90 (MQ=255)
         ttgttgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCtttttttatg                                                                                 >  2:3538248/1‑90 (MQ=255)
           gttgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGAt                                                                              <  2:4287595/93‑1 (MQ=255)
             tgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCa                                                                            <  1:2071478/93‑1 (MQ=255)
             tgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCa                                                                            >  2:2997678/1‑93 (MQ=255)
             tgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCa                                                                            >  1:3241002/1‑93 (MQ=255)
             tgCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCa                                                                            <  1:4723668/93‑1 (MQ=255)
              gCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCAt                                                                           <  2:1670273/93‑1 (MQ=255)
              gCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCAt                                                                           <  2:5757848/93‑1 (MQ=255)
              gCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCAt                                                                           >  2:2407825/1‑93 (MQ=255)
                 aGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTaa                                                                        >  2:1200767/1‑93 (MQ=255)
                   gCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAaca                                                                      <  2:1754337/93‑1 (MQ=255)
                    cGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAacac                                                                     <  1:719304/93‑1 (MQ=255)
                    cGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAacac                                                                     <  1:5608858/93‑1 (MQ=255)
                    cGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAacac                                                                     <  2:1863147/93‑1 (MQ=255)
                       aaGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGAc                                                                  <  1:1182839/93‑1 (MQ=255)
                       aaGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGAc                                                                  <  2:5071035/93‑1 (MQ=255)
                        aGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACt                                                                 >  1:5307036/1‑93 (MQ=255)
                        aGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACt                                                                 <  2:6006520/93‑1 (MQ=255)
                          cTGTGTTGAGAAAAAAGGCACTGATTCCCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTcc                                                               <  2:3425964/93‑1 (MQ=255)
                          cTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTcc                                                               >  1:5481704/1‑93 (MQ=255)
                          cTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTcc                                                               <  1:4376686/93‑1 (MQ=255)
                          cTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTcc                                                               >  2:776696/1‑93 (MQ=255)
                            gtgtTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAg                                                             <  2:738644/93‑1 (MQ=255)
                              tttGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAggg                                                           <  2:267786/92‑1 (MQ=255)
                               ttGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGa                                                          >  1:3636714/1‑93 (MQ=255)
                                tGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGaa                                                         <  1:502477/93‑1 (MQ=255)
                                 taaaaaaaaGGCACTGGTTACCACGGTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAag                                                        <  1:4023124/90‑1 (MQ=255)
                                 gagaAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAag                                                        <  1:2540374/93‑1 (MQ=255)
                                  agaAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTCACACGACTCCAGGGAaga                                                       >  1:5778867/1‑93 (MQ=255)
                                  agaAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAaga                                                       >  2:4645378/1‑93 (MQ=255)
                                   gaAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAagag                                                      <  2:6680725/93‑1 (MQ=255)
                                    aaaaaaGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGc                                                     <  1:665112/93‑1 (MQ=255)
                                     aaaaaGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGcc                                                    >  2:5407229/1‑93 (MQ=255)
                                     aaaaaGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGcc                                                    >  2:961357/1‑93 (MQ=255)
                                     aaaaaGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGcc                                                    >  2:3887730/1‑93 (MQ=255)
                                       aaaGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGa                                                  >  1:217617/1‑93 (MQ=255)
                                        aaGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGaa                                                 >  1:143685/1‑93 (MQ=255)
                                         aGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAAt                                                <  1:6081399/93‑1 (MQ=255)
                                           gCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTc                                              <  2:5041113/93‑1 (MQ=255)
                                            cACTGATTCCCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCg                                             <  2:4757822/93‑1 (MQ=255)
                                             aCTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCCGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCgg                                            >  2:1089000/1‑93 (MQ=255)
                                             aCTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCgg                                            >  1:3815676/1‑93 (MQ=255)
                                             aCTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCgg                                            >  1:5363162/1‑93 (MQ=255)
                                              ctcaTTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGc                                           <  1:593926/90‑1 (MQ=255)
                                               tGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCa                                          <  1:2282556/93‑1 (MQ=255)
                                                gATTACCACGCTCTCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCaa                                         <  2:5433432/93‑1 (MQ=255)
                                                gATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGCCCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGACGAATTCGGCaa                                         <  1:669596/93‑1 (MQ=255)
                                                gATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCaa                                         >  1:3574928/1‑93 (MQ=255)
                                                gATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCaa                                         >  2:4573368/1‑93 (MQ=255)
                                                gATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTAGGCaa                                         >  1:6167411/1‑93 (MQ=255)
                                                  ttACCCCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAata                                       <  2:4676097/93‑1 (MQ=255)
                                                  ttACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAata                                       <  1:3794056/93‑1 (MQ=255)
                                                     ccACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAATACGACTCCAGGGAAGAGCCGAATTCGGCAAtatta                                    <  2:1371169/93‑1 (MQ=255)
                                                     ccACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAAtatta                                    >  1:1915680/1‑93 (MQ=255)
                                                      cACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAAtattat                                   >  1:3473015/1‑93 (MQ=255)
                                                       ccGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                                  <  1:5200181/92‑1 (MQ=255)
                                                       aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                                  <  1:2121691/93‑1 (MQ=255)
                                                       aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                                  >  1:4265422/1‑93 (MQ=255)
                                                       aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                                  >  2:6368923/1‑93 (MQ=255)
                                                       aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                                  >  1:5842892/1‑93 (MQ=255)
                                                        cGCTCGCCAGCAGGCTTAGCTGCCAGCGCTCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTg                                 <  2:2129899/93‑1 (MQ=255)
                                                        cGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTg                                 <  1:1135193/93‑1 (MQ=255)
                                                        cGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTg                                 >  2:2164175/1‑93 (MQ=255)
                                                         gCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGt                                <  2:3854762/93‑1 (MQ=255)
                                                         gCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGt                                <  1:1696726/93‑1 (MQ=255)
                                                          cTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTc                               <  2:3360458/93‑1 (MQ=255)
                                                          cTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTc                               >  1:769337/1‑93 (MQ=255)
                                                             gCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCAtt                            <  1:1542097/93‑1 (MQ=255)
                                                             gCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCAtt                            <  2:671602/93‑1 (MQ=255)
                                                              ccAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTg                           <  2:5800417/93‑1 (MQ=255)
                                                                agcagGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTa                         >  2:3311822/1‑93 (MQ=255)
                                                                agcagGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTa                         >  2:6242137/1‑93 (MQ=255)
                                                                 gcagGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTAt                        >  1:2098288/1‑93 (MQ=255)
                                                                 gcagGCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTAt                        >  2:570265/1‑93 (MQ=255)
                                                                    ggCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGaa                     <  1:4797215/93‑1 (MQ=255)
                                                                    ggCTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGaa                     <  2:2837326/93‑1 (MQ=255)
                                                                      cTTAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAgg                   <  2:1603665/93‑1 (MQ=255)
                                                                       ttAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGa                  <  2:45138/93‑1 (MQ=255)
                                                                        tAGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGat                 <  1:5074125/93‑1 (MQ=255)
                                                                         aGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGata                <  1:297834/93‑1 (MQ=255)
                                                                         aGCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGata                <  2:4107189/93‑1 (MQ=255)
                                                                          gCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGatat               <  1:612354/93‑1 (MQ=255)
                                                                          gCTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGatat               <  1:1360352/93‑1 (MQ=255)
                                                                           cTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATc              <  1:2786164/93‑1 (MQ=255)
                                                                           cTGCCAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATc              <  1:4512040/93‑1 (MQ=255)
                                                                              ccAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCggg           <  1:3143830/93‑1 (MQ=255)
                                                                              ccAGCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCggg           <  1:633755/93‑1 (MQ=255)
                                                                                 gCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAt        <  1:5752522/93‑1 (MQ=255)
                                                                                 gCGACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAt        >  2:1520951/1‑93 (MQ=255)
                                                                                   gACCACCGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAtag      <  1:1016151/93‑1 (MQ=255)
                                                                                       accGCTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCATAGTAGc  >  2:4417184/1‑93 (MQ=255)
                                                                                            |                                                                                        
GATATATTCTTGTTGCCAGGCGTAAGCTGTGTTGAGAAAAAAGGCACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCATAGTAGC  >  NC_000913/3710456‑3710636

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