Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,822 A→T 6.1% intergenic (+123/‑53) ecnB → / → sugE entericidin B membrane lipoprotein/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,8220AT 6.1% 271.7 / 7.6 82intergenic (+123/‑53)ecnB/sugEentericidin B membrane lipoprotein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (41/35);  new base T (2/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGT  >  NC_000913/4376675‑4376970
                                                                                                                                                   |                                                                                                                                                    
cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                                                                                                                                     <  4:254490/149‑1 (MQ=255)
  cagcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCt                                                                                                                                                   <  1:538111/149‑1 (MQ=255)
  cagcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCt                                                                                                                                                   >  4:199426/1‑149 (MQ=255)
       aCGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaa                                                                                                                                              <  3:141345/149‑1 (MQ=255)
       aCGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaaa                                                                                                                                              >  3:51888/1‑149 (MQ=255)
         gAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACg                                                                                                                                            >  1:296900/1‑149 (MQ=255)
           aaGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                                                                                                                          <  2:563458/149‑1 (MQ=255)
            agggcagcaATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCCGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCct                                                                                                                                         <  2:459959/146‑1 (MQ=255)
             gcgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctgat                                                                                                                                        <  1:319775/149‑7 (MQ=255)
             gcgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctgat                                                                                                                                        <  3:52360/149‑7 (MQ=255)
                   caataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCa                                                                                                                                  <  1:89525/149‑1 (MQ=255)
                    aataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAc                                                                                                                                 >  1:289525/1‑149 (MQ=255)
                         gcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                                                                                                                            >  1:311829/1‑149 (MQ=255)
                         gcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                                                                                                                            >  2:542269/1‑149 (MQ=255)
                            ataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACg                                                                                                                         >  1:152427/1‑149 (MQ=255)
                               aCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGcc                                                                                                                      <  2:567969/149‑1 (MQ=255)
                                cGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGccc                                                                                                                     <  2:183007/149‑1 (MQ=255)
                                  gTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCa                                                                                                                   >  4:64422/1‑149 (MQ=255)
                                      gACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtc                                                                                                               >  2:468649/1‑149 (MQ=255)
                                            tgtgTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGa                                                                                                         <  1:191401/149‑1 (MQ=255)
                                                  gTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGTAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                                   <  2:166282/149‑1 (MQ=255)
                                                  gTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                                   >  2:56790/1‑149 (MQ=255)
                                                  gTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                                   <  1:253734/149‑1 (MQ=255)
                                                  gTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                                   >  4:106626/1‑149 (MQ=255)
                                                   tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                                  <  3:447293/149‑1 (MQ=255)
                                                     ccGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGatat                                                                                                <  2:71531/149‑1 (MQ=255)
                                                      cGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATg                                                                                               >  2:353350/1‑149 (MQ=255)
                                                       gtttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGCCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGATCCTGATATGt                                                                                              <  2:282483/149‑1 (MQ=255)
                                                        tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                             >  2:25125/1‑149 (MQ=255)
                                                           gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTg                                                                                          <  3:145728/149‑1 (MQ=255)
                                                           gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTg                                                                                          <  2:152427/149‑1 (MQ=255)
                                                                 tCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTAt                                                                                    <  4:124735/149‑1 (MQ=255)
                                                                  cTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATc                                                                                   >  4:42375/1‑149 (MQ=255)
                                                                   tGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCt                                                                                  <  2:311829/149‑1 (MQ=255)
                                                                   tGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCt                                                                                  >  3:528042/1‑149 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGTAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCGCCGGGGACGGCCCCAATACTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                                <  4:189921/149‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                                <  3:288044/149‑1 (MQ=255)
                                                                      gATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAg                                                                               >  1:437901/1‑149 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCGATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                        >  2:188884/1‑149 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                        >  4:261803/1‑149 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                        >  2:66257/1‑149 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                        >  4:110846/1‑149 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                        >  4:211085/1‑149 (MQ=255)
                                                                              caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCt                                                                       >  1:527554/1‑149 (MQ=255)
                                                                                aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTgg                                                                     <  3:35342/149‑1 (MQ=255)
                                                                                aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGA‑TATCTTAGTTATTGCTGGt                                                                    <  3:288849/149‑1 (MQ=255)
                                                                                     aaTAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc                                                                <  4:528042/149‑1 (MQ=255)
                                                                                      aTAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgct                                                               <  1:225715/149‑1 (MQ=255)
                                                                                         gTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGa                                                            >  1:374245/1‑149 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                           >  1:251996/1‑149 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                           >  1:557615/1‑149 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                           <  3:199426/149‑1 (MQ=255)
                                                                                            gAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGt                                                         >  2:350158/1‑149 (MQ=255)
                                                                                                    ccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATggg                                                 <  3:436125/149‑1 (MQ=255)
                                                                                                                cAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGa                                     >  2:378469/1‑149 (MQ=255)
                                                                                                                 aaaaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGaa                                    <  2:487529/149‑1 (MQ=255)
                                                                                                                    atatCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGCTATGGGCCGTTGGCCTGAAata                                 <  1:353350/149‑1 (MQ=255)
                                                                                                                      atCCGTTCAACCATCAGCTTTGCAGGACGAGCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAAtata                               <  1:378469/149‑1 (MQ=255)
                                                                                                                       tCCGTTCAACCATCGGCTTTGCAGGACGAGCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATAc                              >  2:223042/1‑149 (MQ=255)
                                                                                                                       tCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATAc                              >  4:402366/1‑149 (MQ=255)
                                                                                                                       tCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATAc                              >  2:365494/1‑149 (MQ=255)
                                                                                                                           ttCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCAc                          <  3:261803/149‑1 (MQ=255)
                                                                                                                            tCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACg                         <  2:165471/149‑1 (MQ=255)
                                                                                                                              aaCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGc                       >  3:65391/1‑149 (MQ=255)
                                                                                                                               aCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCt                      >  2:403667/1‑149 (MQ=255)
                                                                                                                               aCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCt                      >  3:534992/1‑149 (MQ=255)
                                                                                                                               aCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCt                      <  3:211085/149‑1 (MQ=255)
                                                                                                                                ccATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCtt                     >  1:415791/1‑149 (MQ=255)
                                                                                                                                ccATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCtt                     <  3:296776/149‑1 (MQ=255)
                                                                                                                                ccAGCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCtt                     <  4:290344/149‑1 (MQ=255)
                                                                                                                                 cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCttt                    >  1:69889/1‑149 (MQ=255)
                                                                                                                                 cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCttt                    >  1:71045/1‑149 (MQ=255)
                                                                                                                                  aTCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTa                   <  3:206489/149‑1 (MQ=255)
                                                                                                                                  aTCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTa                   <  2:527554/149‑1 (MQ=255)
                                                                                                                                   tCAGTTTTGCAGGACGACCTGCAAACGCCTCTTATCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATCTACCCACGGCTTTAg                  <  4:409201/149‑1 (MQ=255)
                                                                                                                                     aGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTc                >  4:493585/1‑149 (MQ=255)
                                                                                                                                       gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt              >  2:512913/2‑149 (MQ=255)
                                                                                                                                       gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGATCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt              >  4:393776/2‑149 (MQ=255)
                                                                                                                                       cTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt              >  4:482158/1‑149 (MQ=255)
                                                                                                                                            cAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGAc         <  1:483524/149‑1 (MQ=255)
                                                                                                                                             aGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACg        >  2:238860/1‑149 (MQ=255)
                                                                                                                                             aGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACg        >  1:54834/1‑149 (MQ=255)
                                                                                                                                              ggACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGc       >  3:39500/1‑149 (MQ=255)
                                                                                                                                                 cgacCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTCTTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGa    <  4:463698/149‑1 (MQ=255)
                                                                                                                                                  gacCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAg   >  3:65394/1‑149 (MQ=255)
                                                                                                                                                  gacCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAg   >  3:520344/1‑149 (MQ=255)
                                                                                                                                                   acCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAgt  <  1:25125/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGT  >  NC_000913/4376675‑4376970

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