Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,898 A→C 15.4% Q113P (CAG→CCG)  alaA → glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,407,8980AC15.4% 487.2 / 81.0 194Q113P (CAG→CCG) alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (80/84);  new base C (18/12);  total (98/96)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01

Only 100 of 197 total aligned reads displayed.

CCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTA  >  NC_000913/2407799‑2407997
                                                                                                   |                                                                                                    
ccATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCa                                                                                                      >  1:2541594/1‑100 (MQ=255)
  atcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                    <  1:2203934/100‑1 (MQ=255)
  atcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                    >  2:162706/1‑100 (MQ=255)
   tcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGc                                                                                                   >  2:3304058/1‑100 (MQ=255)
   tcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGc                                                                                                   >  1:1753339/1‑100 (MQ=255)
    catGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCa                                                                                                  <  2:2011317/100‑1 (MQ=255)
    catGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACAGCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCa                                                                                                  <  1:3298904/100‑1 (MQ=255)
      tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGCGCTTATCGTTCAGGCAATGCAGGCAtt                                                                                                >  1:4201617/1‑100 (MQ=255)
      tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAtt                                                                                                <  2:2442211/100‑1 (MQ=255)
        cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGc                                                                                             >  1:2003525/1‑100 (MQ=255)
         agcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCt                                                                                            >  2:3129539/1‑100 (MQ=255)
         agcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCt                                                                                            >  2:2838101/1‑100 (MQ=255)
         agcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCt                                                                                            >  1:3621419/1‑100 (MQ=255)
          gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTg                                                                                           >  1:1341478/1‑100 (MQ=255)
          gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGcctg                                                                                          <  1:4069986/100‑4 (MQ=255)
           gaggaCCAGGTTCGTGGCATGTGTGATGTTAGCGTGAAAGATATTTACATCGGCAATGGTGTATCGGAGCTGATCGTTCAGGCAATGCAGGCATTGCtn                                                                                           <  1:3508906/95‑2 (MQ=255)
             cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAc                                                                                        <  2:3503884/100‑1 (MQ=38)
             cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAc                                                                                        <  1:720551/100‑1 (MQ=38)
              tACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACa                                                                                       <  2:2706839/100‑1 (MQ=255)
               aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                      >  2:611024/1‑100 (MQ=255)
                 cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCCTGCATTGCCTGAAcgat                                                                                     <  1:699954/100‑4 (MQ=255)
                  aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgg                                                                                   <  1:449532/100‑1 (MQ=255)
                   ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCggg                                                                                  >  2:4069872/1‑100 (MQ=255)
                       cGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACg                                                                              <  2:1552015/100‑1 (MQ=255)
                        gTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGa                                                                             >  1:2470877/1‑100 (MQ=255)
                          ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGaaa                                                                           >  2:3424387/1‑100 (MQ=255)
                            cATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATAGTTCAGGAAATGCAGGCATTGCTGATCAGACGGGTCGAAAt                                                                          >  1:1704569/1‑99 (MQ=255)
                               gCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTAGCGGGGCTGATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTg                                                                      <  1:189244/100‑1 (MQ=255)
                                cGTGATGTTACCTTGGAAGATATTTACATAGGCAATGGTGTATCGGATCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGTGGGGACGAACTGtt                                                                       >  1:177074/1‑98 (MQ=255)
                                cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTgg                                                                     <  2:1228486/100‑1 (MQ=255)
                                cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTAGCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTgg                                                                     <  1:818211/100‑1 (MQ=255)
                                 gTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGt                                                                    <  1:2770606/100‑1 (MQ=255)
                                  tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGtt                                                                   <  2:2426238/100‑1 (MQ=255)
                                       ttACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGc                                                              >  2:538180/1‑100 (MQ=255)
                                        tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCa                                                             <  2:3235902/100‑1 (MQ=255)
                                          ccGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCAcc                                                           >  2:841946/1‑100 (MQ=255)
                                               gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGCTTCCTGCACCAGAtt                                                      >  2:4326548/1‑100 (MQ=255)
                                                aaGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                     >  2:136030/1‑100 (MQ=255)
                                                aaGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                     <  1:2143305/100‑1 (MQ=255)
                                                  gATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAcc                                                   <  2:3987238/100‑1 (MQ=255)
                                                   atatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAccc                                                  <  1:2618652/100‑1 (MQ=255)
                                                   acatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAccc                                                  <  2:2273618/98‑1 (MQ=255)
                                                    tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCa                                                 >  2:2286165/1‑100 (MQ=255)
                                                    tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCa                                                 >  1:294733/1‑100 (MQ=255)
                                                    tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCa                                                 >  2:3592368/1‑100 (MQ=255)
                                                      tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAct                                               >  1:3314734/1‑100 (MQ=255)
                                                      tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAct                                               <  2:2470889/100‑1 (MQ=255)
                                                         aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTg                                            <  2:1205328/100‑1 (MQ=255)
                                                         aCATCGGCAATGGTGTATCGGAGCCTGTCGATCAGGCAATGGCGGCATTGCTGAACAGCGGGGACGAAATGTTGGTGCCTGCACCAGATTACCCACTCTg                                            <  1:3786983/100‑1 (MQ=255)
                                                          cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                           >  1:111414/1‑100 (MQ=255)
                                                          cATCGGCAATGGTGTATCGGAGCTGATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                           <  1:4270108/100‑1 (MQ=255)
                                                              ggCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcg                                       >  2:3141253/1‑100 (MQ=255)
                                                                cAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                     >  2:3372928/1‑100 (MQ=255)
                                                                 aaTGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcgg                                    <  2:4305870/100‑1 (MQ=255)
                                                                  aTGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggc                                   <  2:207/100‑1 (MQ=255)
                                                                   tGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                  <  2:199519/100‑1 (MQ=255)
                                                                         aTCGGAGCTTATCGTTCAGGCAATGGAGGCATGGCTGAACAGTGGGCACGAGATGTTCGTTCCTGCTCCAGATTACCCACTCTGAACCGCGGCGGTTTCg                            <  1:2039788/100‑1 (MQ=255)
                                                                          tCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGc                           <  2:3167168/100‑1 (MQ=255)
                                                                          tCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGc                           <  2:172025/100‑1 (MQ=255)
                                                                           cGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCt                          <  1:3610608/100‑1 (MQ=255)
                                                                            ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                         <  1:841942/100‑1 (MQ=255)
                                                                            ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                         >  2:1610296/1‑100 (MQ=255)
                                                                            ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                         >  1:477846/1‑100 (MQ=255)
                                                                            ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                         >  1:374259/1‑100 (MQ=255)
                                                                            ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                         >  2:4378351/1‑100 (MQ=255)
                                                                               gCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                      <  1:1813936/100‑1 (MQ=255)
                                                                               ccccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                      >  2:2485569/8‑100 (MQ=25)
                                                                                cccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCa                     >  2:2162760/7‑100 (MQ=25)
                                                                                 ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                    >  2:3525402/1‑100 (MQ=255)
                                                                                 ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                    <  2:1105878/100‑1 (MQ=255)
                                                                                 ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                    >  2:1835164/1‑100 (MQ=255)
                                                                                 ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                    >  2:1138006/1‑100 (MQ=255)
                                                                                 ccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                    >  1:564398/6‑100 (MQ=25)
                                                                                  tATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                   >  1:3579726/1‑100 (MQ=255)
                                                                                  tATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                   <  2:1300854/100‑1 (MQ=255)
                                                                                  cgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                   >  2:2645454/5‑100 (MQ=37)
                                                                                   gctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                  >  1:1938995/4‑100 (MQ=37)
                                                                                   gctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                  >  1:88461/4‑100 (MQ=37)
                                                                                   gctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGCTTCGCTTTCCAGCg                  >  2:2845942/4‑100 (MQ=37)
                                                                                   aTCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                  <  1:1272707/100‑1 (MQ=255)
                                                                                    tCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCgg                 <  1:157432/100‑1 (MQ=255)
                                                                                    ctgTTCAGCAATGCCTGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGt                >  2:1771171/3‑100 (MQ=255)
                                                                                     cGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGt                <  2:140283/100‑1 (MQ=255)
                                                                                      gTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa               >  2:670199/1‑100 (MQ=255)
                                                                                      gTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa               >  2:946868/1‑100 (MQ=255)
                                                                                      gTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa               <  1:3115011/100‑1 (MQ=255)
                                                                                        ttcagCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTaa              >  1:49591/5‑100 (MQ=38)
                                                                                        ttcagCAATGCCTGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTaaa             >  2:3385287/5‑100 (MQ=255)
                                                                                        tCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTaaa             >  1:3245318/1‑100 (MQ=255)
                                                                                          aGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc           >  2:3795292/1‑100 (MQ=255)
                                                                                           ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCg          >  1:1199020/1‑100 (MQ=255)
                                                                                            gCAATGCAGGTGTTGCTGAACAGCGGGGACGAAGTGTTGGTTCCTGCACCAGATTACCCGCGCTGGACCGCGGCGGTTTCGCTTTACCGCTGGAAAGCgg         <  1:2655201/100‑1 (MQ=255)
                                                                                            gCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCGACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCgg         <  1:2408976/100‑1 (MQ=255)
                                                                                            gCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCgg         <  2:2601288/100‑1 (MQ=255)
                                                                                             cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt        >  1:3906461/1‑100 (MQ=255)
                                                                                              aaTGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg       <  1:1125272/100‑1 (MQ=255)
                                                                                               aTGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc      >  2:3339974/1‑100 (MQ=255)
                                                                                                tGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa     <  1:2474306/100‑1 (MQ=255)
                                                                                                  cAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTGCCTGCACCAGATTAACCACTCTGGACCGCGGCGGTTTCGCTTTACCGCTGGAAAGCgaaa       <  1:3480995/96‑4 (MQ=255)
                                                                                                   aGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTa  >  2:1318486/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTA  >  NC_000913/2407799‑2407997

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