Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,168,895 C→T intergenic (‑63/‑178) ycfQ ← / → bhsA repressor for bhsA(ycfR)/biofilm, cell surface and signaling protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,168,8950CT96.2% 272.6 / 1.1 79intergenic (‑63/‑178)ycfQ/bhsArepressor for bhsA(ycfR)/biofilm, cell surface and signaling protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (19/57);  total (20/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGC  >  NC_000913/1168749‑1169042
                                                                                                                                                  |                                                                                                                                                   
ggCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGt                                                                                                                                                 <  2:850665/151‑1 (MQ=255)
     tATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTcc                                                                                                                                               >  2:1113325/1‑148 (MQ=255)
        cAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATcc                                                                                                                                         <  1:91256/151‑1 (MQ=255)
         aaGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATccc                                                                                                                                        <  1:765824/151‑1 (MQ=255)
          aGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTCGTGACTGTGAATGATAGGTTATTATAAAAATAATCACCCCCGTTAATCAGTTCAGATcata                                                                                                                                       <  2:1085704/151‑4 (MQ=255)
             gcgGCATCCCTGTCGAACACTTTTGGGCGACCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATaa                                                                                                                                    <  1:481165/151‑1 (MQ=255)
              cgGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATaa                                                                                                                                    >  2:725478/1‑150 (MQ=255)
               ggCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATaaa                                                                                                                                   >  1:422069/1‑150 (MQ=255)
                  aTCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAaat                                                                                                                                >  1:325246/1‑150 (MQ=255)
                    cccTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAaataa                                                                                                                              <  1:1097712/150‑1 (MQ=255)
                        gTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCt                                                                                                                         <  2:746495/151‑1 (MQ=255)
                          cGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCttt                                                                                                                       <  1:1029454/151‑1 (MQ=255)
                              cacTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTAt                                                                                                                   <  1:483142/151‑1 (MQ=255)
                                  tttGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTAttt                                                                                                                 <  1:499407/149‑1 (MQ=255)
                                   ttGGGCGGCCACGGCTTTTTGTTACACATTGTGTTGAGTCAATTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCGCCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATGTa                                                                                                                <  2:647368/149‑1 (MQ=255)
                                     gggCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTg                                                                                                            <  1:498865/151‑1 (MQ=255)
                                       gcggcCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGaa                                                                                                          <  2:12579/151‑1 (MQ=255)
                                         ggcCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAAt                                                                                                        <  2:172232/151‑1 (MQ=255)
                                               ggCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAg                                                                                                  <  2:956515/151‑1 (MQ=255)
                                                gCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTaa                                                                                                    >  1:323868/1‑148 (MQ=255)
                                                 cTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTaa                                                                                                    >  2:73544/1‑147 (MQ=255)
                                                   tttttttACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGAttt                                                                                              >  1:263275/1‑151 (MQ=255)
                                                     tttttACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGAttttt                                                                                            >  2:619528/1‑151 (MQ=255)
                                                      ttttACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGAttttt                                                                                            <  2:764371/150‑1 (MQ=255)
                                                          acccaTTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCCCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAaat                                                                                        <  1:38264/147‑1 (MQ=255)
                                                           cacaTTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAaatt                                                                                       >  1:244811/1‑150 (MQ=255)
                                                             caTTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAaattaat                                                                                    <  1:582390/151‑1 (MQ=255)
                                                             caTTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAaattaat                                                                                    >  2:1008041/1‑151 (MQ=255)
                                                                tgtgTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAattaatta                                                                                  <  2:1036949/150‑1 (MQ=255)
                                                                tgtgTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAattaatta                                                                                  <  2:745634/150‑1 (MQ=255)
                                                                  tgtTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATtaattaat                                                                                <  2:870391/150‑1 (MQ=255)
                                                                  tgtTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATtaattaat                                                                                <  2:977513/150‑1 (MQ=255)
                                                                     tGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGAtt                                                                            >  1:1011743/1‑151 (MQ=255)
                                                                     tGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGAt                                                                             <  1:352779/150‑1 (MQ=255)
                                                                          cAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTtat                                                                       <  2:682568/151‑1 (MQ=255)
                                                                          cAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTtat                                                                       <  2:821875/151‑1 (MQ=255)
                                                                          cAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTta                                                                        <  2:910371/150‑1 (MQ=255)
                                                                           aGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTtata                                                                      <  2:562588/151‑1 (MQ=255)
                                                                           aGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTtata                                                                      <  1:460844/151‑1 (MQ=255)
                                                                            gTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATaa                                                                     <  2:490581/151‑1 (MQ=255)
                                                                            gTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATaa                                                                     <  1:144635/151‑1 (MQ=255)
                                                                               gCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATaaaa                                                                   <  1:35157/150‑1 (MQ=255)
                                                                                ccATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATaaaaaa                                                                 <  1:613965/151‑1 (MQ=255)
                                                                                  atgatgCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAata                                                               <  2:182365/151‑1 (MQ=255)
                                                                                     atgCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCt                                                            >  2:438804/1‑151 (MQ=255)
                                                                                       gCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTg                                                          <  1:327143/151‑1 (MQ=255)
                                                                                       gCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTg                                                          <  1:226059/151‑1 (MQ=255)
                                                                                               aCCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGa                                                   <  1:742909/150‑1 (MQ=255)
                                                                                               aCCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGa                                                   <  1:69172/150‑1 (MQ=255)
                                                                                               aCCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGAt                                                  <  1:619544/151‑1 (MQ=255)
                                                                                                ccTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATc                                                 <  2:225141/151‑1 (MQ=255)
                                                                                                   ggTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAg                                              <  2:971564/151‑1 (MQ=255)
                                                                                                      gACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGAt                                            <  1:614081/150‑1 (MQ=255)
                                                                                                      gACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGAt                                            <  1:769963/150‑1 (MQ=255)
                                                                                                      gACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGAt                                            <  1:1090684/150‑1 (MQ=255)
                                                                                                         tgtgAATGGAGCGTTATTATAAGAATAATCACCTAGGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATcac                                         >  2:782698/1‑150 (MQ=255)
                                                                                                             aatgaaAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTa                                    <  1:126733/151‑1 (MQ=255)
                                                                                                              atgaaAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTAt                                   >  1:56265/1‑151 (MQ=255)
                                                                                                              atgaaAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTAt                                   >  1:419442/1‑151 (MQ=255)
                                                                                                               tgaaAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATc                                  >  1:1044929/1‑151 (MQ=255)
                                                                                                                gaaAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCa                                 <  2:31378/151‑1 (MQ=255)
                                                                                                                gaaAGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCa                                 <  2:54700/151‑1 (MQ=255)
                                                                                                                  aaGGTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCAtt                               <  1:859275/151‑1 (MQ=255)
                                                                                                                     gTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCAttt                              >  1:898160/1‑149 (MQ=255)
                                                                                                                     gTTATTATAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAg                            <  2:91710/151‑1 (MQ=255)
                                                                                                                      ttattaTAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGt                           >  2:766240/1‑151 (MQ=255)
                                                                                                                       tattatAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGt                           <  1:73545/150‑1 (MQ=255)
                                                                                                                        attatAAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGtt                          <  1:1081367/150‑1 (MQ=255)
                                                                                                                           ataAAAATAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcg                      <  1:634400/151‑1 (MQ=255)
                                                                                                                            taAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcg                      >  1:509377/1‑150 (MQ=255)
                                                                                                                             aaaaaTAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcgat                    <  2:725015/151‑1 (MQ=255)
                                                                                                                             aaaaaTAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcga                     <  2:477317/150‑1 (MQ=255)
                                                                                                                              aaaaTAATCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcgatc                   <  1:158776/151‑1 (MQ=255)
                                                                                                                                aataatCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGt                 <  1:255351/151‑1 (MQ=255)
                                                                                                                                 ataatCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGtt                <  1:485210/151‑1 (MQ=255)
                                                                                                                                   aatCACCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTa               <  1:153529/150‑1 (MQ=255)
                                                                                                                                       aCCTCCGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTaa          <  1:10338/151‑1 (MQ=255)
                                                                                                                                           ccGTTCATCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgc      >  1:11233/1‑151 (MQ=255)
                                                                                                                                             gTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgctt    <  1:15227/151‑1 (MQ=255)
                                                                                                                                               tcatcaGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttgc  >  1:1007502/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCCTTATCAAGCGCGGCATCCCTGTCGAACACTTTTGGGCGGCCACGGCTTTTTTTTACACATTGTGTTGAGTCAGTTGCCATGATGCCGTTGTACCTGGTGACTGTGAATGAAAGGTTATTATAAAAATAATCACCTCCGTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGC  >  NC_000913/1168749‑1169042

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