Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 758,405 C→T A159A (GCC→GCT)  sucA → 2‑oxoglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,4050CT98.6% 504.4 / ‑3.8 145A159A (GCC→GCT) sucA2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (62/81);  total (62/83)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 147 total aligned reads displayed.

TCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATC  >  NC_000913/758315‑758492
                                                                                          |                                                                                       
tCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGctct                                                                                       <  2:4665488/93‑1 (MQ=255)
 cTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCtctc                                                                                      <  1:989558/93‑1 (MQ=255)
   gACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCaa                                                                                    <  1:167815/93‑1 (MQ=255)
   gACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCaa                                                                                    >  1:3692123/1‑93 (MQ=255)
   gACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCaa                                                                                    >  2:2899242/1‑93 (MQ=255)
   gACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCAAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCaa                                                                                    <  1:3236977/93‑1 (MQ=255)
    aCCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAg                                                                                   >  1:1912011/1‑93 (MQ=255)
     ccGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGc                                                                                  <  1:1576359/93‑1 (MQ=255)
        aaGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCaaa                                                                               >  1:1390482/1‑93 (MQ=255)
        aaGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCaaa                                                                               <  1:4141127/93‑1 (MQ=255)
         aGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAAc                                                                              >  1:2952340/1‑93 (MQ=255)
          gCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAAcc                                                                             >  1:3252384/1‑93 (MQ=255)
           cAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCt                                                                            <  1:1266786/93‑1 (MQ=255)
           cAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCt                                                                            <  2:260170/93‑1 (MQ=255)
            aGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTa                                                                           <  1:3633151/93‑1 (MQ=255)
             gACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTAc                                                                          >  1:2126210/1‑93 (MQ=255)
               cTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTg                                                                        <  2:4684511/93‑1 (MQ=255)
                ttCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGc                                                                       <  1:1399530/93‑1 (MQ=255)
                ttCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGc                                                                       >  1:2835588/1‑93 (MQ=255)
                 tCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCg                                                                      <  2:903764/93‑1 (MQ=255)
                 tCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCg                                                                      <  1:4749613/93‑1 (MQ=255)
                  ccAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCgg                                                                     <  1:1975382/93‑1 (MQ=255)
                   cAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGc                                                                    >  2:4465480/1‑93 (MQ=255)
                    aGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGcc                                                                   <  1:3894624/93‑1 (MQ=255)
                     ggATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGccc                                                                  >  1:314385/1‑93 (MQ=255)
                     ggAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGccc                                                                  <  1:527996/93‑1 (MQ=255)
                     ggAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGccc                                                                  >  2:1699059/1‑93 (MQ=255)
                       agaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGa                                                                >  2:4320579/1‑93 (MQ=255)
                       agaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGa                                                                >  2:1965862/1‑93 (MQ=255)
                        gaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAt                                                               <  1:1250328/93‑1 (MQ=255)
                        gaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAt                                                               >  1:1213566/1‑93 (MQ=255)
                        gaCATTCAACGTCGGTTCATTTGCAAGCGGCAAATAAACCATGAAATTCGGCGAGCTGCTGGAAGCTCTCAAGCACACCTACTGCGGCCCGAt                                                               <  2:4012107/93‑1 (MQ=255)
                         aCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAtt                                                              <  1:4518343/93‑1 (MQ=255)
                         aCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAtt                                                              <  1:2504825/93‑1 (MQ=255)
                          ccTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTg                                                             <  1:3877520/93‑1 (MQ=255)
                           cTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTgg                                                            <  1:2275462/93‑1 (MQ=255)
                            ttCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAACCCTACTGCGGCCCGATTGGt                                                           <  2:4647188/93‑1 (MQ=255)
                             tCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTg                                                          >  2:4354361/1‑93 (MQ=255)
                              cAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGc                                                         <  2:1896743/93‑1 (MQ=255)
                               aaCGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                        >  2:1762007/1‑93 (MQ=255)
                               aaCGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                        >  2:3995384/1‑93 (MQ=255)
                               aaCGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                        >  1:4554123/1‑93 (MQ=255)
                                aCGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCg                                                       <  2:3901223/93‑1 (MQ=255)
                                aCGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCg                                                       >  1:1506487/1‑93 (MQ=255)
                                  gTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAg                                                     <  1:528077/93‑1 (MQ=255)
                                   tCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGt                                                    >  2:2053451/1‑93 (MQ=255)
                                   tCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGt                                                    <  1:3002237/93‑1 (MQ=255)
                                    cGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTTCCGAGta                                                   <  2:2085946/93‑1 (MQ=255)
                                    cGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGta                                                   <  1:4182684/93‑1 (MQ=255)
                                    cGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGta                                                   >  2:2129019/1‑93 (MQ=255)
                                    cGGTTCATTTGCCAGCGGCAAAGAAACCAAGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGta                                                   >  1:4188690/1‑93 (MQ=255)
                                     ggTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtat                                                  >  1:2317986/1‑93 (MQ=255)
                                      gTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtata                                                 <  2:3450763/93‑1 (MQ=255)
                                      gTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtata                                                 <  1:1695997/93‑1 (MQ=255)
                                       ttCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                                <  2:3013114/93‑1 (MQ=255)
                                       ttCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                                >  1:3489223/1‑93 (MQ=255)
                                        tCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                               >  2:143684/1‑93 (MQ=255)
                                        tCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                               <  2:3366918/93‑1 (MQ=255)
                                         cATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGc                                              <  1:154728/93‑1 (MQ=255)
                                          aTTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGca                                             >  2:2141046/1‑93 (MQ=255)
                                           tttGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGcac                                            >  2:3936962/1‑93 (MQ=255)
                                           tttGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGcac                                            >  2:476297/1‑93 (MQ=255)
                                             tGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACATACTGCGGCCCGATTGGTGCCGAGTATATGCACAt                                          >  2:1361064/1‑93 (MQ=255)
                                              tccAGCGGCAAAGAAACCATGAACCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAtt                                         <  1:3069792/92‑1 (MQ=255)
                                              gCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAtt                                         <  1:1059938/93‑1 (MQ=255)
                                               ccAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                                        >  1:4630067/1‑93 (MQ=255)
                                               ccAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                                        >  1:612025/1‑93 (MQ=255)
                                                 aGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAcc                                      <  1:712728/93‑1 (MQ=255)
                                                     gCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCa                                  >  2:1262604/1‑93 (MQ=255)
                                                       aaaGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                                >  1:551496/1‑93 (MQ=255)
                                                        aaGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCg                               >  1:2215678/1‑93 (MQ=255)
                                                          gacaCCATGCAACTCGGCGAGCTGCTAGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                             <  2:4224048/90‑1 (MQ=255)
                                                          gAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                             >  2:2642797/1‑93 (MQ=255)
                                                            aaCCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                           <  2:4025285/93‑1 (MQ=255)
                                                            aaCCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                           >  2:2258814/1‑93 (MQ=255)
                                                                aTGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                       <  1:225858/93‑1 (MQ=255)
                                                                   aaaCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGc                    <  2:4131180/93‑1 (MQ=255)
                                                                    aaCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAGCCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCt                   <  1:2471506/93‑1 (MQ=255)
                                                                    aaCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCt                   <  2:252583/93‑1 (MQ=255)
                                                                    aaCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCt                   <  1:2081657/93‑1 (MQ=255)
                                                                     aCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg                  <  2:4514132/93‑1 (MQ=255)
                                                                      cTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg                 <  1:2473374/93‑1 (MQ=255)
                                                                      cTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg                 >  1:902071/1‑93 (MQ=255)
                                                                        cGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGAt               <  1:213709/93‑1 (MQ=255)
                                                                          gCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATcc             <  2:3670924/93‑1 (MQ=255)
                                                                          gCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATcc             >  2:1495564/1‑93 (MQ=255)
                                                                            gAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCaa           >  1:1144721/1‑93 (MQ=255)
                                                                             aGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAAc          <  2:2000371/93‑1 (MQ=255)
                                                                             aGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAAc          >  1:3512693/1‑93 (MQ=255)
                                                                             aGCTGCTGGAAGCTCTCAAGCAAAACTACGGCGGCCCGATGGGTGGCGAGCATATGCACATTACCAGCACCGAAGAAAAACGGTGGATCCAAc          >  2:2419000/1‑93 (MQ=255)
                                                                              gctgctGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACa         <  1:2757743/93‑1 (MQ=255)
                                                                               gtgctgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAg        <  1:891465/92‑1 (MQ=255)
                                                                                tgctgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGc       <  2:1174967/93‑1 (MQ=255)
                                                                                tgctgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGc       <  1:4120241/93‑1 (MQ=255)
                                                                                   tgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTa    <  2:3289977/93‑1 (MQ=255)
                                                                                   tgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTa    >  1:2562058/1‑93 (MQ=255)
                                                                                   tgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTa    <  2:2498510/93‑1 (MQ=255)
                                                                                   tgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTa    <  1:2102293/93‑1 (MQ=255)
                                                                                   tgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTa    <  1:3219285/93‑1 (MQ=255)
                                                                                     gAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATc  >  2:504449/1‑93 (MQ=255)
                                                                                          |                                                                                       
TCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATC  >  NC_000913/758315‑758492

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