Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,896 A→G V55A (GTG→GCG)  rhtB ← homoserine, homoserine lactone and S‑methyl‑methionine efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,8960AG97.9% 349.3 / ‑1.5 96V55A (GTG→GCG) rhtBhomoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (48/46);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GAGGCCAGCGATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAA  >  NC_000913/4008749‑4009040
                                                                                                                                                   |                                                                                                                                                
gagGCCAGCGATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGcgca                                                                                                                                                >  1:39894/1‑150 (MQ=255)
    ccAGCGATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATg                                                                                                                                           >  1:695281/1‑151 (MQ=255)
     cAGCGATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATg                                                                                                                                           >  1:1624666/1‑150 (MQ=255)
        cGATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATc                                                                                                                                       >  2:869883/1‑151 (MQ=255)
         gATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCg                                                                                                                                      >  1:357567/1‑151 (MQ=255)
         gATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGAAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATc                                                                                                                                       >  1:1593072/1‑150 (MQ=255)
                    aaTTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                                                                                                           >  1:1445300/1‑151 (MQ=255)
                    aaTTGCACCAGCCGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTc                                                                                                                           >  2:1663306/1‑151 (MQ=255)
                        gCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGa                                                                                                                        >  2:539314/1‑150 (MQ=255)
                         cACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGa                                                                                                                        >  1:1669150/1‑149 (MQ=255)
                            cAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGccc                                                                                                                   >  1:918312/1‑151 (MQ=255)
                            cAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGccc                                                                                                                   >  1:362479/1‑151 (MQ=255)
                               cggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                                                                                                                >  2:1577006/1‑151 (MQ=255)
                                ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                                                                                               <  2:1479614/151‑1 (MQ=255)
                                 gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                                                                                                              <  2:1489324/151‑1 (MQ=255)
                                  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                                                                                                             <  1:525692/151‑1 (MQ=255)
                                    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                                                                                                             >  2:866535/1‑149 (MQ=255)
                                     gcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                                                                                                          <  1:1087575/151‑1 (MQ=255)
                                              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                                                                                                  <  1:1183514/150‑1 (MQ=255)
                                               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGAGGGAAAATAGCGTCCCCAACCCCACGCCACCCAGCGCAATATGAATCGCCCGTCCGGTCTGAAGCCCAGCAATACACGCCatc                                                                                                   <  1:306766/148‑3 (MQ=255)
                                                aTTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                                                                                               <  1:23602/151‑1 (MQ=255)
                                                   cccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg                                                                                             <  1:743768/150‑1 (MQ=255)
                                                    ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCTAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg                                                                                           <  1:172566/151‑1 (MQ=255)
                                                     cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGa                                                                                          <  2:142814/151‑1 (MQ=255)
                                                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt                                                                                   <  2:1097891/151‑1 (MQ=255)
                                                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg                                                                                 >  2:1467732/1‑151 (MQ=255)
                                                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg                                                                                  <  2:884367/150‑1 (MQ=255)
                                                               aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGt                                                                                >  1:1573893/1‑151 (MQ=255)
                                                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGtt                                                                               <  1:1113129/151‑1 (MQ=255)
                                                                  tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTg                                                                              >  2:1643046/1‑150 (MQ=255)
                                                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGc                                                                           >  2:1624933/1‑151 (MQ=255)
                                                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGGTAACCGTGGTTGAGCg                                                                          >  1:510200/1‑150 (MQ=255)
                                                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGa                                                                         <  1:740835/151‑1 (MQ=255)
                                                                       cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCg                                                                          >  1:1074951/1‑149 (MQ=255)
                                                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAgg                                                                       >  2:427576/1‑151 (MQ=255)
                                                                             tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTc                                                                  >  2:887391/1‑151 (MQ=255)
                                                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTc                                                                  >  1:590596/1‑150 (MQ=255)
                                                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCAt                                                                >  1:1539999/1‑151 (MQ=255)
                                                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCAt                                                                <  1:1426724/151‑1 (MQ=255)
                                                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAgt                                                             <  2:1198556/151‑1 (MQ=255)
                                                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAgtgt                                                           >  1:386626/1‑151 (MQ=255)
                                                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGtt                                                          >  1:271248/1‑151 (MQ=255)
                                                                                      caacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTg                                                         >  2:1568083/1‑151 (MQ=255)
                                                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCTATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGa                                                        <  2:76966/151‑1 (MQ=255)
                                                                                        acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGAt                                                       <  2:458917/151‑1 (MQ=255)
                                                                                         cacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGAtt                                                      >  2:916265/1‑151 (MQ=255)
                                                                                          acgtcgaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTg                                                     <  2:281162/145‑1 (MQ=255)
                                                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCa                                                   <  1:1222430/150‑1 (MQ=255)
                                                                                                 aCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCaga                                              <  1:1179955/151‑1 (MQ=255)
                                                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCaga                                              <  1:1535014/150‑1 (MQ=255)
                                                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCag                                               >  2:1388768/1‑149 (MQ=255)
                                                                                                    cAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGcc                                           >  1:561471/1‑151 (MQ=255)
                                                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCt                                          <  1:1386819/151‑1 (MQ=255)
                                                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCt                                          <  1:1428823/151‑1 (MQ=255)
                                                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCt                                          <  1:1631734/151‑1 (MQ=255)
                                                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTg                                         <  1:1356498/151‑1 (MQ=255)
                                                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTg                                         <  1:1049306/151‑1 (MQ=255)
                                                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTg                                         >  2:1702437/1‑151 (MQ=255)
                                                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTg                                         >  1:388185/1‑151 (MQ=255)
                                                                                                        cACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTgg                                        <  1:269830/150‑1 (MQ=255)
                                                                                                            gAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACa                                   <  1:767878/151‑1 (MQ=255)
                                                                                                             aGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAg                                  <  1:149625/151‑1 (MQ=255)
                                                                                                               cGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGc                                >  1:79214/1‑151 (MQ=255)
                                                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCt                               <  2:749560/151‑1 (MQ=255)
                                                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCt                               >  2:1503312/1‑151 (MQ=255)
                                                                                                                 ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCtt                              >  2:359231/1‑151 (MQ=255)
                                                                                                                 ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCtt                              >  1:780417/1‑151 (MQ=255)
                                                                                                                  gAAAATAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCtt                              >  2:589131/1‑150 (MQ=255)
                                                                                                                   aaaaTAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTaa                            <  2:897580/151‑1 (MQ=255)
                                                                                                                    aaaTAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTaaa                           <  1:1203602/151‑1 (MQ=255)
                                                                                                                    aaaTAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTaa                            >  2:1400271/1‑150 (MQ=255)
                                                                                                                    aaaTAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTaa                            >  1:337257/1‑150 (MQ=255)
                                                                                                                      aTAGCGTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAAt                         <  1:477257/151‑1 (MQ=255)
                                                                                                                           gTCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCg                    >  2:221983/1‑151 (MQ=255)
                                                                                                                            tCCCCAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCg                    >  2:134076/1‑150 (MQ=255)
                                                                                                                             ccccAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGa                   <  2:1005943/150‑1 (MQ=255)
                                                                                                                               ccAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGt                <  2:1399794/151‑1 (MQ=255)
                                                                                                                                cAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTc               <  1:1332362/151‑1 (MQ=255)
                                                                                                                                cAACCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTc               >  2:175115/1‑151 (MQ=255)
                                                                                                                                 aaCCCCACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTca              <  1:1116791/151‑1 (MQ=255)
                                                                                                                                   ccccACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagc            >  2:999708/1‑151 (MQ=255)
                                                                                                                                    cccACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagca           <  2:1330247/151‑1 (MQ=255)
                                                                                                                                    cccACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagc            >  1:850342/1‑150 (MQ=255)
                                                                                                                                    cccACGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagc            >  1:853254/1‑150 (MQ=255)
                                                                                                                                       aCGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGt        >  2:903383/1‑151 (MQ=255)
                                                                                                                                        cGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGt        <  1:60038/150‑1 (MQ=255)
                                                                                                                                        cGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGt        >  2:1041399/1‑150 (MQ=255)
                                                                                                                                        cGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTa       <  1:600327/151‑1 (MQ=255)
                                                                                                                                        cGCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTa       <  1:1472584/151‑1 (MQ=255)
                                                                                                                                         gCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAg      <  1:1377659/151‑1 (MQ=255)
                                                                                                                                         gCCAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAg      <  1:1345384/151‑1 (MQ=255)
                                                                                                                                          ccAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAgg     >  1:1602119/1‑151 (MQ=255)
                                                                                                                                          ccAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAgg     >  1:1175769/1‑151 (MQ=255)
                                                                                                                                          ccAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAgg     >  1:106725/1‑151 (MQ=255)
                                                                                                                                           cAACCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGc    <  1:382647/151‑1 (MQ=255)
                                                                                                                                             aCCAGCGCAATATGAATCGCCAGTTCGGTCTGAAGCCCAGCAATAGACGCCCCCGCGCCGCGATAACCGTGGTTGAGCGAGATGGTCATAGTGTTGAATGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGCTAGGCaa  <  1:427541/151‑1 (MQ=255)
                                                                                                                                             aCCAGCGCAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCaa  <  2:304705/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GAGGCCAGCGATTTAAGGTCAATTGCACCAGCGGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAA  >  NC_000913/4008749‑4009040

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