Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,899 G→T 15.1% Q113H (CAG→CAT 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,8990GT15.1% 261.3 / 33.7 93Q113H (CAG→CATalaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base T (9/5);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

CATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCA  >  NC_000913/2407800‑2407994
                                                                                                   |                                                                                                 
catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAg                                                                                                   <  2:1240528/100‑1 (MQ=255)
 atcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                  <  2:970858/100‑1 (MQ=255)
     tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAtt                                                                                              >  2:1552713/1‑100 (MQ=255)
     tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAtt                                                                                              <  2:148348/100‑1 (MQ=255)
       cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGc                                                                                           <  2:2075674/100‑1 (MQ=255)
       cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGc                                                                                           <  2:2133186/100‑1 (MQ=255)
         gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTg                                                                                         <  1:479933/100‑1 (MQ=255)
            cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAAc                                                                                      <  1:945803/100‑1 (MQ=255)
            cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAc                                                                                      <  1:2276328/100‑1 (MQ=38)
            cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAc                                                                                      <  1:52104/100‑1 (MQ=38)
              aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                    <  2:1911959/100‑1 (MQ=255)
              aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                    >  2:1364426/1‑100 (MQ=255)
               ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGc                                                                                   >  2:1279277/1‑100 (MQ=255)
               ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGc                                                                                   >  2:2464756/1‑100 (MQ=255)
               ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAAAAGc                                                                                   <  2:807344/100‑1 (MQ=255)
                cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAcgata                                                                                  <  2:1341257/100‑5 (MQ=37)
                cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAcgata                                                                                  <  2:1341093/100‑5 (MQ=37)
                 aGGCTCGTGGCATGCGTTATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgg                                                                                 <  2:2240623/100‑1 (MQ=255)
                  ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCggg                                                                                >  1:1900313/1‑100 (MQ=255)
                  ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCggg                                                                                >  2:313978/1‑100 (MQ=255)
                    cTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                               >  1:1971616/1‑99 (MQ=255)
                     tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGAc                                                                            <  1:1072022/100‑1 (MQ=255)
                     tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGAc                                                                            <  2:1122693/100‑1 (MQ=255)
                     tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGCCTGAAcgataagctc                                                                             <  2:1366325/100‑10 (MQ=21)
                       gTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGa                                                                          >  2:945965/1‑100 (MQ=255)
                        tGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGCCGaa                                                                         >  1:6399/1‑100 (MQ=255)
                         ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGaaa                                                                        >  1:1607301/1‑100 (MQ=255)
                         ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGaa                                                                         <  2:1491842/99‑1 (MQ=255)
                           cATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATg                                                                      >  2:1936089/1‑100 (MQ=255)
                             tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                    >  2:917272/1‑100 (MQ=255)
                             tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                    <  1:5172/100‑1 (MQ=255)
                             ctcGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                    <  1:784220/98‑1 (MQ=255)
                              gagtgatgtgatCGTGGAAGATATTTACCGCGTCAATGGTGTAGCGGAGCGTAGCGTTCAGTCAACGCAGGCATTGCTGAACAGCGGGGACGAAATGTTg                                                                   <  1:2406374/98‑1 (MQ=255)
                              gCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTg                                                                   <  1:558049/100‑1 (MQ=255)
                              gCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTg                                                                   >  2:1026547/1‑100 (MQ=255)
                               cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGttag                                                                  >  2:784223/1‑98 (MQ=255)
                               cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAACGtt                                                                    >  1:207028/1‑98 (MQ=255)
                                gTGATGTTGCCTTACAAGATATGTCCATCGGCAATGTTCTATCGGAGCTTATCGTTCAGGCAATGCAGCCATTTCTGAACAGCGGTTACGAAATGTTg                                                                   <  1:1953560/98‑1 (MQ=255)
                                   aTGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTcc                                                              <  2:1069181/100‑1 (MQ=255)
                                   aTGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTcc                                                              <  1:2518092/100‑1 (MQ=255)
                                     gTTACCGTGGAAGCTATATACATCGGCAATGGTGTATCGGAGTTAATCGTTCAGGCAATGCAGGCGTTGCTGAACAGCGGGGACGAAATGTGAGTACCTg                                                            >  1:183138/1‑100 (MQ=255)
                                       tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCa                                                          >  2:1398693/1‑100 (MQ=255)
                                       tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCa                                                          <  2:2308371/100‑1 (MQ=255)
                                            tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGa                                                     <  2:1722406/100‑1 (MQ=255)
                                              gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGAtt                                                   >  2:1216741/1‑100 (MQ=255)
                                                aGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                 <  1:1040627/100‑1 (MQ=255)
                                                 gATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAcc                                                <  2:254571/100‑1 (MQ=255)
                                                  atatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAccc                                               >  1:1335615/1‑100 (MQ=255)
                                                   tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCa                                              >  1:2667/1‑100 (MQ=255)
                                                    atTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCGGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAc                                             <  1:62052/100‑1 (MQ=255)
                                                     tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGACATGTTGGTTCCTGCACCAGATTACCCAc                                             >  1:44799/1‑99 (MQ=255)
                                                     tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAct                                            >  2:687763/1‑100 (MQ=255)
                                                       tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                          >  1:2122928/1‑100 (MQ=255)
                                                       tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                          <  1:1694957/100‑1 (MQ=255)
                                                        aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTg                                         >  2:42336/1‑100 (MQ=255)
                                                         cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                        >  2:818989/1‑100 (MQ=255)
                                                         cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                        <  1:782359/100‑1 (MQ=255)
                                                           tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGAc                                      >  2:1662918/1‑100 (MQ=255)
                                                            cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGCTTACCCACTCTGGAcc                                     >  1:1515441/1‑100 (MQ=255)
                                                               cAATGGTGTCTCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                  >  1:2557664/1‑100 (MQ=255)
                                                                    gtgtATCGGAGCTTATCGTTCAGGCAATGCAGGCAGTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGcggcgg                              <  2:33848/99‑1 (MQ=255)
                                                                        aTCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                         <  2:2213647/100‑1 (MQ=255)
                                                                           ggAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                      >  1:698005/1‑100 (MQ=255)
                                                                            gtccccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     >  2:1100310/10‑100 (MQ=21)
                                                                            gtccccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     >  2:978071/10‑100 (MQ=21)
                                                                            gAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     <  2:1231143/100‑1 (MQ=255)
                                                                              gCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGAGCGGGGCGGTTTCGCTTTcc                   >  1:2139008/1‑100 (MQ=255)
                                                                              gCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                   <  1:421499/100‑1 (MQ=255)
                                                                              ccccgctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                   >  2:2281967/8‑100 (MQ=25)
                                                                               cgcggggtacAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     <  1:1466382/88‑1 (MQ=255)
                                                                                ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                 <  2:469112/100‑1 (MQ=255)
                                                                                ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                 >  2:2389104/1‑100 (MQ=255)
                                                                                 tATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                >  1:1888160/1‑100 (MQ=255)
                                                                                 tATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                >  2:773334/1‑100 (MQ=255)
                                                                                  gctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg               >  2:735361/4‑100 (MQ=37)
                                                                                  gctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg               >  2:1578196/4‑100 (MQ=37)
                                                                                  gctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg               >  1:1971589/4‑100 (MQ=37)
                                                                                   tCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCgg              <  2:123148/100‑1 (MQ=255)
                                                                                     gTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa            >  2:1671239/1‑100 (MQ=255)
                                                                                     gTTCAGGCAATGCAGGCATTGCTCACACGGGGGACCCCAATGATGGTTACTAAAACAGATTACCCACTCTGAACGGCGGCGAGTTCGGTTTACAGCCGt             >  1:483254/1‑99 (MQ=255)
                                                                                     gTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa            >  2:2315693/1‑100 (MQ=38)
                                                                                     gTTCAGCAATGCCTGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTaa           >  2:896085/1‑100 (MQ=255)
                                                                                        tcagCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGATTTCCAGCGGTaaa          >  2:2356006/4‑100 (MQ=255)
                                                                                         aGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc        <  2:758508/100‑1 (MQ=255)
                                                                                         aGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc        >  2:343516/1‑100 (MQ=255)
                                                                                         aGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc        <  1:1026544/100‑1 (MQ=255)
                                                                                          ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCg       >  2:189822/1‑100 (MQ=255)
                                                                                          ggCAATGCAGGCATTCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCg       <  2:2557683/100‑1 (MQ=255)
                                                                                            cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt     <  1:465966/100‑1 (MQ=255)
                                                                                             aaTGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg    <  1:623423/100‑1 (MQ=255)
                                                                                              aTGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc   <  2:2214550/100‑1 (MQ=255)
                                                                                              aTGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc   <  1:1189176/100‑1 (MQ=255)
                                                                                               tGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa  <  2:350988/100‑1 (MQ=255)
                                                                                                gCAGGCAATGGTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTTGCTTTACCGCTGGAAAGCgaaac   <  1:825849/98‑5 (MQ=255)
                                                                                                 cAGGCACTGCTGAACAGCGGGGGCGAAATGTTGGTTCATGCACCAGATTACACAGTTTGAACCGCGGCGATTTCGCCTTCCAGCGGTAAAGCGGTGc   >  1:1718804/1‑97 (MQ=255)
                                                                                                   |                                                                                                 
CATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCA  >  NC_000913/2407800‑2407994

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