Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,898 A→C 6.3% 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,8980AC 6.3% 484.7 / 23.0 144Q113P (CAG→CCG) ‡alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (69/66);  new base C (5/4);  total (74/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 148 total aligned reads displayed.

ACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGC  >  NC_000913/2407750‑2408046
                                                                                                                                                    |                                                                                                                                                     
  cGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                                                     <  3:925397/149‑1 (MQ=255)
     tCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAt                                                                                                                                                  >  2:1395154/1‑149 (MQ=255)
       aaGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTg                                                                                                                                                >  3:255682/1‑149 (MQ=255)
       aaGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTg                                                                                                                                                >  4:123492/1‑149 (MQ=255)
          ggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTg                                                                                                                                            >  1:1050082/1‑149 (MQ=255)
            tATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTG‑cctgaa                                                                                                                                         <  2:974154/149‑6 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                                    >  2:1039481/1‑149 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                                    >  3:302822/1‑149 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                                    >  1:259964/1‑149 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                                    <  2:327392/149‑1 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                                    <  2:298660/149‑1 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgataa                                                                                                                                    <  1:230159/149‑6 (MQ=255)
                   aTTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCggg                                                                                                                                   >  3:544955/1‑149 (MQ=255)
                        aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCTTGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGa                                                                                                                              >  4:718897/1‑149 (MQ=255)
                          aGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGaaa                                                                                                                            <  2:130153/149‑1 (MQ=255)
                              cTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGtt                                                                                                                        <  1:138241/149‑1 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                       <  4:301527/149‑1 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                       <  3:428407/149‑1 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                       >  1:702495/1‑149 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                       <  3:1377579/149‑1 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                       >  2:349652/1‑149 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                       >  2:730538/1‑149 (MQ=255)
                                ttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTgg                                                                                                                      <  3:951346/149‑1 (MQ=255)
                                 tACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGt                                                                                                                     <  1:1051546/149‑1 (MQ=255)
                                  aCTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGtt                                                                                                                    >  1:935509/1‑149 (MQ=255)
                                    tCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGt                                                                                                                         >  3:1365938/1‑142 (MQ=255)
                                    tCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGt                                                                                                                         <  4:1365938/142‑1 (MQ=255)
                                    tCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                                                  >  1:978132/1‑149 (MQ=255)
                                     ccGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCt                                                                                                                 <  4:325076/149‑1 (MQ=255)
                                                  catcatGCAGCACTACCAGGCTCGTGGTATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAg                                                                                                                                                      <  2:399913/100‑1 (MQ=255)
                                                  catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTTAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAcc                                                                                                    >  3:410903/1‑149 (MQ=255)
                                                     catGCAGCACTACCAGGCCCGTGGCATGCTTGATGTTACCGTGAAAGATATTTACTTCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCTGATTACCCAc                                                                                                 <  2:834623/149‑1 (MQ=255)
                                                      atGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAct                                                                                                >  2:169993/1‑149 (MQ=255)
                                                        gcagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                                              >  4:98302/1‑149 (MQ=255)
                                                        gcagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                                              >  1:359554/1‑149 (MQ=255)
                                                           gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGa                                                                                           <  3:351299/149‑1 (MQ=255)
                                                           gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGa                                                                                           <  1:709862/149‑1 (MQ=255)
                                                           gcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGa                                                                                           <  4:323343/149‑1 (MQ=255)
                                                               tACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgc                                                                                       >  1:1476216/1‑149 (MQ=255)
                                                                aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                                                                      >  1:748165/1‑149 (MQ=255)
                                                                   aGGCTCGTGGCATTCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                                                   >  3:1101636/1‑149 (MQ=255)
                                                                       tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                                                  <  3:33966/114‑1 (MQ=255)
                                                                       tCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                                                  >  4:33966/1‑114 (MQ=255)
                                                                                gCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCa                                                                      <  1:169993/149‑1 (MQ=255)
                                                                                   tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCAGCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                   >  1:739839/1‑149 (MQ=255)
                                                                                    gATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCgg                                                                  >  3:521329/1‑149 (MQ=255)
                                                                                     aTGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGt                                                                 <  4:521329/149‑1 (MQ=255)
                                                                                         tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAg                                                             <  3:467622/149‑1 (MQ=255)
                                                                                         tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAg                                                             <  4:446664/149‑1 (MQ=255)
                                                                                          aCCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc                                                            >  1:1456272/1‑149 (MQ=255)
                                                                                             gTGGAAGATATTTACGTCGGCAATGGTGTAGCGGAGCTTATCGTTCAGGCAACGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCGCTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt                                                         <  4:437408/149‑1 (MQ=255)
                                                                                             gTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt                                                         <  2:359554/149‑1 (MQ=255)
                                                                                              tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActc                                                                                               >  3:1385423/1‑110 (MQ=255)
                                                                                              tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActc                                                                                               <  4:1385423/110‑1 (MQ=255)
                                                                                              tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg                                                        >  3:1253090/1‑149 (MQ=255)
                                                                                              tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg                                                        >  3:155538/1‑149 (MQ=255)
                                                                                              tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATACAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg                                                        >  3:671801/1‑149 (MQ=255)
                                                                                                gacgATATTTACATCGGCAATGGAGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa                                                      >  2:151604/4‑149 (MQ=255)
                                                                                                gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCggg                                                                                                                                   >  1:548263/1‑72 (MQ=255)
                                                                                                gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa                                                      >  4:90089/1‑149 (MQ=255)
                                                                                                 aaGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGa                                                                                                         >  2:668591/1‑97 (MQ=255)
                                                                                                 aaGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGa                                                                                                         <  1:668591/97‑1 (MQ=255)
                                                                                                    atatTTACATCGGCAATGGGGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACGAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTAt                                                  <  2:105541/149‑1 (MQ=255)
                                                                                                     tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATc                                                 >  2:775350/1‑149 (MQ=255)
                                                                                                     tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATc                                                 <  1:609306/149‑1 (MQ=255)
                                                                                                        ttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCttt                                              <  2:1486233/149‑1 (MQ=255)
                                                                                                         tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                                             <  2:54388/149‑1 (MQ=255)
                                                                                                          aCATCGGCAATGGTGTATCGGGGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAACGGTAAAGCGGTGCATTATCTTTGc                                            <  4:1384278/149‑1 (MQ=255)
                                                                                                          aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAACCGGGGCTTTATTTTTGc                                            >  4:1280179/1‑149 (MQ=255)
                                                                                                            aTCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                          >  1:115909/1‑149 (MQ=255)
                                                                                                             tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGAt                                         >  2:920006/1‑149 (MQ=255)
                                                                                                                gCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACCCTGGACCGCGGCGGATTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGag                                      >  3:131505/1‑148 (MQ=255)
                                                                                                                  aaTGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATc                                    <  4:18807/149‑1 (MQ=255)
                                                                                                                   aTGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATcc                                   <  2:56729/149‑1 (MQ=255)
                                                                                                                     ggTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCctc                                 <  4:27285/149‑1 (MQ=255)
                                                                                                                     ggTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCctc                                 <  4:813718/149‑1 (MQ=255)
                                                                                                                      gtgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCctct                                <  2:1378832/149‑1 (MQ=255)
                                                                                                                       tgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTg                               >  1:832775/1‑149 (MQ=255)
                                                                                                                       tgtAACGGAGCTTATCGTTCAGGCAAGGCAGGCCTTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGCCCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGCATCCTCTg                               >  2:1099221/1‑149 (MQ=255)
                                                                                                                         tATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGAc                             <  3:449739/149‑1 (MQ=255)
                                                                                                                             ggAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGCTGAATCCTCTGACTGGt                         >  2:312726/1‑149 (MQ=255)
                                                                                                                              gAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt                                                         >  3:1200892/1‑116 (MQ=255)
                                                                                                                                 cTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGc                                                                    <  4:1316783/102‑1 (MQ=255)
                                                                                                                                    aTCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGa                  <  2:859679/149‑1 (MQ=255)
                                                                                                                                      tgTTCA‑GCAATGCCTGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCt               >  2:1058701/2‑149 (MQ=255)
                                                                                                                                       gTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCt               >  4:407446/1‑149 (MQ=255)
                                                                                                                                        ttCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGAt                                         >  3:988244/1‑122 (MQ=255)
                                                                                                                                        ttCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTc              >  2:733008/1‑149 (MQ=255)
                                                                                                                                         ttcagCAATGCCTGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCg             >  3:495594/5‑149 (MQ=255)
                                                                                                                                          cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgt            <  4:255682/149‑2 (MQ=255)
                                                                                                                                          cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga            <  1:1118387/149‑1 (MQ=255)
                                                                                                                                          cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga            <  2:1230733/149‑1 (MQ=255)
                                                                                                                                          cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga            <  2:936080/149‑1 (MQ=255)
                                                                                                                                           cagCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga            >  3:563375/3‑149 (MQ=255)
                                                                                                                                           aGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgat           <  3:718897/149‑1 (MQ=255)
                                                                                                                                           aGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgat           <  1:37471/149‑1 (MQ=255)
                                                                                                                                             gCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgatga         >  1:504466/1‑149 (MQ=255)
                                                                                                                                               aaTGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGata       >  1:1445531/1‑149 (MQ=255)
                                                                                                                                                  gCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTc    <  4:469807/149‑1 (MQ=255)
                                                                                                                                                   cAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGTACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTcg   >  2:748228/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGC  >  NC_000913/2407750‑2408046

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