Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,936,278 A→G 5.5% M1M (ATG→GTG) † rbsB → D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,936,2780AG 5.5% 231.2 / 3.2 73M1M (ATG→GTG) †rbsBD‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

GCATTAATTCTTGGCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCA  >  NC_000913/3936131‑3936423
                                                                                                                                                   |                                                                                                                                                 
gcATTAATTCTTGGCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAatat                                                                                                                                                  <  1:389751/149‑1 (MQ=255)
   ttAATTCTTGGCTTCCTTAATAATTGATTGAATTTGTTAGGTGTTTACTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGaa                                                                                                                                               >  5:192003/1‑149 (MQ=255)
     aaTTCTTGGCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACa                                                                                                                                             >  8:224582/1‑149 (MQ=255)
      aTTCTTGGCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACAt                                                                                                                                            >  3:8527/1‑149 (MQ=255)
        tCTTGGCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATg                                                                                                                                           >  6:400400/1‑148 (MQ=255)
             gCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGaaaaaa                                                                                                                                     <  1:202366/149‑1 (MQ=255)
               ttCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACt                                                                                                                                   >  7:363387/1‑149 (MQ=255)
                tCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTg                                                                                                                                  >  2:111738/1‑149 (MQ=255)
                 ccTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACt                                                                                                                                   <  6:315707/147‑1 (MQ=255)
                 ccTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACt                                                                                                                                   <  5:301294/147‑1 (MQ=255)
                  cTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGc                                                                                                                                >  5:227567/1‑149 (MQ=255)
                   ttAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAGGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGc                                                                                                                                >  6:209253/1‑148 (MQ=255)
                   ttAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGc                                                                                                                                >  5:202483/1‑148 (MQ=255)
                         atGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTg                                                                                                                         <  4:8527/149‑1 (MQ=255)
                              ttGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTc                                                                                                                    <  2:261513/149‑1 (MQ=255)
                               tGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTcc                                                                                                                   >  2:233436/1‑149 (MQ=255)
                               tGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTcc                                                                                                                   >  3:178703/1‑149 (MQ=255)
                                  aTTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCt                                                                                                                <  7:224582/149‑1 (MQ=255)
                                           gtgtTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCt                                                                                                        >  2:211642/1‑148 (MQ=255)
                                             gtTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCg                                                                                                                  >  2:178119/1‑136 (MQ=255)
                                             gtTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCg                                                                                                                  <  1:178119/136‑1 (MQ=255)
                                                     ctATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGATGGCGGTGCTGGTATACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTTCCCTGGGTTCCGCTGTTGCGCTAATCGCCACCg                                                                                             >  8:183310/1‑149 (MQ=255)
                                                      tATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGt                                                                                            >  5:232397/1‑149 (MQ=255)
                                                      tATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGt                                                                                            >  6:392593/1‑149 (MQ=255)
                                                        ttACCAGATGATCGGCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCa                                                                                          <  2:389231/149‑1 (MQ=255)
                                                        ttACCAGATGATAGTCAAAGCGGTTGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCa                                                                                          <  2:297976/149‑1 (MQ=255)
                                                            cAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGc                                                                                      <  1:208563/149‑1 (MQ=255)
                                                              gatgatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCg                                                                                     >  7:135/1‑148 (MQ=255)
                                                              gatgatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCg                                                                                     >  6:60398/1‑148 (MQ=255)
                                                              gatgatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCg                                                                                     >  3:207343/1‑148 (MQ=255)
                                                               atgatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGaa                                                                                   >  1:29198/1‑149 (MQ=255)
                                                               atgatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGaa                                                                                   <  2:29198/149‑1 (MQ=255)
                                                                tgatCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAACGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAAt                                                                                  <  5:209253/149‑1 (MQ=255)
                                                                     gTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGa                                                                              <  1:211642/148‑1 (MQ=255)
                                                                            cGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATTAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCaaaa                                                                      <  7:183310/149‑1 (MQ=255)
                                                                              gtggtgATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGa                                                                    <  6:227567/149‑1 (MQ=255)
                                                                               tggtgATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGa                                                                    >  7:137085/1‑148 (MQ=255)
                                                                                 gtgATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGacac                                                                 <  1:111738/149‑1 (MQ=255)
                                                                                 gtgATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGacac                                                                 >  2:209087/1‑149 (MQ=255)
                                                                                    aTTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCa                                                               <  7:153304/148‑1 (MQ=255)
                                                                                    aTTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCAt                                                              >  5:206312/1‑149 (MQ=255)
                                                                                     ttttGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGa                                                                              >  5:14371/1‑132 (MQ=255)
                                                                                     ttttGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGa                                                                              <  6:14371/132‑1 (MQ=255)
                                                                                             gCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCtgg                                                      <  7:276201/148‑1 (MQ=255)
                                                                                                 tGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGtc                                                 <  6:206312/149‑1 (MQ=255)
                                                                                                    tGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGtctc                                               >  8:196877/1‑148 (MQ=255)
                                                                                                      gTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCAc                                            >  3:223206/1‑149 (MQ=255)
                                                                                                        aGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGc                                          >  7:25783/1‑149 (MQ=255)
                                                                                                          acaacaAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGCACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCtt                                        <  6:34112/149‑1 (MQ=255)
                                                                                                           caacaaAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCtt                                        <  5:60398/148‑1 (MQ=255)
                                                                                                            aacaaAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGc                                                                                      >  5:350994/1‑101 (MQ=255)
                                                                                                            aacaaAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGc                                                                                      <  6:350994/101‑1 (MQ=255)
                                                                                                               aaaaaGCAGTAATAACGACTACAGGACATCTTGAATGTGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATgg                                                                           <  8:334558/109‑1 (MQ=255)
                                                                                                               aaaaaGCAGTAATAACGACTACAGGACATCTTGAATGTGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATgg                                                                           <  8:328943/109‑1 (MQ=255)
                                                                                                               aaaaaGCAGTAATAACGACTACAGGACATCTTGAATGTGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATgg                                                                           >  7:334558/1‑109 (MQ=255)
                                                                                                               aaaaaGCAGTAATAACGACTACAGGACATCTTGAATGTGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATgg                                                                           >  7:328943/1‑109 (MQ=255)
                                                                                                                        taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAGGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCtt                          >  6:232187/1‑149 (MQ=255)
                                                                                                                          ataaCGACTACAGGACATCTTGAATATGACCATGAAATAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCttt                         <  5:376713/148‑1 (MQ=255)
                                                                                                                          ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCttt                         <  5:400400/148‑1 (MQ=255)
                                                                                                                          ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTg                        >  3:204555/1‑149 (MQ=255)
                                                                                                                          ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTg                        <  7:196877/149‑1 (MQ=255)
                                                                                                                           taaCGACTACAGGACATCTTGAAAATGAACATGAAAAAACTGGCTCCCCTGGTCTCCTCTGTTGCGCTAAGCGCCACTGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCTCTGGTTGTCACCACGCTTAACAAACCGTGCTTTGt                       <  6:192003/149‑1 (MQ=255)
                                                                                                                                 cTACCGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTTGTTGCCGCTGGTGCGCTAAGCGCCACCGTCAGTGCGAATGCGAGGGCAAAAGACACCATCGCGCTGGGGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCt                 <  4:57259/149‑1 (MQ=255)
                                                                                                                                  tACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCt                 >  1:171608/1‑148 (MQ=255)
                                                                                                                                    cAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTg                <  4:178703/147‑1 (MQ=255)
                                                                                                                                       gACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGa               >  6:82484/1‑145 (MQ=255)
                                                                                                                                         cATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATg         <  4:204555/149‑1 (MQ=255)
                                                                                                                                           tCTTGGATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc       >  5:368481/1‑149 (MQ=255)
                                                                                                                                           tCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc       >  6:139347/1‑149 (MQ=255)
                                                                                                                                           tCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc       >  8:267626/1‑149 (MQ=255)
                                                                                                                                           tCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc       >  2:252147/1‑149 (MQ=255)
                                                                                                                                            cTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcg      >  3:288976/1‑149 (MQ=255)
                                                                                                                                               gAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcgcg    >  3:124127/1‑148 (MQ=255)
                                                                                                                                                aaTATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATTCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCa  <  5:392593/149‑1 (MQ=255)
                                                                                                                                                 atatGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCa  <  8:25783/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCATTAATTCTTGGCTTCCTTAATAATGGATTGAATTTGTTAGGTGTTTCCTCCTATTACCAGATGATCGTCAAAGCGGTGGTGATTTTGCTGGCGGTGCTGGTAGACAACAAAAAGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCA  >  NC_000913/3936131‑3936423

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