Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,199,049 T→C 5.5% L87P (CTG→CCG) ‡ nfi → endonuclease V; deoxyinosine 3' endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,199,0490TC 5.5% 238.3 / 3.4 72L87P (CTG→CCG) ‡nfiendonuclease V; deoxyinosine 3' endonuclease
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (2/2);  total (35/37)
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.76e-01

GCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAACG  >  NC_000913/4198906‑4199190
                                                                                                                                               |                                                                                                                                             
gCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAg                                                                                                                                          >  1:265301/1‑149 (MQ=255)
   gggCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGt                                                                                                                                       <  2:265301/149‑1 (MQ=255)
    ggCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCCAGCAGCgca                                                                                                                                      <  2:241350/149‑3 (MQ=255)
     gCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTgg                                                                                                                                     <  3:252319/149‑1 (MQ=255)
       ggAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGga                                                                                                                                   <  2:278028/149‑1 (MQ=255)
          gaAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGAt                                                                                                                                >  3:203671/1‑149 (MQ=255)
                aCGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTc                                                                                                                          <  2:51587/149‑1 (MQ=255)
                aCGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGCTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGGGGGAGATGCTGTc                                                                                                                          >  3:46880/1‑149 (MQ=255)
                      gcggcgATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTggg                                                                                                                                    <  2:288116/133‑1 (MQ=255)
                      gcggcgATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTggg                                                                                                                                    >  1:288116/1‑133 (MQ=255)
                      gcggcgATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCaaaa                                                                                                                    <  4:140602/149‑1 (MQ=255)
                                tgctgcTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTc                                                                                                                          <  4:15465/133‑1 (MQ=255)
                                tgctgcTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTc                                                                                                                          >  3:15465/1‑133 (MQ=255)
                                      tgATATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTACCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGc                                                                                                                  <  1:25010/135‑1 (MQ=255)
                                      tgAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTGCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGc                                                                                                                  >  2:25010/1‑135 (MQ=255)
                                      tgAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGttt                                                                                                    >  2:453639/1‑149 (MQ=255)
                                        aaaTATCCCTCGCTTGAGCCGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGt                                                                                                  <  1:41290/149‑1 (MQ=255)
                                             tCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATg                                                                                             <  4:115574/149‑1 (MQ=255)
                                               ccTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGa                                                                                               <  2:111082/145‑1 (MQ=255)
                                               ccTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGa                                                                                               >  1:111082/1‑145 (MQ=255)
                                                cTCGCTTGAGCTGGTTGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTc                                                                                          >  3:204440/1‑149 (MQ=255)
                                                cTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCGGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTc                                                                                          <  4:217755/149‑1 (MQ=255)
                                                 tCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCa                                                                                         <  1:43152/149‑1 (MQ=255)
                                                       gAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGAt                                                                                   <  4:203671/149‑1 (MQ=255)
                                                         gCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAgt                                                                                                        >  4:230065/1‑126 (MQ=255)
                                                         gCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAgt                                                                                                        <  3:230065/126‑1 (MQ=255)
                                                            ggTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGc                                                                              <  2:119841/149‑1 (MQ=255)
                                                              tCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATg                                                                                       <  3:447509/138‑1 (MQ=255)
                                                              tCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATg                                                                                       >  4:447509/1‑138 (MQ=255)
                                                                gAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTc                                                                                          >  4:115611/1‑133 (MQ=255)
                                                                gAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTc                                                                                          <  3:115611/133‑1 (MQ=255)
                                                                   tataAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCg                                                                       >  3:312007/1‑149 (MQ=255)
                                                                     taAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGa                                                                                                              <  3:34785/108‑1 (MQ=255)
                                                                     taAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGa                                                                                                              >  4:34785/1‑108 (MQ=255)
                                                                        aGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGgaga                                                                                                                                 >  2:14542/1‑86 (MQ=255)
                                                                        aGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGgaga                                                                                                                                 <  1:14542/86‑1 (MQ=255)
                                                                              ccGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCg                                                            >  1:380032/1‑149 (MQ=255)
                                                                              ccGCATCGCCACCACCATGCCTTACATTCCAGGTTTTATTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCg                                                            >  4:48818/1‑149 (MQ=255)
                                                                                gCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTg                                                               >  1:471589/1‑144 (MQ=255)
                                                                                gCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTg                                                               <  2:471589/144‑1 (MQ=255)
                                                                                      ccaccaccaTGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGAt                                                                                   <  1:232208/118‑1 (MQ=255)
                                                                                      ccaccaccaTGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGAt                                                                                   >  2:232208/1‑118 (MQ=255)
                                                                                            ccaTGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCAtttt                                              <  1:46496/149‑1 (MQ=255)
                                                                                             caTGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGCTGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTg                                             >  3:436283/1‑149 (MQ=255)
                                                                                               tGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGc                                           <  2:183570/149‑1 (MQ=255)
                                                                                               tGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGc                                           <  2:224396/149‑1 (MQ=255)
                                                                                                  cTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTa                                        >  1:341891/1‑149 (MQ=255)
                                                                                                   ttACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTAt                                       >  2:441986/1‑149 (MQ=255)
                                                                                                    cccATCCCAGGTTTTCTTTCCTTGCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGGTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCCTTTTGTCTTAtt                                      <  4:46880/147‑1 (MQ=255)
                                                                                                       aTTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCgc                                                        <  1:58658/128‑1 (MQ=255)
                                                                                                       aTTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCgc                                                        >  2:58658/1‑128 (MQ=255)
                                                                                                         tCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATtggtgg                                 >  4:37942/1‑149 (MQ=255)
                                                                                                          ccAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGa                                >  3:50924/1‑149 (MQ=255)
                                                                                                          ccAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGa                                <  4:312007/149‑1 (MQ=255)
                                                                                                           cAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGAt                               <  1:441986/149‑1 (MQ=255)
                                                                                                            aGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGAtg                              >  2:393935/1‑149 (MQ=255)
                                                                                                                 ttctttcCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCAt                                                                                        <  2:300626/86‑1 (MQ=255)
                                                                                                                 ttctttcCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCAt                                                                                        >  1:300626/1‑86 (MQ=255)
                                                                                                                      tccttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCAt                    >  2:212172/1‑149 (MQ=255)
                                                                                                                        cttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATtg                  >  4:3981/1‑149 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCg                                                                    <  1:243548/98‑1 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCg                                                                    >  2:243548/1‑98 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTAtt                                      >  3:157959/1‑128 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTAtt                                      <  4:157959/128‑1 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGAcc                      >  4:107665/1‑144 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGAcc                      <  3:107665/144‑1 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATtgg                 >  3:112345/1‑149 (MQ=255)
                                                                                                                                   atCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGaaa       <  4:204440/149‑1 (MQ=255)
                                                                                                                                       tGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAAc   <  4:436283/149‑1 (MQ=255)
                                                                                                                                       cacgcTGCCAGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATtg                  >  4:386731/3‑134 (MQ=255)
                                                                                                                                       cacgcTGCCAGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATtg                  <  3:386731/132‑1 (MQ=255)
                                                                                                                                        acgcTGCCAGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAACg  >  3:208263/2‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAACG  >  NC_000913/4198906‑4199190

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