Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,199,050 G→A 5.5% L87L (CTG→CTA) ‡ 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,0500GA 5.5% 231.8 / 3.6 72L87L (CTG→CTA) ‡nfiendonuclease V; deoxyinosine 3' endonuclease
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base A (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 = 7.26e-01

GCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAACGGCTCTGCG  >  NC_000913/4198906‑4199198
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                  >  3:15465/1‑133 (MQ=255)
                                tgctgcTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTc                                                                                                                                  <  4:15465/133‑1 (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                                                                                               >  4:447509/1‑138 (MQ=255)
                                                              tCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATg                                                                                               <  3:447509/138‑1 (MQ=255)
                                                                gAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTc                                                                                                  <  3:115611/133‑1 (MQ=255)
                                                                gAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTc                                                                                                  >  4:115611/1‑133 (MQ=255)
                                                                   tataAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCg                                                                               >  3:312007/1‑149 (MQ=255)
                                                                     taAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGa                                                                                                                      <  3:34785/108‑1 (MQ=255)
                                                                     taAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGa                                                                                                                      >  4:34785/1‑108 (MQ=255)
                                                                        aGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGgaga                                                                                                                                         <  1:14542/86‑1 (MQ=255)
                                                                        aGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGgaga                                                                                                                                         >  2:14542/1‑86 (MQ=255)
                                                                              ccGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCg                                                                    >  1:380032/1‑149 (MQ=255)
                                                                              ccGCATCGCCACCACCATGCCTTACATTCCAGGTTTTATTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCg                                                                    >  4:48818/1‑149 (MQ=255)
                                                                                gCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTg                                                                       <  2:471589/144‑1 (MQ=255)
                                                                                gCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTg                                                                       >  1:471589/1‑144 (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:224396/149‑1 (MQ=255)
                                                                                               tGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGc                                                   <  2:183570/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                                                                >  2:58658/1‑128 (MQ=255)
                                                                                                       aTTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCgc                                                                <  1:58658/128‑1 (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                                                                                                >  1:300626/1‑86 (MQ=255)
                                                                                                                 ttctttcCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCAt                                                                                                <  2:300626/86‑1 (MQ=255)
                                                                                                                      tccttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCAt                            >  2:212172/1‑149 (MQ=255)
                                                                                                                        cttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATtg                          >  4:3981/1‑149 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCg                                                                            >  2:243548/1‑98 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCg                                                                            <  1:243548/98‑1 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTAtt                                              >  3:157959/1‑128 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTAtt                                              <  4:157959/128‑1 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGAcc                              <  3:107665/144‑1 (MQ=255)
                                                                                                                         ttccGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGAcc                              >  4:107665/1‑144 (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)
                                                                                                                                                ggCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAACGGCTCTGCg  <  4:50924/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGCTGGCAGCGTGGGAGATGCTGTCGCAAAAGCCGGATTTAGTGTTTGTCGATGGTCATGGGATCTCGCATCCTCGCCGTCTTGGCGTCGCCAGCCATTTTGGCTTATTGGTGGATGTGCCGACCATTGGCGTGGCGAAAAAACGGCTCTGCG  >  NC_000913/4198906‑4199198

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: