Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,898 A→C 8.0% 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 8.0% 204.8 / 8.2 62Q113P (CAG→CCG) ‡alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base C (3/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCG  >  NC_000913/2407750‑2408045
                                                                                                                                                    |                                                                                                                                                    
aCCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCa                                                                                                                                                      <  1:497950/149‑1 (MQ=255)
         gggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑Ct                                                                                                                                            >  3:457680/1‑149 (MQ=255)
          ggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTg                                                                                                                                           >  2:283987/1‑149 (MQ=255)
             aTTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                                                    <  4:71799/138‑1 (MQ=255)
             aTTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                                                    >  3:71799/1‑138 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTACAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCCTGCATTG‑cctgaa                                                                                                                                        <  3:68468/144‑6 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCCTGCATTG‑cctgaa                                                                                                                                        >  4:68468/1‑139 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgg                                                                                                                                   <  3:50338/149‑1 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgataa                                                                                                                                   <  2:469373/149‑6 (MQ=255)
                  gATTCCAAAGGTCTTTACTCCGCGCGGAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGGTGGTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCTTTG‑CTGAACAGCgg                                                                                                                                   <  3:51738/149‑1 (MQ=255)
                    ttCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCgggg                                                                                                                                 >  3:114091/1‑149 (MQ=255)
                     tCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGa                                                                                                                                <  2:374935/149‑1 (MQ=255)
                        aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑c                                                                                                                                             <  4:322742/133‑1 (MQ=255)
                        aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑c                                                                                                                                             >  3:322742/1‑133 (MQ=255)
                        aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGa                                                                                                                             >  4:388953/1‑149 (MQ=255)
                         aaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGaa                                                                                                                            >  2:303074/1‑149 (MQ=255)
                          aGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGaaa                                                                                                                           <  2:363821/149‑1 (MQ=255)
                            gTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATg                                                                                                                         >  4:655577/1‑149 (MQ=255)
                              cTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGtt                                                                                                                       <  4:325355/149‑1 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                      <  4:592420/149‑1 (MQ=255)
                               tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                      >  4:351039/1‑149 (MQ=255)
                                ttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTgg                                                                                                                     <  3:388753/149‑1 (MQ=255)
                                    tCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                                                 <  2:635152/149‑1 (MQ=255)
                                             aaaGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGa                                                                                                        <  1:427492/149‑1 (MQ=255)
                                             aaaGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGa                                                                                                        <  1:303074/149‑1 (MQ=255)
                                                   atcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAccc                                                                                                  <  3:655577/149‑1 (MQ=255)
                                                       tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGTTTACCCActc                                                                                              >  3:432660/1‑149 (MQ=255)
                                                          agcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTgg                                                                                           <  3:23572/149‑1 (MQ=255)
                                                            caCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGAc                                                                                         <  3:3984/149‑1 (MQ=255)
                                                               tACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgc                                                                                      >  2:653900/1‑149 (MQ=255)
                                                                          tGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGc                                                                           <  2:142837/149‑1 (MQ=255)
                                                                              aTGCGTGATGTTACCGTGGAAGATATTTACATCGGCAAAGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTc                                                                       >  1:539682/1‑149 (MQ=255)
                                                                                   tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                  <  2:232679/149‑1 (MQ=255)
                                                                                         tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAg                                                            >  2:356722/1‑149 (MQ=255)
                                                                                            cGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                  <  4:249090/140‑1 (MQ=255)
                                                                                            cGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                  >  3:249090/1‑140 (MQ=255)
                                                                                             gTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt                                                        <  2:111660/149‑1 (MQ=255)
                                                                                               ggAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGc                                                      <  3:601638/149‑1 (MQ=255)
                                                                                                gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGt                                                                <  1:4874/138‑1 (MQ=255)
                                                                                                gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGt                                                                >  2:4874/1‑138 (MQ=255)
                                                                                                gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa                                                     >  1:388955/1‑149 (MQ=255)
                                                                                                   gATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTa                                                  >  1:522704/1‑149 (MQ=255)
                                                                                                       tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCtt                                              >  4:521955/1‑149 (MQ=255)
                                                                                                         tACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                                            <  3:38600/149‑1 (MQ=255)
                                                                                                          aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                           <  4:671205/149‑1 (MQ=255)
                                                                                                          aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                           <  4:432660/149‑1 (MQ=255)
                                                                                                          aCAGCGGCAATGGTGTATCGGAGCATATCGTACAGGCAATGCAGGCGTTG‑CTGAACAGCGTGGACGAAATGTTGGTTCCTGCATCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTACAGCGGTGCATTATCGTTGc                                           <  2:571310/149‑1 (MQ=255)
                                                                                                            aTCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                         >  1:277695/1‑149 (MQ=255)
                                                                                                            aTCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                         <  2:106814/149‑1 (MQ=255)
                                                                                                            aTCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                         >  1:683308/1‑149 (MQ=255)
                                                                                                             tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGAt                                        >  3:526532/1‑149 (MQ=255)
                                                                                                                  aaTGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATc                                   <  1:283987/149‑1 (MQ=255)
                                                                                                                           tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTg                          <  1:546073/149‑1 (MQ=255)
                                                                                                                                 cTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTccc                    >  4:433618/1‑149 (MQ=255)
                                                                                                                                  ttATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCg                   >  3:251142/1‑149 (MQ=255)
                                                                                                                                    aTCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGa                 >  4:547275/1‑149 (MQ=255)
                                                                                                                                         ttcagCAATGCCTGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCg            >  4:115882/5‑149 (MQ=255)
                                                                                                                                          cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga           <  2:521847/149‑1 (MQ=255)
                                                                                                                                           aGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgat          <  3:521955/149‑1 (MQ=255)
                                                                                                                                            agCAATGCCTGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgatg         >  2:63030/2‑149 (MQ=255)
                                                                                                                                                aTGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGatat     >  2:693538/1‑149 (MQ=255)
                                                                                                                                                   cAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTcg  >  4:610941/1‑149 (MQ=255)
                                                                                                                                                   cAGGCATTG‑CTCAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTTGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGTTTCCCGGACCTCGATGATATTcg  <  2:539682/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCG  >  NC_000913/2407750‑2408045

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: ‑