Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,050 A→T 5.2% intergenic (+18/+78) trpT → / ← hdfR tRNA‑Trp/flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,0500AT 5.2% 833.8 / 24.0 232intergenic (+18/+78)trpT/hdfRtRNA‑Trp/flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (128/90);  new base T (4/8);  total (134/98)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 241 total aligned reads displayed.

TGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCA  >  NC_000913/3946901‑3947198
                                                                                                                                                     |                                                                                                                                                      
tgcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGGCCGTTCCGCCACCCTAATTAGGGGCGTAGTGCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGaa                                                                                                                                                        >  2:1365245/1‑150 (MQ=255)
 gcctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGaaa                                                                                                                                                       >  1:346499/1‑150 (MQ=255)
   ctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCCCCCGAATTAGGGGCGTAGGTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGGTGGGAGTTCGCGTCTCTCCGCCCCGGCCAGAAAGCATCCGTAGCGAAAGc                                                                                                                                                     >  1:1020039/1‑150 (MQ=38)
            aGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑ttt                                                                                                                                          >  2:1949807/1‑150 (MQ=255)
               gggTCGCGGGTTCGAGTGCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCTGAAATCATCCTTAGCGAAAGCTAA‑GGA‑tttttt                                                                                                                                       <  2:1549521/150‑1 (MQ=255)
                ggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCATTAGCTTTCGCTAA‑GGA‑tgatttc                                                                                                                                      <  2:326303/150‑7 (MQ=255)
                ggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCATTAGCTTTCGCTAA‑GGA‑tgatttc                                                                                                                                      <  2:326242/150‑7 (MQ=255)
                   cgcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTAt                                                                                                                                   <  2:2885637/150‑1 (MQ=255)
                    gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTGGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATc                                                                                                                                  >  1:2916661/1‑150 (MQ=255)
                    gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATc                                                                                                                                  <  2:2405209/150‑1 (MQ=255)
                    gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAAATGGTAGAGAACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATAATACATAGCGAAAGCTAA‑AGA‑ATTTTTTTAta                                                                                                                                  >  1:2644651/1‑149 (MQ=255)
                    gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGGAGGTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCGAA‑GGA‑TTTTTTTTATc                                                                                                                                  >  1:357746/1‑150 (MQ=255)
                       ggTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGa                                                                                                                               <  2:455812/150‑1 (MQ=255)
                       ggTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGa                                                                                                                               <  2:570373/150‑1 (MQ=255)
                         ttCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGaaa                                                                                                                             >  1:1906723/1‑150 (MQ=255)
                         ttCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGaaa                                                                                                                             <  2:489829/150‑1 (MQ=255)
                           cGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATa                                                                                                                           >  1:1703905/1‑150 (MQ=255)
                           cGAGTCCCGTCCGTTCCGCCACCCGAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAAACGGGTGTTGGGAGTTCGAGTCACTCCGCCCCTGCCAGAAATCAACCATAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATa                                                                                                                           >  1:1048962/1‑150 (MQ=255)
                             aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAAc                                                                                                                         >  1:1970169/1‑150 (MQ=255)
                             aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAAc                                                                                                                         <  2:1165986/150‑1 (MQ=255)
                             aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAAc                                                                                                                         <  2:1124986/150‑1 (MQ=255)
                               tCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAAccc                                                                                                                       >  2:174003/1‑150 (MQ=255)
                               tCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAAccc                                                                                                                       >  1:2764087/1‑150 (MQ=255)
                                 ccgtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCctc                                                                                                                     >  2:2478627/1‑150 (MQ=255)
                                 ccgtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCctc                                                                                                                     >  1:1114704/1‑150 (MQ=255)
                                   gtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCtctc                                                                                                                   <  1:1775698/150‑1 (MQ=255)
                                    tccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTcc                                                                                                                  <  1:2425576/150‑1 (MQ=255)
                                       gtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGaa                                                                                                               <  2:129009/150‑1 (MQ=255)
                                       gtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGaa                                                                                                               <  2:2134110/150‑1 (MQ=255)
                                       gtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGaa                                                                                                               <  1:142852/150‑1 (MQ=255)
                                          ccGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGCGTTCGAGTCTCTCCGCCCCTGCCAGACATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCGCTCCGAAGTa                                                                                                            >  2:1244636/1‑150 (MQ=255)
                                            gCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTaaa                                                                                                          <  2:942274/150‑1 (MQ=255)
                                                cccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                      >  1:801374/1‑150 (MQ=255)
                                                cccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                      >  1:898364/1‑150 (MQ=255)
                                                 ccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCtt                                                                                                     <  1:604671/150‑1 (MQ=255)
                                                 ccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCtt                                                                                                     >  1:1586222/1‑150 (MQ=255)
                                                        aGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                                                                                              >  1:2028507/1‑150 (MQ=255)
                                                          gggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCCGCa                                                                                            >  1:1826609/1‑150 (MQ=255)
                                                              gTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGATTTTTTTTTTTATTAAATAACCAACACCGAGGTAATTCCTGCTACCGGCATcc                                                                                         >  2:650165/1‑150 (MQ=255)
                                                               tAGTTCAATTTGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGATATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGATGTAAATCCTTCTACCGGCATCCtt                                                                                       <  1:2723130/150‑1 (MQ=255)
                                                                  ttCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcc                                                                                    <  1:217487/150‑1 (MQ=255)
                                                                    cAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccag                                                                                  >  1:1456289/1‑150 (MQ=255)
                                                                     aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGACGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                 >  1:450193/1‑150 (MQ=255)
                                                                     aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                 >  2:2694573/1‑150 (MQ=255)
                                                                     aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                 <  2:1859032/150‑1 (MQ=255)
                                                                     aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                 <  2:329380/150‑1 (MQ=255)
                                                                       ttGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagcca                                                                               <  2:2180541/150‑1 (MQ=255)
                                                                           tAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                           >  1:1084730/1‑150 (MQ=255)
                                                                            agagCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCat                                                                          >  2:146526/1‑150 (MQ=255)
                                                                              agCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                        <  1:1314688/150‑1 (MQ=255)
                                                                              agCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAGTCatat                                                                        >  1:2583156/1‑150 (MQ=255)
                                                                               gCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCCTTCATAtt                                                                       >  1:527651/1‑150 (MQ=255)
                                                                                 aCCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTaa                                                                     <  1:1076021/150‑1 (MQ=255)
                                                                                   cGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTCGCGACAGCTAA‑GGA‑TTTTTTTTATCTGAATTACCCCTCTCCGAAGTAAATCCTTCTACCGGCAGCCTTGCCCGCCATGCCTATTactc                                                                   >  1:707668/1‑147 (MQ=255)
                                                                                     gTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATaca                                                                 <  1:2077/150‑1 (MQ=255)
                                                                                        tcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACtt                                                              >  1:1295134/1‑150 (MQ=255)
                                                                                         ccAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCCTTCATATTAATACACTTc                                                             >  2:366407/1‑150 (MQ=255)
                                                                                            aaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATc                                                          <  2:2789558/150‑1 (MQ=255)
                                                                                            aaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATc                                                          <  1:2698900/150‑1 (MQ=255)
                                                                                             aaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATcc                                                         <  2:2307319/150‑1 (MQ=255)
                                                                                              aCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCa                                                        <  2:750454/150‑1 (MQ=255)
                                                                                                 gggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCa                                                     >  1:382337/1‑150 (MQ=255)
                                                                                                 gggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTGCTACCGGCATCCTTGCCAGCCCTTCATATGAATACACTTCATCCAGcc                                                     >  2:89953/1‑149 (MQ=255)
                                                                                                  ggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCAc                                                    >  1:878766/1‑150 (MQ=255)
                                                                                                   gtgtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACg                                                   >  1:1390017/1‑150 (MQ=255)
                                                                                                   gtgtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAGTCATATTAATACACTTCATCCAGCACg                                                   >  1:2521700/1‑150 (MQ=255)
                                                                                                         ggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAAtt                                             >  1:1863916/1‑150 (MQ=255)
                                                                                                            gTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTc                                          <  2:1168393/150‑1 (MQ=255)
                                                                                                             ttCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCa                                         >  1:628881/1‑150 (MQ=255)
                                                                                                              tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAATCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaa                                        <  1:2897897/150‑1 (MQ=255)
                                                                                                              tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaa                                        <  2:1995317/150‑1 (MQ=255)
                                                                                                                   tctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                   >  1:167494/1‑150 (MQ=255)
                                                                                                                   tctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                   >  1:1563639/1‑150 (MQ=255)
                                                                                                                     tctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcg                                 >  1:1463704/1‑150 (MQ=255)
                                                                                                                      ctcCGCCCCTGCCATTAATAATCCATAGCGTAAGCTAATGTT‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcg                                 <  1:650165/150‑1 (MQ=255)
                                                                                                                       tcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcgcg                               >  1:1515905/1‑150 (MQ=255)
                                                                                                                       tcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGATAATTTTCAAAAGATcgcg                               >  1:309543/1‑150 (MQ=255)
                                                                                                                        ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                              <  1:2772385/150‑1 (MQ=255)
                                                                                                                          gCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAAt                            <  2:2487648/150‑1 (MQ=255)
                                                                                                                           ccccTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATc                           >  1:1532038/1‑150 (MQ=255)
                                                                                                                           ccccTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATc                           >  1:2203959/1‑150 (MQ=255)
                                                                                                                             ccTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCaa                         >  1:400957/1‑150 (MQ=255)
                                                                                                                              cTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAAc                        >  1:2761953/1‑150 (MQ=255)
                                                                                                                              cTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAAc                        <  2:1639942/150‑1 (MQ=255)
                                                                                                                               tGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACg                       >  2:182604/1‑150 (MQ=255)
                                                                                                                                      aaTCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttttt                <  2:1407411/150‑1 (MQ=255)
                                                                                                                                      aaTCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttttt                <  1:1949807/150‑1 (MQ=255)
                                                                                                                                       atcatcCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTa               >  1:2097763/1‑150 (MQ=255)
                                                                                                                                       aaaatcCTTAGCTTTCGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTa               >  2:2549839/4‑150 (MQ=255)
                                                                                                                                       aaaatcCTTAGCTTTCGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTa               >  1:942864/4‑150 (MQ=255)
                                                                                                                                          atcCTTAGCTTTCGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                  >  1:485810/1‑144 (MQ=255)
                                                                                                                                           tcCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGc           >  1:992501/1‑150 (MQ=255)
                                                                                                                                            ccTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCt          >  1:1363836/1‑150 (MQ=255)
                                                                                                                                              ttAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAt        >  1:2337828/1‑150 (MQ=255)
                                                                                                                                              ttAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAt        >  1:2337840/1‑150 (MQ=255)
                                                                                                                                              ttAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCAAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAt        >  2:784801/1‑150 (MQ=255)
                                                                                                                                                aGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAttt      <  1:2254375/150‑1 (MQ=255)
                                                                                                                                                  cGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTg    >  2:2254676/1‑150 (MQ=255)
                                                                                                                                                    aaaGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCa  <  1:340653/150‑1 (MQ=255)
                                                                                                                                                    aaaGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCa  >  2:2794490/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAA‑GGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCA  >  NC_000913/3946901‑3947198

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