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→CAT) alaAglutamate‑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

CATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCA  >  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)
       cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑c                                                                                           <  2:2075674/100‑1 (MQ=255)
       cagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑c                                                                                           <  2:2133186/100‑1 (MQ=255)
         gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTg                                                                                         <  1:479933/100‑1 (MQ=255)
            cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAAc                                                                                      <  1:945803/100‑1 (MQ=255)
            cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAc                                                                                      <  1:2276328/100‑1 (MQ=38)
            cTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAc                                                                                      <  1:52104/100‑1 (MQ=38)
              aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAg                                                                                    <  2:1911959/100‑1 (MQ=255)
              aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAg                                                                                    >  2:1364426/1‑100 (MQ=255)
               ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGc                                                                                   >  2:1279277/1‑100 (MQ=255)
               ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGc                                                                                   >  2:2464756/1‑100 (MQ=255)
               ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAAAAGc                                                                                   <  2:807344/100‑1 (MQ=255)
                cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgata                                                                                  <  2:1341257/100‑5 (MQ=37)
                cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgata                                                                                  <  2:1341093/100‑5 (MQ=37)
                 aGGCTCGTGGCATGCGTTATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                 <  2:2240623/100‑1 (MQ=255)
                  ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCggg                                                                                >  1:1900313/1‑100 (MQ=255)
                  ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCggg                                                                                >  2:313978/1‑100 (MQ=255)
                    cTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgggg                                                                               >  1:1971616/1‑99 (MQ=255)
                     tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑c                                                                            <  1:1072022/100‑1 (MQ=255)
                     tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑c                                                                            <  2:1122693/100‑1 (MQ=255)
                     tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgataagctc                                                                             <  2:1366325/100‑10 (MQ=21)
                       gTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGa                                                                          >  2:945965/1‑100 (MQ=255)
                        tGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGC‑CGaa                                                                         >  1:6399/1‑100 (MQ=255)
                         ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGaaa                                                                        >  1:1607301/1‑100 (MQ=255)
                         ggCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGaa                                                                         <  2:1491842/99‑1 (MQ=255)
                           cATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATg                                                                      >  2:1936089/1‑100 (MQ=255)
                             tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGtt                                                                    >  2:917272/1‑100 (MQ=255)
                             tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGtt                                                                    <  1:5172/100‑1 (MQ=255)
                             ctcGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGtt                                                                    <  1:784220/98‑1 (MQ=255)
                              gagtgatgtgatCGTGGAAGATATTTACCGCGTCAATGGTGTAGCGGAGCGTAGCGTTCAGTCAACGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTg                                                                   <  1:2406374/98‑1 (MQ=255)
                              gCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTg                                                                   <  1:558049/100‑1 (MQ=255)
                              gCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTg                                                                   >  2:1026547/1‑100 (MQ=255)
                               cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGttag                                                                  >  2:784223/1‑98 (MQ=255)
                               cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAACGtt                                                                    >  1:207028/1‑98 (MQ=255)
                                gTGATGTTGCCTTACAAGATATGTCCATCGGCAATGTTCTATCGGAGCTTATCGTTCAGGCAATGCAGCCATTT‑CTGAACAGCGGTTA‑CGAAATGTTg                                                                   <  1:1953560/98‑1 (MQ=255)
                                   aTGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTcc                                                              <  2:1069181/100‑1 (MQ=255)
                                   aTGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTcc                                                              <  1:2518092/100‑1 (MQ=255)
                                     gTTACCGTGGAAGCTATATACATCGGCAATGGTGTATCGGAGTTAATCGTTCAGGCAATGCAGGCGTTG‑CTGAACAGCGGGGA‑CGAAATGTGAGTACCTg                                                            >  1:183138/1‑100 (MQ=255)
                                       tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCa                                                          >  2:1398693/1‑100 (MQ=255)
                                       tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCa                                                          <  2:2308371/100‑1 (MQ=255)
                                            tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGa                                                     <  2:1722406/100‑1 (MQ=255)
                                              gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGAtt                                                   >  2:1216741/1‑100 (MQ=255)
                                                aGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTAc                                                 <  1:1040627/100‑1 (MQ=255)
                                                 gATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTAcc                                                <  2:254571/100‑1 (MQ=255)
                                                  atatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTAccc                                               >  1:1335615/1‑100 (MQ=255)
                                                   tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCa                                              >  1:2667/1‑100 (MQ=255)
                                                    atTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CGGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCAc                                             <  1:62052/100‑1 (MQ=255)
                                                     tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGACATGTTGGTTCCTGCACCAGATTACCCAc                                             >  1:44799/1‑99 (MQ=255)
                                                     tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCAct                                            >  2:687763/1‑100 (MQ=255)
                                                       tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCActct                                          >  1:2122928/1‑100 (MQ=255)
                                                       tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCActct                                          <  1:1694957/100‑1 (MQ=255)
                                                        aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTg                                         >  2:42336/1‑100 (MQ=255)
                                                         cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                        >  2:818989/1‑100 (MQ=255)
                                                         cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                        <  1:782359/100‑1 (MQ=255)
                                                           tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGAc                                      >  2:1662918/1‑100 (MQ=255)
                                                            cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGCTTACCCACTCTGGAcc                                     >  1:1515441/1‑100 (MQ=255)
                                                               cAATGGTGTCTCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                  >  1:2557664/1‑100 (MQ=255)
                                                                    gtgtATCGGAGCTTATCGTTCAGGCAATGCAGGCAGTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGcggcgg                              <  2:33848/99‑1 (MQ=255)
                                                                        aTCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                         <  2:2213647/100‑1 (MQ=255)
                                                                           ggAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCtt                      >  1:698005/1‑100 (MQ=255)
                                                                            gtccccgctgTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     >  2:1100310/10‑100 (MQ=21)
                                                                            gtccccgctgTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     >  2:978071/10‑100 (MQ=21)
                                                                            gAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     <  2:1231143/100‑1 (MQ=255)
                                                                              gCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGAGCGGGGCGGTTTCGCTTTcc                   >  1:2139008/1‑100 (MQ=255)
                                                                              gCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                   <  1:421499/100‑1 (MQ=255)
                                                                              ccccgctgTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTcc                   >  2:2281967/8‑100 (MQ=25)
                                                                               cgcggggtacAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCttt                     <  1:1466382/88‑1 (MQ=255)
                                                                                ttATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                 <  2:469112/100‑1 (MQ=255)
                                                                                ttATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAg                 >  2:2389104/1‑100 (MQ=255)
                                                                                 tATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                >  1:1888160/1‑100 (MQ=255)
                                                                                 tATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                >  2:773334/1‑100 (MQ=255)
                                                                                  gctgTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg               >  2:735361/4‑100 (MQ=37)
                                                                                  gctgTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg               >  2:1578196/4‑100 (MQ=37)
                                                                                  gctgTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg               >  1:1971589/4‑100 (MQ=37)
                                                                                   tCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCgg              <  2:123148/100‑1 (MQ=255)
                                                                                     gTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa            >  2:1671239/1‑100 (MQ=255)
                                                                                     gTTCAGGCAATGCAGGCATTG‑CT‑CACACGGGGGACCCCAATGATGGTTACTAAAACAGATTACCCACTCTGAACGGCGGCGAGTTCGGTTTACAGCCGt             >  1:483254/1‑99 (MQ=255)
                                                                                     gTTCA‑GCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTa            >  2:2315693/1‑100 (MQ=38)
                                                                                     gTTCA‑GCAATGCCTGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTaa           >  2:896085/1‑100 (MQ=255)
                                                                                        tcagCAATGCCTGCATTGCCTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGATTTCCAGCGGTaaa          >  2:2356006/4‑100 (MQ=255)
                                                                                         aGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc        <  2:758508/100‑1 (MQ=255)
                                                                                         aGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc        >  2:343516/1‑100 (MQ=255)
                                                                                         aGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc        <  1:1026544/100‑1 (MQ=255)
                                                                                          ggCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCg       >  2:189822/1‑100 (MQ=255)
                                                                                          ggCAATGCAGGCATTC‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCg       <  2:2557683/100‑1 (MQ=255)
                                                                                            cAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt     <  1:465966/100‑1 (MQ=255)
                                                                                             aaTGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg    <  1:623423/100‑1 (MQ=255)
                                                                                              aTGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc   <  2:2214550/100‑1 (MQ=255)
                                                                                              aTGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc   <  1:1189176/100‑1 (MQ=255)
                                                                                               tGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa  <  2:350988/100‑1 (MQ=255)
                                                                                                gCAGGCAATG‑GTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTTGCTTTACCGCTGGAAAGCgaaac   <  1:825849/98‑5 (MQ=255)
                                                                                                 cAGGCACTG‑CTGAACAGCGGGGG‑CGAAATGTTGGTTCATGCACCAGATTACACAGTTTGAACCGCGGCGATTTCGCCTTCCAGCGGTAAAGCGGTGc   >  1:1718804/1‑97 (MQ=255)
                                                                                                   |                                                                                                 
CATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGA‑CGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCA  >  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: ‑