Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,362,528 G→A 5.3% intergenic (‑175/+23) yjdQ ← / ← pheU pseudogene, P4‑like integrase remnant/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,362,5280GA 5.3% 317.3 / 5.1 95intergenic (‑175/+23)yjdQ/pheUpseudogene, P4‑like integrase remnant/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base A (2/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

AGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAG  >  NC_000913/4362380‑4362676
                                                                                                                                                    |                                                                                                                                                    
aGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCg                                                                                                                                                      >  3:319664/1‑149 (MQ=255)
  tGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTg                                                                                                                                                    >  1:282265/1‑149 (MQ=255)
    cTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGa                                                                                                                                                  <  2:539529/149‑1 (MQ=255)
    cTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGa                                                                                                                                                  >  2:387126/1‑149 (MQ=255)
     tataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAg                                                                                                                                                 <  2:865796/149‑1 (MQ=255)
     tataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAg                                                                                                                                                 <  1:778263/149‑1 (MQ=255)
      ataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAgg                                                                                                                                                <  4:815242/149‑1 (MQ=255)
      ataGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGCTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAgg                                                                                                                                                >  3:208224/1‑149 (MQ=255)
         gCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTcc                                                                                                                                             <  2:528657/149‑1 (MQ=255)
         gCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTcc                                                                                                                                             <  4:362585/149‑1 (MQ=255)
          cTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAACGGACCTCCACGGAGGTCCGt                                                                                                                                           <  4:786899/149‑1 (MQ=255)
           tCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGt                                                                                                                                           <  4:72563/149‑1 (MQ=255)
            cAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTATATCGTGGCGTCCTGAAACGAAAAACGGTCCTCCGTGGAGGTCCGta                                                                                                                                          <  2:745472/149‑2 (MQ=255)
            cAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGtt                                                                                                                                          <  2:104963/149‑1 (MQ=255)
            cAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCAGTGGAGGTCCGtt                                                                                                                                          >  2:578432/1‑149 (MQ=255)
             aGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTACATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGttt                                                                                                                                         >  3:731657/1‑149 (MQ=255)
             aGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGttt                                                                                                                                         <  2:778919/149‑1 (MQ=255)
                aaCAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTtata                                                                                                                                      <  1:397421/149‑1 (MQ=255)
                aaCAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTtata                                                                                                                                      <  2:125169/149‑1 (MQ=255)
                 aCAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTatat                                                                                                                                     <  4:454965/149‑1 (MQ=255)
                   aaaTTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGa                                                                                                                                   >  2:239329/1‑149 (MQ=255)
                   aaaTTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGa                                                                                                                                   >  1:338609/1‑149 (MQ=255)
                     aatttCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAAt                                                                                                                                 >  3:553343/3‑149 (MQ=255)
                     aTTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAAt                                                                                                                                 <  4:319664/149‑1 (MQ=255)
                     aTTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAAt                                                                                                                                 <  2:177615/149‑1 (MQ=255)
                       tttCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAAttt                                                                                                                               <  2:111164/149‑1 (MQ=255)
                        ttCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTg                                                                                                                              <  3:679922/149‑1 (MQ=255)
                        ttCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTg                                                                                                                              <  3:411370/149‑1 (MQ=255)
                           ccGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTgg                                                                                                                                                   <  1:638329/125‑1 (MQ=255)
                           ccGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTgg                                                                                                                                                   >  2:638329/1‑125 (MQ=255)
                              aaTTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGccc                                                                                                                        <  3:810649/149‑1 (MQ=255)
                                        tGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGaa                                                                                                              <  3:182381/149‑1 (MQ=255)
                                           ccGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCg                                                                                                           >  2:768291/1‑149 (MQ=255)
                                            cGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGa                                                                                                          >  4:496591/1‑149 (MQ=255)
                                                          gATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAg                                                                                                    >  3:230923/1‑141 (MQ=255)
                                                          gATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAg                                                                                                    <  4:230923/141‑1 (MQ=255)
                                                           aTGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACgggg                                                                                           >  2:569968/1‑149 (MQ=255)
                                                              tgGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAtt                                                                                        <  2:338609/149‑1 (MQ=255)
                                                               ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAttt                                                                                       <  1:244499/149‑1 (MQ=255)
                                                                gggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAtttt                                                                                      >  2:262928/1‑149 (MQ=255)
                                                                gggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAtttt                                                                                      <  1:239329/149‑1 (MQ=255)
                                                                 ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTc                                                                                     >  4:60042/1‑149 (MQ=255)
                                                                 ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTc                                                                                     <  4:22516/149‑1 (MQ=255)
                                                                        cTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATcccc                                                                              <  2:612263/149‑1 (MQ=255)
                                                                          ggTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGAGCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTg                                                                            <  2:774962/149‑1 (MQ=255)
                                                                           gctGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGc                                                                           <  4:208224/147‑1 (MQ=255)
                                                                                   cGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGa                                                                   <  2:422890/149‑1 (MQ=255)
                                                                                     gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAAcc                                                                                                       <  1:27940/111‑1 (MQ=255)
                                                                                     gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAAcc                                                                                                       >  2:27940/1‑111 (MQ=255)
                                                                                     gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACt                                                                 >  2:384340/1‑149 (MQ=255)
                                                                                     gTAGAATTGACTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGAAACGGGGATTTTCAATCCCCTGCTCTACCGACt                                                                 <  4:553343/149‑1 (MQ=255)
                                                                                      tAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACATCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTg                                                                >  1:752563/1‑149 (MQ=255)
                                                                                         aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGc                                                             >  1:680539/1‑149 (MQ=255)
                                                                                          aTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCt                                                            >  3:354105/1‑149 (MQ=255)
                                                                                           ttGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTa                                                           >  1:378393/1‑149 (MQ=255)
                                                                                           ttGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTa                                                           >  4:348956/1‑149 (MQ=255)
                                                                                           ttGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTa                                                           >  2:119576/1‑149 (MQ=255)
                                                                                               tttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCg                                                       >  2:81382/1‑149 (MQ=255)
                                                                                               tttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCg                                                       >  3:104800/1‑149 (MQ=255)
                                                                                                      gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACg                                                >  1:720783/1‑149 (MQ=255)
                                                                                                        aGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACggg                                              >  2:380428/1‑149 (MQ=255)
                                                                                                          accGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGGATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGAATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGgc                                            >  2:304346/3‑149 (MQ=255)
                                                                                                          aaCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCACGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGgc                                            <  3:245698/149‑1 (MQ=255)
                                                                                                               ttAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                       >  1:303638/1‑149 (MQ=255)
                                                                                                               ttAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                       >  1:507145/1‑149 (MQ=255)
                                                                                                               ttAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAtt                                       >  4:680867/1‑149 (MQ=255)
                                                                                                                   aTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAc                                   >  4:50531/1‑149 (MQ=255)
                                                                                                                   aTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAc                                   <  4:774059/149‑1 (MQ=255)
                                                                                                                    tCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAcc                                  >  4:138433/1‑149 (MQ=255)
                                                                                                                         gCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAAt                             <  3:252522/149‑1 (MQ=255)
                                                                                                                          cGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATc                            <  1:316316/149‑1 (MQ=255)
                                                                                                                            tCCTGAAACGAGAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACCCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCAACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTg                          >  1:716350/1‑149 (MQ=255)
                                                                                                                            tCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                  >  3:386958/1‑109 (MQ=255)
                                                                                                                            tCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                  <  4:386958/109‑1 (MQ=255)
                                                                                                                              cTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGc                                                    >  1:477796/1‑121 (MQ=255)
                                                                                                                              cTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGc                                                    <  2:477796/121‑1 (MQ=255)
                                                                                                                                    cGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAtctc                  >  3:197488/1‑149 (MQ=255)
                                                                                                                                     tataaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                 >  3:127009/4‑149 (MQ=255)
                                                                                                                                     gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                 >  1:206381/1‑149 (MQ=255)
                                                                                                                                      aaaaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGt                <  1:304346/149‑1 (MQ=255)
                                                                                                                                       aaaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTc               >  2:129994/1‑149 (MQ=255)
                                                                                                                                        aaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCa              >  2:167873/1‑149 (MQ=255)
                                                                                                                                         aaCGGTCCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAc                                   <  3:157263/127‑1 (MQ=255)
                                                                                                                                         aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAc                                   >  4:157263/1‑127 (MQ=255)
                                                                                                                                         aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCaa             >  1:299261/1‑149 (MQ=255)
                                                                                                                                         aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCaa             >  1:872064/1‑149 (MQ=255)
                                                                                                                                          aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc            >  4:452516/1‑149 (MQ=255)
                                                                                                                                          aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAAc            <  3:786161/149‑1 (MQ=255)
                                                                                                                                             gACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAt                                                                                         <  1:47033/69‑1 (MQ=21)
                                                                                                                                             gACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAt                                                                                         >  2:47033/1‑69 (MQ=21)
                                                                                                                                                cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAtt      >  2:682463/1‑149 (MQ=255)
                                                                                                                                                cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAAGCTGCACATCTCGTCAACCTAAtt      >  1:712526/1‑149 (MQ=255)
                                                                                                                                                   cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa   <  4:347183/149‑1 (MQ=255)
                                                                                                                                                   cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa   <  2:282265/149‑1 (MQ=255)
                                                                                                                                                   cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa   <  1:174090/149‑1 (MQ=255)
                                                                                                                                                   cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa   <  2:507145/149‑1 (MQ=255)
                                                                                                                                                    gTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAg  <  2:680539/149‑1 (MQ=255)
                                                                                                                                                    gTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAg  <  1:768291/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAG  >  NC_000913/4362380‑4362676

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: