Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,317,042 C→T 5.1% G61S (GGC→AGC)  trpA ← tryptophan synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,317,0420CT 5.1% 282.1 / 3.8 79G61S (GGC→AGC) trpAtryptophan synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (43/32);  new base T (2/2);  total (45/34)
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.17e-01

GATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCT  >  NC_000913/1316910‑1317175
                                                                                                                                    |                                                                                                                                     
gATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGcc                                                                                                                                >  1:1120942/1‑140 (MQ=255)
gATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGcc                                                                                                                                >  1:800923/1‑140 (MQ=255)
      cATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGTGTGGCTTTTTGAATCGTCGGGCCATCCGCAAGTggg                                                                                                                          <  1:98536/140‑1 (MQ=255)
      cATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTggg                                                                                                                          >  2:357998/1‑140 (MQ=255)
      cATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTggg                                                                                                                          <  1:236529/140‑1 (MQ=255)
      cATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTggg                                                                                                                          >  2:553102/1‑140 (MQ=255)
      cATACATCAACAGGCCAATGGGAATGGTAGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTggg                                                                                                                          <  1:368862/140‑1 (MQ=255)
       aTACATGAAAAGGCAAATGGGAATGGTCGGGTGTTTCTGGAGAATCAGTGCCAGCATTTCAAAACATTGTGGAGGAGTCACACCTGCCGCATAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCTGTGGGt                                                                                                                         <  1:1643915/140‑1 (MQ=255)
       aTACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGt                                                                                                                         <  1:1643881/140‑1 (MQ=255)
        tACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTc                                                                                                                        >  1:1342810/1‑140 (MQ=255)
         aCATCCACAGGCCAATGGGAATGGTCGGGTGTATCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCg                                                                                                                       <  2:140578/140‑1 (MQ=255)
         aCATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGACGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCg                                                                                                                       <  2:140618/140‑1 (MQ=255)
               aCAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCTATCCGCCAGTGGGTCGGAGaa                                                                                                                  <  2:1581859/139‑1 (MQ=255)
               aCAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCTATCCGCCAGTGGGTCGGAGaa                                                                                                                  <  2:1584917/139‑1 (MQ=255)
               aCAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCTATCCGCCAGTGGGTCGGAGaa                                                                                                                  >  1:1581859/1‑139 (MQ=255)
               aCAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCTATCCGCCAGTGGGTCGGAGaa                                                                                                                  >  1:1584917/1‑139 (MQ=255)
               aCAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAg                                                                                                                 >  2:332581/1‑140 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             <  2:1449102/140‑1 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             <  1:1612665/140‑1 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             <  1:306315/140‑1 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             >  2:1421652/1‑140 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             <  2:1448981/140‑1 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             >  2:1655963/1‑140 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             >  2:289621/1‑140 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             <  2:416845/140‑1 (MQ=255)
                   gCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGa                                                                                                             <  1:961027/140‑1 (MQ=255)
                     cAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATa                                                                                                           >  1:146090/1‑140 (MQ=255)
                     cAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATa                                                                                                           >  2:756992/1‑140 (MQ=255)
                             aTGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTc                                                                                                   >  2:252305/1‑140 (MQ=255)
                             aTGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTc                                                                                                   >  2:252324/1‑140 (MQ=255)
                               ggTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCa                                                                                                 <  1:1046183/140‑1 (MQ=255)
                                gTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAg                                                                                                <  2:14397/140‑1 (MQ=255)
                                gTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAg                                                                                                <  1:248999/140‑1 (MQ=255)
                                gTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAg                                                                                                <  1:1653530/140‑1 (MQ=255)
                                gTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAg                                                                                                <  1:1420680/140‑1 (MQ=255)
                                   gggTGTTTCTGGCGAATCAGTGCCAGCATATCAAATCATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAgcgc                                                                                             >  2:790937/1‑140 (MQ=255)
                                               cGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGc                                                                                 <  2:80131/140‑1 (MQ=255)
                                               cGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGc                                                                                 >  1:1084310/1‑140 (MQ=255)
                                               cGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGc                                                                                 <  1:1689489/140‑1 (MQ=255)
                                                        gCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAg                                                                        <  2:1037821/140‑1 (MQ=255)
                                                        gCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAg                                                                        >  1:498764/1‑140 (MQ=255)
                                                                 tCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCCGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGAt                                                               >  2:96007/1‑140 (MQ=255)
                                                                 tCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGAt                                                               >  1:1394269/1‑140 (MQ=255)
                                                                 tCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGAt                                                               >  1:1400489/1‑140 (MQ=255)
                                                                 tCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGAt                                                               >  1:194872/1‑140 (MQ=255)
                                                                 tCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTAGAATCGTCGGGCCATCCGCCAGTGGGTCGAAGTCGGGGATACCTAACTCCAGCGCGACAGCACCGGCTTCAATTAGCGGAGCGAt                                                               >  1:568607/1‑140 (MQ=255)
                                                                   aaaaCATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATaa                                                             <  1:1645943/140‑1 (MQ=255)
                                                                       cATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAAtttt                                                         >  2:717640/1‑140 (MQ=255)
                                                                       cATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAAtttt                                                         >  2:716300/1‑140 (MQ=255)
                                                                        aTTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCAGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTc                                                        <  1:876202/140‑1 (MQ=255)
                                                                          tgtgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCaa                                                      >  2:1285455/1‑140 (MQ=255)
                                                                          tgtgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCaa                                                      >  1:556374/1‑140 (MQ=255)
                                                                          tgtgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCaa                                                      >  1:442504/1‑140 (MQ=255)
                                                                            tgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGGATCGATAATGTTCAATg                                                    >  2:471891/1‑140 (MQ=255)
                                                                            tgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGACGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATg                                                    >  1:1386421/1‑140 (MQ=255)
                                                                             gCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGa                                                   <  2:695267/140‑1 (MQ=255)
                                                                             gCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGa                                                   <  2:695192/140‑1 (MQ=255)
                                                                             gCCGGAGTCAAACATGCAGCAAAGGCGCTCATAGTGGAGTATTGAATCGTCGGGCCATCCGCCATTGGGACGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTATCAATGa                                                   >  2:1272695/1‑140 (MQ=255)
                                                                              ccGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGAc                                                  <  1:1408235/140‑1 (MQ=255)
                                                                                  aGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCt                                              >  1:188712/1‑140 (MQ=255)
                                                                                  aGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCt                                              >  2:1181107/1‑140 (MQ=255)
                                                                                      acacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAAt                                          >  1:729842/1‑140 (MQ=255)
                                                                                        acCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGc                                        <  1:935930/140‑1 (MQ=255)
                                                                                        acCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGc                                        <  1:935696/140‑1 (MQ=255)
                                                                                        acCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGc                                        <  1:935492/140‑1 (MQ=255)
                                                                                            gCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCgg                                    >  1:1583138/1‑140 (MQ=255)
                                                                                            gCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCgg                                    >  1:1398924/1‑140 (MQ=255)
                                                                                             ccGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGa                                   >  2:254645/1‑140 (MQ=255)
                                                                                                      cgcgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGc                          <  2:1582150/140‑1 (MQ=255)
                                                                                                        cgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGt                        >  2:732626/1‑140 (MQ=255)
                                                                                                        cgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGt                        >  2:733240/1‑140 (MQ=255)
                                                                                                              gTGGCGTTGTGACTCGTCGGGACATCCGCCAGTGGGTCGGAGAAGGGGATACCTATCTCCCGCGCGTCAGCACCGGCTTCAATTAGCGGATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGaa                  <  2:1108656/140‑1 (MQ=255)
                                                                                                               tGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGaaa                 >  2:297118/1‑140 (MQ=255)
                                                                                                                         aaTCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATg       >  2:1293239/1‑140 (MQ=255)
                                                                                                                            cgtcgGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATgcgc    <  2:260665/140‑1 (MQ=255)
                                                                                                                              tcgGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCct  >  2:1430725/1‑140 (MQ=255)
                                                                                                                              tcgGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCct  >  2:1430851/1‑140 (MQ=255)
                                                                                                                              tcgGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCct  >  1:685426/1‑140 (MQ=255)
                                                                                                                              tcgGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCct  >  1:153109/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                     
GATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCT  >  NC_000913/1316910‑1317175

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