Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,567,472 G→C 7.4% intergenic (‑253/+32) yddY ← / ← yddW protein YddY/putative lipoprotein YddW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,567,4720GC 7.4% 347.3 / 14.8 108intergenic (‑253/+32)yddY/yddWprotein YddY/putative lipoprotein YddW
Reads supporting (aligned to +/- strand):  ref base G (48/52);  new base C (4/4);  total (52/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 111 total aligned reads displayed.

TTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCG  >  NC_000913/1567323‑1567621
                                                                                                                                                     |                                                                                                                                                     
tttGGGTGAGAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAg                                                                                                                                                       <  1:626718/150‑1 (MQ=255)
tttCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAg                                                                                                                                                       <  2:939921/150‑1 (MQ=255)
    ggTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCAcc                                                                                                                                                   <  1:24438/150‑1 (MQ=255)
     gTGACAGCGAATATCGTATTGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCg                                                                                                                                                  <  1:960007/150‑1 (MQ=255)
     gTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCg                                                                                                                                                  <  1:509845/150‑1 (MQ=255)
          aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTg                                                                                                                                             <  1:372186/150‑1 (MQ=255)
          aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTg                                                                                                                                             <  1:1557332/150‑1 (MQ=255)
           tcGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                            <  1:637680/149‑1 (MQ=255)
           gCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                            <  1:1108865/150‑1 (MQ=255)
           gCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                            >  2:1293917/1‑150 (MQ=255)
           gCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCAGGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                            <  2:24772/150‑1 (MQ=255)
            cGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGtt                                                                                                                                           <  2:247037/150‑1 (MQ=255)
             gAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTa                                                                                                                                          <  1:707196/150‑1 (MQ=255)
                tatCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACaa                                                                                                                                       >  2:1410195/1‑150 (MQ=255)
                tatCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACaa                                                                                                                                       >  2:437299/1‑150 (MQ=255)
                tatCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTGTTAccc                                                                                                                                       <  2:222961/150‑3 (MQ=39)
                  tCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAg                                                                                                                                     >  1:966526/1‑150 (MQ=255)
                    gTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTa                                                                                                                                   <  2:874357/150‑1 (MQ=255)
                       tGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGa                                                                                                                                >  1:823102/1‑150 (MQ=255)
                       tGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGa                                                                                                                                <  1:635365/150‑1 (MQ=255)
                        ggTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAg                                                                                                                                 <  2:750790/148‑1 (MQ=255)
                        ggTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAg                                                                                                                                 >  1:750790/1‑148 (MQ=255)
                        ggTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGaa                                                                                                                               >  1:879986/1‑150 (MQ=255)
                        ggTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGaa                                                                                                                               >  1:1327249/1‑150 (MQ=255)
                         gTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAAt                                                                                                                              >  2:921259/1‑150 (MQ=255)
                          tttttCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATg                                                                                                                             <  1:82007/150‑1 (MQ=255)
                             ttcatattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGa                                                                                                                          >  1:1516120/1‑150 (MQ=255)
                              tcatattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAg                                                                                                                         <  2:542037/150‑1 (MQ=255)
                              tcatattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAg                                                                                                                         <  1:583464/150‑1 (MQ=255)
                              tcatattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAg                                                                                                                         <  2:654987/150‑1 (MQ=255)
                               catattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGt                                                                                                                        >  1:862436/1‑150 (MQ=255)
                                 tattcataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGttt                                                                                                                      <  1:928772/150‑1 (MQ=255)
                                      ataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTg                                                                                                                 >  1:150992/1‑150 (MQ=255)
                                      ataCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGAACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTg                                                                                                                 >  2:1385886/1‑150 (MQ=255)
                                        aCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGcc                                                                                                               <  2:151969/150‑1 (MQ=255)
                                         cATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGccc                                                                                                              >  1:658079/1‑150 (MQ=255)
                                          aTTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGcccc                                                                                                             <  2:1420335/150‑1 (MQ=255)
                                            ttttATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAg                                                                                                           <  1:921259/150‑1 (MQ=255)
                                                 ttaGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACt                                                                                                      >  2:65798/1‑150 (MQ=255)
                                                   aGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                                    >  1:1022009/1‑150 (MQ=255)
                                                   aGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGAATTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAAAAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                                    >  2:767863/1‑150 (MQ=255)
                                                    gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTg                                                                                                   <  1:784589/150‑1 (MQ=255)
                                                    gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTg                                                                                                   >  1:812374/1‑150 (MQ=255)
                                                     ggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGc                                                                                                  <  1:1337447/150‑1 (MQ=255)
                                                      taTTTATGTCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCa                                                                                                 <  1:732535/149‑1 (MQ=255)
                                                      gATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCa                                                                                                 >  2:1044450/1‑150 (MQ=255)
                                                            tGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCACAGCGACTTTGCAGATAGc                                                                                           <  1:827212/150‑1 (MQ=255)
                                                             ggCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCt                                                                                          >  1:1359556/1‑150 (MQ=255)
                                                              gCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTg                                                                                         >  1:811333/1‑150 (MQ=255)
                                                                 gTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGAcc                                                                                      <  1:780474/150‑1 (MQ=255)
                                                                    tAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCt                                                                                   <  1:778539/150‑1 (MQ=255)
                                                                     aaCTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCtt                                                                                  >  2:12682/1‑150 (MQ=255)
                                                                      aCTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACAGCttg                                                                                 >  1:633124/1‑150 (MQ=255)
                                                                         aaGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCttgttg                                                                              <  1:12682/150‑1 (MQ=255)
                                                                          aGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                             >  2:1598603/1‑150 (MQ=255)
                                                                           gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAg                                                                            >  1:1128957/1‑150 (MQ=255)
                                                                           gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAg                                                                            <  2:639970/150‑1 (MQ=255)
                                                                               ggTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTg                                                                        >  2:504084/1‑150 (MQ=255)
                                                                                gTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGc                                                                       >  1:491716/1‑150 (MQ=255)
                                                                                 ttAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCg                                                                      <  2:1112260/150‑1 (MQ=255)
                                                                                    aTTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGtt                                                                   >  1:71371/1‑150 (MQ=255)
                                                                                    aTTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTATAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGtt                                                                   >  2:839623/1‑150 (MQ=255)
                                                                                         aCTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTc                                                              >  1:353326/1‑150 (MQ=255)
                                                                                          cTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCa                                                             <  2:680166/150‑1 (MQ=255)
                                                                                          cTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTCGTTGAGTCTGCGGTTTATTCa                                                             <  1:1437364/150‑1 (MQ=255)
                                                                                           ttAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAg                                                            <  1:37405/150‑1 (MQ=255)
                                                                                           ttAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAg                                                            <  2:862436/150‑1 (MQ=255)
                                                                                             aagtaagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGAt                                                          <  1:77658/150‑1 (MQ=255)
                                                                                                  agCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTcc                                                     <  1:1175621/150‑1 (MQ=255)
                                                                                                    cATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTc                                                   >  1:805602/1‑150 (MQ=255)
                                                                                                     aTGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCa                                                  >  1:1085868/1‑150 (MQ=255)
                                                                                                      tGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCAc                                                 >  2:215852/1‑150 (MQ=255)
                                                                                                      tGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCAc                                                 >  2:1298407/1‑150 (MQ=255)
                                                                                                            ttaGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACa                                           <  1:1516981/150‑1 (MQ=255)
                                                                                                              aGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATTGTCCTCACGGAACAAg                                         <  2:1432315/150‑1 (MQ=255)
                                                                                                                tGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGAt                                       >  2:943873/1‑150 (MQ=255)
                                                                                                                tGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGAt                                       >  1:1031849/1‑150 (MQ=255)
                                                                                                                     tAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGc                                  >  2:1580452/1‑150 (MQ=255)
                                                                                                                     tAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGc                                  <  1:1598603/150‑1 (MQ=255)
                                                                                                                      aGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGcc                                 <  2:150992/150‑1 (MQ=255)
                                                                                                                      aGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGcc                                 <  2:1516120/150‑1 (MQ=255)
                                                                                                                          tAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTa                             >  2:407856/1‑150 (MQ=255)
                                                                                                                               gggTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGGGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTc                        <  1:1520298/150‑1 (MQ=255)
                                                                                                                                ggTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCg                       <  1:489498/150‑1 (MQ=255)
                                                                                                                                ggTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCg                       >  1:1074195/1‑150 (MQ=255)
                                                                                                                                  tAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCggg                     <  2:879986/150‑1 (MQ=255)
                                                                                                                                   aaCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGc                    <  2:805602/150‑1 (MQ=255)
                                                                                                                                   aaCAAGCCGGTGCTTTACCCACCGGCTTGTTACCAAGTAAGAATGGGAGTTTAACTGCCCCAGCCACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGc                    >  1:681038/1‑150 (MQ=255)
                                                                                                                                    aCAAGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGca                   >  1:617879/1‑150 (MQ=255)
                                                                                                                                     cAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcac                  >  1:811957/1‑150 (MQ=255)
                                                                                                                                      aaGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcaca                 >  2:175198/1‑150 (MQ=255)
                                                                                                                                       aGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAg                <  1:504084/150‑1 (MQ=255)
                                                                                                                                       aGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAg                >  2:1435071/1‑150 (MQ=255)
                                                                                                                                       aGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAg                >  1:1066232/1‑150 (MQ=255)
                                                                                                                                           ggTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACTAGATGGTGCCGCTAATTTCGGGCACAGCATc            >  2:1511742/1‑150 (MQ=255)
                                                                                                                                               ggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTa        <  2:658079/150‑1 (MQ=255)
                                                                                                                                                   aaGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGAt    <  1:4922/150‑1 (MQ=255)
                                                                                                                                                   aaGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGAt    <  2:823102/150‑1 (MQ=255)
                                                                                                                                                    aGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATc   <  1:1108343/150‑1 (MQ=255)
                                                                                                                                                     gCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCg  >  1:1482681/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCGGTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAAGATCG  >  NC_000913/1567323‑1567621

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