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,4050CT99.3% 507.2 / ‑5.2 142A159A (GCC→GCT) sucA2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (59/81);  total (60/82)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 149 total aligned reads displayed.

TCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGA  >  NC_000913/758315‑758494
                                                                                          |                                                                                         
 cTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCtctc                                                                                        <  1:316455/93‑1 (MQ=255)
     ccGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGc                                                                                    <  2:2452191/93‑1 (MQ=255)
     ccGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGc                                                                                    >  1:2994397/1‑93 (MQ=255)
      cGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCa                                                                                   <  2:510163/93‑1 (MQ=255)
      cGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCa                                                                                   <  2:3799733/93‑1 (MQ=255)
       gAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCaa                                                                                  >  1:4576093/1‑93 (MQ=255)
       gAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCaa                                                                                  >  1:3680785/1‑93 (MQ=255)
        aaGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCaaa                                                                                 <  1:1990857/93‑1 (MQ=255)
         aGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAAc                                                                                >  2:3730013/1‑93 (MQ=255)
          gCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAAcc                                                                               >  2:2284438/1‑93 (MQ=255)
           cAGACTTCCAGGACACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCt                                                                              >  2:4502860/1‑93 (MQ=255)
            aGACTTCCAGGAGACCTTCAACGGCGGGTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTa                                                                             >  1:3323008/1‑93 (MQ=255)
             gACTTCCAGGAGACCTTCCACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTAc                                                                            >  1:218791/1‑93 (MQ=255)
             gACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTAc                                                                            <  2:885884/93‑1 (MQ=255)
              aCTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACt                                                                           <  2:328162/93‑1 (MQ=255)
              aCTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACt                                                                           >  1:3061585/1‑93 (MQ=255)
              aCTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACt                                                                           <  1:2003520/93‑1 (MQ=255)
              aCTTCCAGGAGACCTTCAACGTCGGTTCATCTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACt                                                                           <  1:1914224/93‑1 (MQ=255)
               cTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTg                                                                          >  2:2868709/1‑93 (MQ=255)
               cTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTg                                                                          <  1:3907592/93‑1 (MQ=255)
                ttCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGc                                                                         >  1:1597107/1‑93 (MQ=255)
                   cAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGc                                                                      >  1:871369/1‑93 (MQ=255)
                     ggAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGccc                                                                    >  1:13037/1‑93 (MQ=255)
                      gagaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCg                                                                   >  1:1287067/1‑93 (MQ=255)
                       agaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGa                                                                  <  2:1368981/93‑1 (MQ=255)
                       agaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGACCGa                                                                  >  2:1919055/1‑93 (MQ=255)
                        gaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAt                                                                 <  2:3954996/93‑1 (MQ=255)
                        gaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAt                                                                 >  2:180872/1‑93 (MQ=255)
                        gaCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAt                                                                 >  1:1174445/1‑93 (MQ=255)
                         aCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAtt                                                                >  2:2630845/1‑93 (MQ=255)
                         aCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAtt                                                                <  2:1373024/93‑1 (MQ=255)
                         aCCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGAtt                                                                <  1:2046461/93‑1 (MQ=255)
                          ccTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTg                                                               >  1:3786516/1‑93 (MQ=255)
                           cTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTgg                                                              <  2:1757028/93‑1 (MQ=255)
                            ttCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGt                                                             <  2:4336955/93‑1 (MQ=255)
                            ttCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGt                                                             >  2:4357959/1‑93 (MQ=255)
                              cAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGc                                                           <  1:2007300/93‑1 (MQ=255)
                               aaCGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGcc                                                          >  1:689057/1‑93 (MQ=255)
                                 cGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                                        <  1:2586731/93‑1 (MQ=255)
                                 cGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                                        <  1:3396845/93‑1 (MQ=255)
                                 cGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                                        <  2:1285172/93‑1 (MQ=255)
                                  gTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAg                                                       <  2:3011891/93‑1 (MQ=255)
                                  gTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAg                                                       <  2:726617/93‑1 (MQ=255)
                                      gTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtata                                                   <  1:489248/93‑1 (MQ=255)
                                      gTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtata                                                   >  1:4593924/1‑93 (MQ=255)
                                      gTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtata                                                   <  2:1639793/93‑1 (MQ=255)
                                      gTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGtata                                                   <  2:3854758/93‑1 (MQ=255)
                                        tCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                                 <  2:3415163/93‑1 (MQ=255)
                                            ttGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGATTATATGcaca                                             <  2:2869156/93‑1 (MQ=255)
                                             tGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAt                                            <  1:3961883/93‑1 (MQ=255)
                                              gCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACAtt                                           <  2:4536546/93‑1 (MQ=255)
                                                cAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                                         <  2:1019209/93‑1 (MQ=255)
                                                  gCGGCAAAGAAACCATGAAACTCAGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                                       <  1:832485/93‑1 (MQ=255)
                                                    ggCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGc                                     <  2:1942481/93‑1 (MQ=255)
                                                     gCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCa                                    >  1:381071/1‑93 (MQ=255)
                                                     gCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCa                                    <  2:2816372/93‑1 (MQ=255)
                                                     gCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCa                                    <  2:1130614/93‑1 (MQ=255)
                                                      cAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAc                                   <  2:1994505/93‑1 (MQ=255)
                                                       aaaaaaACCATGAAACTCGGCGAGCTGCTGGAAGCACTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                                  <  2:2041865/93‑1 (MQ=255)
                                                       aaaGAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                                  >  1:3808175/1‑93 (MQ=255)
                                                          gAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                               >  1:3467359/1‑93 (MQ=255)
                                                          gAAACCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                               >  1:2099791/1‑93 (MQ=255)
                                                          aaaaCCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                               <  2:2604874/92‑1 (MQ=255)
                                                           aaaCCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                              <  2:3527290/93‑1 (MQ=255)
                                                             aCCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa                            >  2:4407554/1‑93 (MQ=255)
                                                             aCCATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa                            >  2:1603693/1‑93 (MQ=255)
                                                              ccATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaa                           <  2:401627/93‑1 (MQ=255)
                                                               cATGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa                          >  1:1724967/1‑93 (MQ=255)
                                                                aTGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                         <  2:839889/93‑1 (MQ=255)
                                                                aTGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                         <  1:3070802/93‑1 (MQ=255)
                                                                aTGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                         >  1:530589/1‑93 (MQ=255)
                                                                aTGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa                         >  1:3699395/1‑93 (MQ=255)
                                                                 tGAAACTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAAc                        <  1:2770730/93‑1 (MQ=255)
                                                                   aaaCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGc                      <  1:4195794/93‑1 (MQ=255)
                                                                     aCTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg                    >  2:87025/1‑93 (MQ=255)
                                                                      cTCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg                   <  1:3414262/93‑1 (MQ=255)
                                                                       tCGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGa                  >  2:1312945/1‑93 (MQ=255)
                                                                       ccgccGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGa                  <  2:1782436/89‑1 (MQ=255)
                                                                        cGGCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGAt                 <  1:2601768/93‑1 (MQ=255)
                                                                         ggCGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc                <  2:3693172/93‑1 (MQ=255)
                                                                           cGAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCa              <  2:3031882/93‑1 (MQ=255)
                                                                            gAGCTGCTGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCaa             >  2:4581811/1‑93 (MQ=255)
                                                                              gctgctGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACa           >  2:437854/1‑93 (MQ=255)
                                                                              gctgctGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACa           >  1:378293/1‑93 (MQ=255)
                                                                              gctgctGGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACa           <  2:2165151/93‑1 (MQ=255)
                                                                               ctgctgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAg          >  1:111282/1‑93 (MQ=255)
                                                                               ctgctgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAg          <  2:2296955/93‑1 (MQ=255)
                                                                                tgctgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGc         >  2:579911/1‑93 (MQ=255)
                                                                                   tgGAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTa      <  2:257103/93‑1 (MQ=255)
                                                                                    ggAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTAt     >  2:4259392/1‑93 (MQ=255)
                                                                                    ggAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTAt     <  2:2477078/93‑1 (MQ=255)
                                                                                    ggAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTAt     >  2:4673100/1‑93 (MQ=255)
                                                                                    ggAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTAt     <  2:2000707/93‑1 (MQ=255)
                                                                                    ggAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTAt     >  1:1048379/1‑93 (MQ=255)
                                                                                    ggAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTAt     >  2:1155070/1‑93 (MQ=255)
                                                                                     gAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATc    <  2:878828/93‑1 (MQ=255)
                                                                                     gAAGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATc    <  1:4115966/93‑1 (MQ=255)
                                                                                      aaGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCg   >  2:3228983/1‑93 (MQ=255)
                                                                                       aGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGa  <  1:2318917/93‑1 (MQ=255)
                                                                                       aGCTCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGa  >  1:3467830/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGA  >  NC_000913/758315‑758494

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