Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,780,128 T→C 5.9% I101T (ATT→ACT)  lldD → L‑lactate dehydrogenase, FMN‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,780,1280TC 5.9% 376.4 / 7.2 114I101T (ATT→ACT) lldDL‑lactate dehydrogenase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base T (50/57);  new base C (3/4);  total (53/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TGAAAAACATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGGTTGTTCG  >  NC_000913/3779981‑3780276
                                                                                                                                                   |                                                                                                                                                    
tGAAAAACATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTAtt                                                                                                                                                     >  2:281506/1‑149 (MQ=255)
 gAAAAACATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTc                                                                                                                                                    <  3:901725/149‑1 (MQ=255)
  aaaaaTATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGGTGCTGGTGGTACTGACTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGCAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTcc                                                                                                                                                   <  2:497522/149‑1 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGt                                                                                                                                                 <  1:448252/149‑1 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGt                                                                                                                                                 >  3:451503/1‑149 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGt                                                                                                                                                 <  2:782772/149‑1 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGt                                                                                                                                                 >  1:782772/1‑149 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGt                                                                                                                                                 >  3:856782/1‑149 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGt                                                                                                                                                 <  4:856782/149‑1 (MQ=255)
    aaaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGt                                                                                                                                                 <  4:451503/149‑1 (MQ=255)
     aaCATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGtt                                                                                                                                                <  3:923509/149‑1 (MQ=255)
        aTGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTAc                                                                                                                                             >  2:512222/1‑149 (MQ=255)
           tCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTActct                                                                                                                                          >  3:830082/1‑149 (MQ=255)
             cGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTActctct                                                                                                                                        >  1:702938/1‑149 (MQ=255)
                cTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGa                                                                                                                                     >  3:360315/1‑149 (MQ=255)
                 ttAAGCCTGGAAACGACGCTGCTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGAGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGCTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGAc                                                                                                                                    >  4:555819/1‑149 (MQ=255)
                       cTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTc                                                                                                                              >  1:935933/1‑149 (MQ=255)
                          gAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGt                                                                                                                           >  3:59654/1‑149 (MQ=255)
                           aaaCGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGtt                                                                                                                          <  4:553911/149‑1 (MQ=255)
                           aaaCGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGtt                                                                                                                          <  3:313545/149‑1 (MQ=255)
                            aaCGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGttt                                                                                                                         <  3:169090/149‑1 (MQ=255)
                            aaCGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGttt                                                                                                                         <  2:156729/149‑1 (MQ=255)
                                  gCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGa                                                                                                                   >  4:718308/1‑149 (MQ=255)
                                   cTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGAt                                                                                                                  <  1:837286/149‑1 (MQ=255)
                                   cTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGAt                                                                                                                  <  2:91317/149‑1 (MQ=255)
                                   cTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGAt                                                                                                                  <  4:449118/149‑1 (MQ=255)
                                    tGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGAtt                                                                                                                 <  1:920240/149‑1 (MQ=255)
                                      tttAATGAGAACTTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCACGGTATTCCGTTTACTCTCTCGACGGTTTcc                                                                                                                             <  4:911447/135‑1 (MQ=255)
                                      tttAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTcc                                                                                                                             >  3:911447/1‑135 (MQ=255)
                                       ttAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              >  4:863098/1‑149 (MQ=255)
                                       ttAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              <  4:830082/149‑1 (MQ=255)
                                       ttAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              <  2:826715/149‑1 (MQ=255)
                                       ttAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              <  2:702938/149‑1 (MQ=255)
                                          aTGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaagaa                                                                                                           <  2:467056/149‑1 (MQ=255)
                                          aTGAGAAATTGTAGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaagaa                                                                                                           <  1:86024/149‑1 (MQ=255)
                                           tGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGaagaag                                                                                                          >  1:41012/1‑149 (MQ=255)
                                               aaaTTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTcgc                                                                                                      >  1:471599/1‑149 (MQ=255)
                                               aaaTTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTGCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCAGTTTACTCTCTAGACGGTTTCCGTTTGCCGGATTGAAGTAGTcgc                                                                                                      >  1:724429/1‑149 (MQ=255)
                                                  ttGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGa                                                                                                                                     >  3:136129/1‑115 (MQ=255)
                                                  ttGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGa                                                                                                                                     <  4:136129/115‑1 (MQ=255)
                                                  ttGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTATGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCgcc                                                                                                   <  4:602445/149‑1 (MQ=255)
                                                    gTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCgccag                                                                                                 >  2:800942/1‑149 (MQ=255)
                                                     tCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCgccagc                                                                                                >  4:49155/1‑149 (MQ=255)
                                                       gATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCgccagcca                                                                                              >  2:663087/1‑149 (MQ=255)
                                                       gATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGCCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCgccagcca                                                                                              >  3:148015/1‑149 (MQ=255)
                                                           ccGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCaa                                                                                          <  1:116962/149‑1 (MQ=255)
                                                                  cACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCa                                                                                   <  2:338878/149‑1 (MQ=255)
                                                                   aCTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGaagaag                                                                                                          <  3:537430/125‑1 (MQ=255)
                                                                    cTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAAt                                                                                 >  1:817692/1‑149 (MQ=255)
                                                                    cTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAAt                                                                                 >  4:235686/1‑149 (MQ=255)
                                                                    cTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTACGTTTACTCTCTCGACGGTTTCCGATTGCCCGAtt                                                                                                                 <  3:213708/117‑1 (MQ=255)
                                                                     tGGTTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              <  4:610137/119‑1 (MQ=255)
                                                                     tGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAAtg                                                                                <  3:235686/149‑1 (MQ=255)
                                                                         tCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGttt                                                                                                                         <  4:218060/104‑1 (MQ=255)
                                                                         tCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGt                                                                            >  3:198607/1‑149 (MQ=255)
                                                                             gTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCa                                                                        <  4:702880/149‑1 (MQ=255)
                                                                             gTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCa                                                                        <  2:79574/149‑1 (MQ=255)
                                                                              tGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAg                                                                       >  2:483739/1‑149 (MQ=255)
                                                                              tGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAg                                                                       <  4:198607/149‑1 (MQ=255)
                                                                                ggcttGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCt                                                                     >  3:498775/4‑149 (MQ=255)
                                                                                ggTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCt                                                                     <  4:74794/149‑1 (MQ=255)
                                                                                 gTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTcc                                                                                                                             >  4:663845/1‑92 (MQ=255)
                                                                                 gTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTcc                                                                                                                             <  3:663845/92‑1 (MQ=255)
                                                                                  tttGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCttt                                                                   <  4:498775/149‑1 (MQ=255)
                                                                                  tttGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCttt                                                                   <  1:756516/149‑1 (MQ=255)
                                                                                    tgtgtgGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTACTCCGTTTACTCTCTCGACGCTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTCCCAGCTTTAt                                                                 <  3:569302/149‑1 (MQ=255)
                                                                                        tgGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTAc                                                             <  1:663087/149‑1 (MQ=255)
                                                                                             tGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTgcgc                                                        <  2:521392/149‑1 (MQ=255)
                                                                                             tGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGCTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTgcgc                                                        >  3:782935/1‑149 (MQ=255)
                                                                                                      gtcgtGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGc                                               <  4:59654/149‑1 (MQ=255)
                                                                                                         gtGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCttt                                            <  1:281506/149‑1 (MQ=255)
                                                                                                            gCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATg                                         <  4:493884/149‑1 (MQ=255)
                                                                                                            gCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATg                                         >  4:466928/1‑149 (MQ=255)
                                                                                                             cGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGc                                        <  4:360315/149‑1 (MQ=255)
                                                                                                              gAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCg                                       >  1:576525/1‑149 (MQ=255)
                                                                                                              gAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCg                                       >  2:176814/1‑149 (MQ=255)
                                                                                                               aaGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTAt                                          >  4:7846/1‑145 (MQ=255)
                                                                                                                 gTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTaa                                    >  1:522258/1‑149 (MQ=255)
                                                                                                                 gTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTaa                                    >  3:190161/1‑149 (MQ=255)
                                                                                                                     aGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              >  3:182558/1‑71 (MQ=255)
                                                                                                                     aGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTgaa                                                                                                              <  4:182558/71‑1 (MQ=255)
                                                                                                                     aGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAAcgcg                                >  1:555736/1‑149 (MQ=255)
                                                                                                                       gCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCt                              <  1:662909/149‑1 (MQ=255)
                                                                                                                       gCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCt                              <  3:49155/149‑1 (MQ=255)
                                                                                                                           ccAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGgag                          <  3:688605/149‑1 (MQ=255)
                                                                                                                             aaaGCGGCGGATGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGgagcg                        >  3:598547/1‑149 (MQ=255)
                                                                                                                             aaaGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGgagcg                        >  3:632746/1‑149 (MQ=255)
                                                                                                                              aaGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGgagcga                       >  1:841361/1‑149 (MQ=255)
                                                                                                                              aaGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGgagcga                       >  4:750042/1‑149 (MQ=255)
                                                                                                                                         gcgcATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAagcagc            <  2:935933/149‑1 (MQ=255)
                                                                                                                                         gcgcATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAagcagc            <  4:190161/149‑1 (MQ=255)
                                                                                                                                             aTGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGgt        >  3:781932/1‑149 (MQ=255)
                                                                                                                                              tGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGgtt       >  4:344896/1‑149 (MQ=255)
                                                                                                                                                gTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGgttgt     >  1:530976/1‑149 (MQ=255)
                                                                                                                                                 tATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGgttgtt    <  3:718308/149‑1 (MQ=255)
                                                                                                                                                 tATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGgttgtt    >  1:559773/1‑149 (MQ=255)
                                                                                                                                                  aTTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGGTTGTTc   <  1:483739/149‑1 (MQ=255)
                                                                                                                                                   ttCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAagcagc            >  4:144854/1‑139 (MQ=255)
                                                                                                                                                   ttCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAagcagc            <  3:144854/139‑1 (MQ=255)
                                                                                                                                                   ttCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGGTTGTTCg  >  4:456595/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGAAAAACATGTCCGACTTAAGCCTGGAAACGACGCTGTTTAATGAGAAATTGTCGATGCCGGTGGCACTGGCTCCGGTGGGTTTGTGTGGCATGTATGCGCGTCGTGGCGAAGTTCAGGCAGCCAAAGCGGCGGACGCGCATGGTATTCCGTTTACTCTCTCGACGGTTTCCGTTTGCCCGATTGAAGAAGTCGCGCCAGCCATCAAGCGCCCAATGTGGTTCCAGCTTTATGTACTGCGCGATCGCGGCTTTATGCGTAACGCGCTGGAGCGAGCAAAAGCAGCGGGTTGTTCG  >  NC_000913/3779981‑3780276

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