Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,007,054 C→T 6.4% C397C (TGC→TGT)  recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,0540CT 6.4% 310.9 / 7.6 94C397C (TGC→TGT) recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (48/40);  new base T (3/3);  total (51/43)
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.65e-01

CGCCGTTGTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAC  >  NC_000913/4006908‑4007195
                                                                                                                                                  |                                                                                                                                             
cgcCGTTGTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCgg                                                                                                                                             <  1:825853/149‑1 (MQ=255)
 gcCGTTGTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCggg                                                                                                                                            >  1:252668/1‑149 (MQ=255)
 gcCGTTGTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCggg                                                                                                                                            >  4:519136/1‑149 (MQ=255)
   cGTTGTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGaa                                                                                                                                          >  4:289775/1‑149 (MQ=255)
       gTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGc                                                                                                                                      <  1:331265/149‑1 (MQ=255)
         cTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGa                                                                                                                                    <  3:838184/149‑1 (MQ=255)
                agaAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGc                                                                                                                             >  3:23706/1‑149 (MQ=255)
                 gaAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGcc                                                                                                                            >  1:692080/1‑149 (MQ=255)
                  aaGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCt                                                                                                                           >  2:267338/1‑149 (MQ=255)
                  aaGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCt                                                                                                                           >  2:708097/1‑149 (MQ=255)
                  aaGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCt                                                                                                                           >  3:151335/1‑149 (MQ=255)
                  aaGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCt                                                                                                                           >  3:388571/1‑149 (MQ=255)
                    gCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCg                                                                                                                         >  4:427504/1‑149 (MQ=255)
                     ccGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGa                                                                                                                        >  4:428418/1‑149 (MQ=255)
                          ggggCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATccgc                                                                                                                   >  3:79698/1‑149 (MQ=255)
                           gggCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATccgcc                                                                                                                  >  2:399043/1‑149 (MQ=255)
                             gCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGa                                                                                                                >  1:296019/1‑149 (MQ=255)
                             gCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGa                                                                                                                >  4:611642/1‑149 (MQ=255)
                                gTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAAc                                                                                                             >  2:829880/1‑149 (MQ=255)
                                  tGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAg                                                                                                           <  3:90889/149‑1 (MQ=255)
                                   gCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCAGAAACAGt                                                                                                          <  4:23706/149‑1 (MQ=255)
                                      ggATATCGAGCGCCACGAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTacg                                                                                                       >  2:400924/1‑149 (MQ=255)
                                      ggATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTacg                                                                                                       >  3:354734/1‑149 (MQ=255)
                                      ggATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTacg                                                                                                       >  4:832279/1‑149 (MQ=255)
                                        atatCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGTGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTacgac                                                                                                     >  3:617932/1‑149 (MQ=255)
                                        atatCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGTGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTacgac                                                                                                     >  1:142709/1‑149 (MQ=255)
                                          atCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACgg                                                                                                   >  2:725412/1‑149 (MQ=255)
                                            cGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGTGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGtt                                                                                                 <  2:142709/149‑1 (MQ=255)
                                            cGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGTGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGtt                                                                                                 <  4:617932/149‑1 (MQ=255)
                                                 gcCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAAcc                                                                                            >  2:640332/1‑149 (MQ=255)
                                                   cacaAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGa                                                                                          >  3:582177/1‑149 (MQ=255)
                                                    acaAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGAt                                                                                         >  4:214679/1‑149 (MQ=255)
                                                    acaAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGAt                                                                                         <  1:484418/149‑1 (MQ=255)
                                                      aaaCTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGc                                                                                       <  1:737447/149‑1 (MQ=255)
                                                      aaaCTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGc                                                                                       >  3:286124/1‑149 (MQ=255)
                                                        aCTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTc                                                                                     <  3:140327/149‑1 (MQ=255)
                                                          tCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAg                                                                                   <  3:611642/149‑1 (MQ=255)
                                                            aaTGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGAt                                                                                 >  2:117049/1‑149 (MQ=255)
                                                            aaTGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGAt                                                                                 >  1:17199/1‑149 (MQ=255)
                                                             aTGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGAtt                                                                                >  4:428450/1‑149 (MQ=255)
                                                              tGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTg                                                                               <  3:519136/149‑1 (MQ=255)
                                                              tGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTg                                                                               >  4:448916/1‑149 (MQ=255)
                                                               gCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGc                                                                              <  4:434287/149‑1 (MQ=255)
                                                                cGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGcc                                                                             <  4:286124/149‑1 (MQ=255)
                                                                     ggCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTc                                                                        <  3:181884/149‑1 (MQ=255)
                                                                      gcgcGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTcc                                                                       >  4:302804/1‑149 (MQ=255)
                                                                           tttGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCAt                                                                  >  3:819641/1‑149 (MQ=255)
                                                                             tGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTg                                                                <  4:582177/149‑1 (MQ=255)
                                                                                 gAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCg                                                            >  3:389436/1‑149 (MQ=255)
                                                                                 gAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCg                                                            >  4:34581/1‑149 (MQ=255)
                                                                                 gAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCg                                                            >  1:577159/1‑149 (MQ=255)
                                                                                    gcgcAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGt                                                                                                          <  3:672597/100‑1 (MQ=255)
                                                                                    gcgcAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGt                                                                                                          >  4:672597/1‑100 (MQ=255)
                                                                                       cAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGaa                                                      <  2:296019/149‑1 (MQ=255)
                                                                                       cAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGaa                                                      >  4:710684/1‑149 (MQ=255)
                                                                                        aaaCTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAAt                                                     <  3:497463/149‑1 (MQ=255)
                                                                                                cgtcgtCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATAGGTCGTGTGAATCAGCGGtt                                             >  1:388932/1‑149 (MQ=255)
                                                                                                 gtcgtcTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGttt                                            <  4:256273/149‑1 (MQ=255)
                                                                                                     tcTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGGGCCGTGCg                                                                                                                                              <  1:598955/47‑1 (MQ=255)
                                                                                                     tcTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGGGCCGTGCg                                                                                                                                              >  2:598955/1‑47 (MQ=255)
                                                                                                     tcTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGa                                        <  3:508711/149‑1 (MQ=255)
                                                                                                      cTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGAt                                       <  1:400924/149‑1 (MQ=255)
                                                                                                       tGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGtgtg                                                        >  1:767415/1‑131 (MQ=255)
                                                                                                       tGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGtgtg                                                        <  2:767415/131‑1 (MQ=255)
                                                                                                        ggTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATgg                                     <  2:692080/149‑1 (MQ=255)
                                                                                                          tATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGt                                   <  3:44543/149‑1 (MQ=255)
                                                                                                               cTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTAtgt                              <  4:343071/149‑1 (MQ=255)
                                                                                                                tGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTAtgtg                             >  1:88761/1‑149 (MQ=255)
                                                                                                                tGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTAtgtg                             >  3:359330/1‑149 (MQ=255)
                                                                                                                tGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTAtgtg                             >  1:838557/1‑149 (MQ=255)
                                                                                                                tGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGCTAtgtg                             >  2:74334/1‑149 (MQ=255)
                                                                                                                 gAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATgtgg                            <  1:812208/149‑1 (MQ=255)
                                                                                                                  aaCTATTTTGGCGAAGGGCGTCAGGAGCCGTGTGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATgg                                     <  1:796407/139‑1 (MQ=255)
                                                                                                                  aaCTATTTTGGCGAAGGGCGTCAGGAGCCGTGTGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATgg                                     >  2:796407/1‑139 (MQ=255)
                                                                                                                  aaCTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATgtggt                           >  1:106478/1‑149 (MQ=255)
                                                                                                                  aaCTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATgtggt                           <  1:74334/149‑1 (MQ=255)
                                                                                                                   aCTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATgtggtg                          <  1:760877/149‑1 (MQ=255)
                                                                                                                    cTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGtggtgg                         <  1:871839/149‑1 (MQ=255)
                                                                                                                     tATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGa                        >  2:487799/1‑149 (MQ=255)
                                                                                                                      aTTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGaa                       <  3:427504/149‑1 (MQ=255)
                                                                                                                      aTTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGaa                       <  3:289775/149‑1 (MQ=255)
                                                                                                                       ttttGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAg                      <  4:151335/149‑1 (MQ=255)
                                                                                                                       ttttGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAg                      <  4:388571/149‑1 (MQ=255)
                                                                                                                          tGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGa                   <  3:302804/149‑1 (MQ=255)
                                                                                                                             cGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTCCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTc                >  4:60408/1‑149 (MQ=255)
                                                                                                                              gAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCg               >  3:838849/1‑149 (MQ=255)
                                                                                                                                aGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCgtg             >  3:459461/1‑149 (MQ=255)
                                                                                                                                  ggCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCgtggt           <  3:448916/149‑1 (MQ=255)
                                                                                                                                    cGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGc         <  2:620033/149‑1 (MQ=255)
                                                                                                                                    cGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGc         <  2:252668/149‑1 (MQ=255)
                                                                                                                                      tCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCta       >  3:807905/1‑149 (MQ=255)
                                                                                                                                        aGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCtaat     <  4:293600/149‑1 (MQ=255)
                                                                                                                                          gAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCtaataa   <  1:639009/149‑1 (MQ=255)
                                                                                                                                           aGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAc  <  4:325599/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGCCGTTGTCTGGAAGAGAAGCCGCAGGGGCAGTTGCAGGATATCGAGCGCCACAAACTCAATGCGATGGGCGCGTTTGCCGAAGCGCAAACTTGCCGTCGTCTGGTATTGCTGAACTATTTTGGCGAAGGGCGTCAGGAGCCGTGCGGGAACTGCGATATCTGCCTCGATCCGCCGAAACAGTACGACGGTTCAACCGATGCTCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAAC  >  NC_000913/4006908‑4007195

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