Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,142,433 A→C 5.3% intergenic (‑212/‑97) frwA ← / → frwC putative PTS enzyme: Hpr, enzyme I and II components/putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,142,4330AC 5.3% 285.4 / 3.1 96intergenic (‑212/‑97)frwA/frwCputative PTS enzyme: Hpr, enzyme I and II components/putative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base A (37/53);  new base C (4/2);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAATCTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGG  >  NC_000913/4142287‑4142581
                                                                                                                                                  |                                                                                                                                                    
accaagcTGCCAACAGCTGGACAAATGAAATGTAACCGTCAATGTGAGACGCGGCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATgg                                                                                                                                                   <  2:787246/144‑1 (MQ=255)
aCAAATCTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATgg                                                                                                                                                   <  2:740007/150‑1 (MQ=255)
aCAAATCTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATgg                                                                                                                                                   <  2:740012/150‑1 (MQ=255)
aCAAATCTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTACGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATgg                                                                                                                                                   >  1:1061963/1‑150 (MQ=255)
    aTCTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCtt                                                                                                                                               <  2:963329/150‑1 (MQ=255)
      cTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCtttt                                                                                                                                             <  1:1309493/150‑1 (MQ=255)
       tcccAAAAGCTGGCCACATGTAATCTAACCGGCAATTAGCGACGCGTCTCACAAGACGCTGGTTTGCGGCAAGCTGCCGGTTTATCGCAAGTTATGAGGCGGACCGCACTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCttttt                                                                                                                                            <  2:45680/148‑1 (MQ=255)
       tGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCttttt                                                                                                                                            <  1:390621/150‑1 (MQ=255)
         ccAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCATGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGGCCATTTACTGGACttttgta                                                                                                                                          <  2:60714/150‑4 (MQ=255)
         ccAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTAt                                                                                                                                          >  2:99635/1‑150 (MQ=255)
         ccAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTAt                                                                                                                                          >  1:939967/1‑150 (MQ=255)
           aaaaGCTGGACAAATGTAATGTAACCGTCAATTTGCGACTCGTCTCACAAGACGCTGGTTTGCGGCATGCTTCCGGTTTATCTCAAGTTATGAGGCTGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATcc                                                                                                                                        <  1:1182707/150‑1 (MQ=255)
           aaaaGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATcc                                                                                                                                        >  2:650314/1‑150 (MQ=255)
           aaaaGCTGGACAAATGTAATGTAACAGTCAATTTGCGACGCGTCACACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATAATACAAAAGTCCAGTAAATGGCCTTTTTATcc                                                                                                                                        <  1:622430/150‑1 (MQ=255)
            aaaGCTGGACAAAAGTAATGTAACCGTCAATATGCGACCCGTCTCACAAGACGCGGTCTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCTCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCa                                                                                                                                       <  2:59453/150‑1 (MQ=255)
                 tGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGtt                                                                                                                                  <  2:593600/150‑1 (MQ=255)
                 tGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGtt                                                                                                                                  >  1:1155707/1‑150 (MQ=255)
                  tgACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGttt                                                                                                                                 <  1:743338/149‑1 (MQ=255)
                  ggACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGttt                                                                                                                                 >  1:290309/1‑150 (MQ=255)
                     cAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCt                                                                                                                              >  1:378494/1‑150 (MQ=255)
                                  acagTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGTCATGCTTCCGGTATATCGCAATTTATGAGGCTGATCGCAGTTTTGTACTGAGATTACAAACGTCCAGTAAATGGCCTGTTTATCCACTGTTTGCTCCGCCGGCGCTTg                                                                                                                 >  2:212538/4‑150 (MQ=255)
                                  aCCGTCAAGTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTACGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCATATTATCCACTGTTTGCTCCGCCTGCGATTg                                                                                                                 <  1:806978/150‑1 (MQ=255)
                                   ccGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGc                                                                                                                <  2:534928/150‑1 (MQ=255)
                                   ccGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGc                                                                                                                >  1:857980/1‑150 (MQ=255)
                                    cGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGcc                                                                                                               <  2:122923/150‑1 (MQ=255)
                                    cGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGACTTTTTATCAACTGTTTTCTCCGCCTGCGATTTcc                                                                                                               >  2:330159/1‑150 (MQ=255)
                                    cGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGGTTCTCGCAAGTGATGAGGCGGCTCGCATTTTTGTCCTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGcc                                                                                                               <  1:629538/150‑1 (MQ=255)
                                    cGTCAATTTGCGACGCGTCTCACAAGACGCTGGTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGcc                                                                                                               <  1:572680/150‑1 (MQ=255)
                                    cGTCAATTTGCGACGCGTCTCACAAGACGCTGGTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGcc                                                                                                               <  1:648570/150‑1 (MQ=255)
                                     gTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCt                                                                                                              >  1:1210654/1‑150 (MQ=255)
                                      tCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGCGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTAGCTCCGCCTGCGATTGCCTa                                                                                                             <  2:671403/150‑1 (MQ=255)
                                            tGCGACGCGTCTCACAAGTCGCTGTTTTGCGGCATGCTTCCGGTTTAACGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGttt                                                                                                       <  2:924635/150‑1 (MQ=255)
                                              cGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTtcc                                                                                                     >  1:409962/1‑150 (MQ=255)
                                                 cgcgTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTtcctcc                                                                                                  <  2:463420/150‑1 (MQ=255)
                                                  gcgTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCa                                                                                                 >  2:638724/1‑150 (MQ=255)
                                                    gTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGCGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGAAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATc                                                                                               >  2:1035404/1‑150 (MQ=255)
                                                    gTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCTCATGCGATTGCCTATTGTTTCCTCCATc                                                                                               >  1:693822/1‑150 (MQ=255)
                                                       tcACAAGACGCTGTTTTGCGGTATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCaaa                                                                                            <  1:1179769/150‑1 (MQ=255)
                                                       tcACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCaaa                                                                                            >  2:399275/1‑150 (MQ=255)
                                                           aaGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAata                                                                                        >  2:461157/1‑150 (MQ=255)
                                                             gACGCTGTTTTGCGGCATGCTTCGGGTTTATCGCAAGTTATGAGTCGGATCGCATTTATGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATg                                                                                      <  2:670933/150‑1 (MQ=255)
                                                             gACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATg                                                                                      >  1:1036720/1‑150 (MQ=255)
                                                               cGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGa                                                                                    >  2:210266/1‑150 (MQ=255)
                                                                gCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGAc                                                                                   >  1:904906/1‑150 (MQ=255)
                                                                   gTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATg                                                                                >  2:600946/1‑150 (MQ=255)
                                                                     tttGCGGCATGCTTCCGGTTTATCTCAAGTTATGAGGCGGATCGCATTTTAGTACTGATATTACAAAAGTCCACTAAATTGCCTTTTTATCGACTGTTTGCTCCGCCTGCGATTTCCTATTGTTTCCTCCATCAAAAATATGGACATggg                                                                              <  1:805960/150‑1 (MQ=255)
                                                                      ttGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGc                                                                             <  1:499906/150‑1 (MQ=255)
                                                                       tGCGGCATGCTGCCGGTTTATCGCAAGTTATGAGGCGTATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCt                                                                            <  1:526154/150‑1 (MQ=255)
                                                                         cGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTTCTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCCCTGTTTGCTCCGCCTGCGATTGACTATTGTTTCCTCCATCAAAAATATGGACATGGGCttt                                                                          >  2:1039386/1‑150 (MQ=255)
                                                                         cGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCttt                                                                          <  1:1158237/150‑1 (MQ=255)
                                                                         cGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCttt                                                                          <  2:576322/150‑1 (MQ=255)
                                                                          ggCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTc                                                                         >  1:557388/1‑150 (MQ=255)
                                                                           gcctGCTTCCGGGTACTCGCAAGTTTTGCGGCGGTCCGCATTTTGGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTAGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATTGGCTATct                                                                        <  2:642172/147‑2 (MQ=255)
                                                                           gCATGCTTGCTGTTTATCGCAAGTTATGAGGCTGATCGCATTTTTGTACTGATTTTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTTCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCg                                                                        <  2:554591/150‑1 (MQ=255)
                                                                             aTGCTTCCTGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCa                                                                      <  2:1061963/150‑1 (MQ=255)
                                                                             aTGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCa                                                                      <  1:546668/150‑1 (MQ=255)
                                                                                cTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGGCATGGGCTTTCGCAAGc                                                                   >  2:1176592/1‑150 (MQ=255)
                                                                                   ccGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCg                                                                <  2:68375/150‑1 (MQ=255)
                                                                                     ggattATCGCAAGTATTGCGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGGCTTTTTATCCACTGTTTGCTCCTCCTGCGATTGCCTCTTGTTTCCTCCATCAAAAATATGGACATGGGCTTGCGAAAGcccatg                                                              <  2:193239/147‑4 (MQ=255)
                                                                                     ggTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCGTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCggg                                                              <  1:901546/150‑1 (MQ=255)
                                                                                        ttATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGagc                                                           <  1:1148391/150‑1 (MQ=255)
                                                                                        ttATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGagc                                                           <  2:939967/150‑1 (MQ=255)
                                                                                          aTCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGagcag                                                         >  2:275908/1‑150 (MQ=255)
                                                                                            cGCAAGTTCTGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCt                                                       <  1:1091686/150‑1 (MQ=255)
                                                                                            cGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCt                                                       <  1:1231538/150‑1 (MQ=255)
                                                                                                aGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGAAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTatg                                                   <  2:1157233/150‑1 (MQ=255)
                                                                                                        ggCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGtt                                           <  1:1103309/150‑1 (MQ=255)
                                                                                                          cGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCAATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTgg                                         <  1:1080038/150‑1 (MQ=255)
                                                                                                           ggATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTGTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGt                                        <  1:148988/150‑1 (MQ=255)
                                                                                                           ggATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATACACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGt                                        <  1:1277703/150‑1 (MQ=255)
                                                                                                               ggCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAg                                    <  1:1033651/149‑1 (MQ=255)
                                                                                                                 cATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGAt                                  >  1:374252/1‑150 (MQ=255)
                                                                                                                  aTTTTTGTACTGATATTACAAAAGTCCATTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATc                                 <  2:500108/150‑1 (MQ=255)
                                                                                                                  aTTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATc                                 <  2:298745/150‑1 (MQ=255)
                                                                                                                   tttttGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGGTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTAAGAATGAGTTGGTGCAGATcc                                >  1:1057188/1‑150 (MQ=255)
                                                                                                                    ttttGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCt                               <  1:602769/150‑1 (MQ=255)
                                                                                                                      ttGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGa                             <  1:331686/150‑1 (MQ=255)
                                                                                                                       tGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGaa                            >  2:1024953/1‑150 (MQ=255)
                                                                                                                       agTACTGATATTACAAAAGAACACTAAATGGCCTTTTTATCCACTGTTTGCTCAGCCTCCGATTGCCAATTGTTTCCTCCATCAAAAATATGGACCTGGGCTTTCTCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCATATCCTGaa                            >  2:688125/2‑150 (MQ=255)
                                                                                                                        gTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGGGCAGCTTATGAATGAGTTGGTGCAGATCCTGaaa                           >  2:813857/1‑150 (MQ=255)
                                                                                                                          aCTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGaaaaa                         >  1:1091908/1‑150 (MQ=255)
                                                                                                                          aCTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGaaaaa                         >  2:1254241/1‑150 (MQ=255)
                                                                                                                              ttatTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATAcc                     <  1:602982/149‑1 (MQ=255)
                                                                                                                                ggataaaaaGGCCATTTACTGGACTTTTTATCCACTGTTTGCTCCGCCTGAGATAGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCg                   >  2:711042/6‑150 (MQ=255)
                                                                                                                                  ataaaaaGGCCATTTACTGGACTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGATTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTc                 >  1:925269/4‑150 (MQ=255)
                                                                                                                                    cAACAGTCCACTAACTGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTAATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGGCAg               <  2:496328/150‑1 (MQ=255)
                                                                                                                                    cAAAACTCCACTAAATGGCCATTTAATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAg               <  1:815926/150‑1 (MQ=255)
                                                                                                                                     aaaaGGCCATTTACTGGACTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGc              >  1:566951/1‑150 (MQ=255)
                                                                                                                                     aaaaGGCCATTTACTGGACTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATACTGAAAAATACCCGTCAGc              >  1:595632/1‑150 (MQ=255)
                                                                                                                                      aaaGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCa             <  1:1203175/150‑1 (MQ=255)
                                                                                                                                         gTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGTGCAGCTTATGAATGAGTTGGTGCAGCTCCTGAAAAATACCCGTCAGCAttt          >  2:573525/1‑150 (MQ=255)
                                                                                                                                         gTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCAttt          >  1:624075/1‑150 (MQ=255)
                                                                                                                                          tCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGACTGCGATTGCCTATTGTTTCCTACATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGAGCAGCTCCTGAAAAAAACCCGTCAGCATTTa         >  2:766786/1‑150 (MQ=255)
                                                                                                                                             aGTAAATGGCCTTTTTATCCTCTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCATGAAAAATACCCGTCAGCATTTAATg      >  2:974068/1‑150 (MQ=255)
                                                                                                                                             aGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCACGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCAGTCAGCATTTAATg      >  1:100372/1‑150 (MQ=255)
                                                                                                                                              gTAAATGGCCTTTTCATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAATAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGa     <  1:290/150‑1 (MQ=255)
                                                                                                                                                 aaTGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATTGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATAACCGTCAGCATTTAATGACgg  >  1:513769/1‑150 (MQ=255)
                                                                                                                                                 aaTGGCCTTTTTATCCACTGTTTGATCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGAATTTAATGACgg  >  2:739746/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACAAATCTGCCAAAAGCTGGACAAATGTAATGTAACCGTCAATTTGCGACGCGTCTCACAAGACGCTGTTTTGCGGCATGCTTCCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGG  >  NC_000913/4142287‑4142581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: