Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,899 G→T 6.2% 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,8990GT 6.2% 218.6 / 4.8 65Q113H (CAG→CAT) ‡alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base T (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

AAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGA  >  NC_000913/2407757‑2408036
                                                                                                                                              |                                                                                                                                          
aaGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTg                                                                                                                                      <  1:170613/149‑1 (MQ=255)
  gggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCCTGTATTG‑cc                                                                                                                                   <  2:226056/149‑2 (MQ=255)
  gggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑Ct                                                                                                                                   <  3:637118/149‑1 (MQ=255)
    gTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGa                                                                                                                                 >  1:556075/1‑149 (MQ=255)
         gCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGc                                                                                                                            <  3:350103/149‑1 (MQ=255)
           gATTCCAGAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGATCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCGATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                          >  4:585664/1‑149 (MQ=255)
           gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCGGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                          <  3:610360/149‑1 (MQ=255)
           gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                          <  2:598596/149‑1 (MQ=255)
            aTTCCAAAGGTCTTTACTCCGCGCGTAAAACCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCggg                                                                                                                         <  1:410164/149‑1 (MQ=255)
               ccAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGAc                                                                                                                      <  2:107441/149‑1 (MQ=255)
                 aaaGGTCTATACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGa                                                                                                                    >  3:486313/1‑149 (MQ=255)
                   aGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGaaa                                                                                                                  <  2:556075/149‑1 (MQ=255)
                      tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGt                                                                                                               <  4:522922/149‑1 (MQ=255)
                        tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                             <  1:522098/149‑1 (MQ=255)
                                    gTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGaaa                                                                                                                  <  1:16219/132‑1 (MQ=255)
                                    gTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGaaa                                                                                                                  >  2:16219/1‑132 (MQ=255)
                                       aaGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGAt                                                                                              <  3:585664/149‑1 (MQ=255)
                                           catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTgg                                                                                                            >  1:125059/1‑131 (MQ=255)
                                           catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTgg                                                                                                            <  2:125059/131‑1 (MQ=255)
                                           catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAcc                                                                                          >  4:231260/1‑149 (MQ=255)
                                           catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAcc                                                                                          >  1:563697/1‑149 (MQ=255)
                                           catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgata                                                                                                                           <  4:227315/116‑5 (MQ=38)
                                           catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgata                                                                                                                           >  3:227315/1‑112 (MQ=38)
                                              catGCAGCACTACCAGGCTCGTGGCATGCGTGATTTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAc                                                                                       <  4:486313/149‑1 (MQ=255)
                                                 gcagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                                    >  4:291505/1‑149 (MQ=255)
                                                      aCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGAcc                                                                               >  3:176943/1‑149 (MQ=255)
                                                            aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                                         <  3:231260/149‑1 (MQ=255)
                                                            aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                                         >  4:673067/1‑149 (MQ=255)
                                                            aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                                         >  4:90889/1‑149 (MQ=255)
                                                             ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcgg                                                                           <  2:382081/146‑1 (MQ=255)
                                                             ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcgg                                                                           >  1:382081/1‑146 (MQ=255)
                                                              gCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGt                                                                       <  2:450632/149‑1 (MQ=255)
                                                               cTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGtt                                                                      <  2:52465/149‑1 (MQ=255)
                                                                 cGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTc                                                                    <  2:367045/149‑1 (MQ=255)
                                                                            tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                         <  3:152879/149‑1 (MQ=255)
                                                                                 ttACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTaaa                                                    <  1:182044/149‑1 (MQ=255)
                                                                                   aCCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc                                                  >  4:476607/1‑149 (MQ=255)
                                                                                   aCCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc                                                  <  3:476607/149‑1 (MQ=255)
                                                                                      gTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt                                               >  1:101389/1‑149 (MQ=255)
                                                                                         gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggc                                                                          <  4:488336/119‑1 (MQ=255)
                                                                                         gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggc                                                                          >  3:488336/1‑119 (MQ=255)
                                                                                          aaGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCAt                                           <  4:3898/149‑1 (MQ=255)
                                                                                            gATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTa                                         <  1:646030/149‑1 (MQ=255)
                                                                                            gATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTa                                         >  3:55438/1‑149 (MQ=255)
                                                                                              tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATc                                       <  3:291505/149‑1 (MQ=255)
                                                                                                     aTCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                <  3:629893/149‑1 (MQ=255)
                                                                                                       cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATg                              >  1:208893/1‑149 (MQ=255)
                                                                                                         gCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGaa                            >  4:497917/1‑149 (MQ=255)
                                                                                                                   aTCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACt                  <  4:176943/149‑1 (MQ=255)
                                                                                                                    tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCctc                       >  3:568136/1‑143 (MQ=255)
                                                                                                                    tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCctc                       <  4:568136/143‑1 (MQ=255)
                                                                                                                    tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTg                 >  2:258964/1‑149 (MQ=255)
                                                                                                                    tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTg                 >  2:686264/1‑149 (MQ=255)
                                                                                                                         gCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTcc            >  2:262102/1‑149 (MQ=255)
                                                                                                                           ttATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCg          <  4:55438/149‑1 (MQ=255)
                                                                                                                             aTCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTGGCGATGAATCCTCTGACTGGTTCCCAGa        >  2:373761/1‑149 (MQ=255)
                                                                                                                              tCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGAc       <  2:101389/149‑1 (MQ=255)
                                                                                                                                gTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCt     >  2:535226/1‑149 (MQ=255)
                                                                                                                                gTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCt     >  3:484807/1‑149 (MQ=255)
                                                                                                                                   cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga  >  1:331529/1‑149 (MQ=255)
                                                                                                                                    cagCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga  >  2:439502/3‑149 (MQ=255)
                                                                                                                                    aGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTc                                                              >  1:673394/1‑88 (MQ=255)
                                                                                                                                    aGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTc                                                              <  2:673394/88‑1 (MQ=255)
                                                                                                                                     ggCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATc                                       >  3:78647/1‑110 (MQ=255)
                                                                                                                                     ggCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATc                                       <  4:78647/110‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                          
AAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGA  >  NC_000913/2407757‑2408036

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