Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,373,280 T→A 5.1% intergenic (+118/+85) yjeI → / ← yjeJ DUF4156 family lipoprotein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,373,2800TA 5.1% 336.6 / 5.8 97intergenic (+118/+85)yjeI/yjeJDUF4156 family lipoprotein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (47/45);  new base A (2/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCG  >  NC_000913/4373137‑4373424
                                                                                                                                               |                                                                                                                                                
gTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                             >  3:585446/1‑149 (MQ=255)
gTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                             <  3:283237/149‑1 (MQ=255)
   aGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                          >  3:279213/1‑149 (MQ=255)
      tttATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                       >  3:317052/1‑149 (MQ=255)
        tataAGTGACCGAACTGATGTCATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGGTTCGCTTATCAGGCCTACATAAGGCATAGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                     <  3:223603/149‑1 (MQ=255)
         ataAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGGCGCATCCgg                                                                                                                                    >  2:80134/1‑149 (MQ=255)
          taAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                   <  4:43760/149‑1 (MQ=255)
           aaGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                  <  2:647433/149‑2 (MQ=38)
            aGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                                 >  3:306559/1‑149 (MQ=255)
             gTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                                                                                                                <  4:415836/149‑1 (MQ=255)
             gTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                                                                                                                >  4:111591/1‑149 (MQ=255)
               gCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAc                                                                                                                              >  1:764520/1‑149 (MQ=255)
                       tGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATc                                                                                                                      >  4:363151/1‑149 (MQ=255)
                       tGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATc                                                                                                                      >  2:730493/1‑149 (MQ=255)
                         aTGTGATCCAGAGGCTGTGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                    <  4:474645/149‑1 (MQ=255)
                         aTGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                    <  2:513115/149‑1 (MQ=255)
                             gATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                   <  3:551647/130‑1 (MQ=38)
                             gATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                   >  4:551647/1‑130 (MQ=38)
                              aTCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTAt                                                                                                               >  2:697424/1‑149 (MQ=255)
                               tCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATc                                                                                                              >  2:673956/1‑149 (MQ=255)
                               tCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATc                                                                                                              <  1:80134/149‑1 (MQ=255)
                                 cAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTAt                                                                                                               >  3:82954/1‑146 (MQ=255)
                                 cAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTAt                                                                                                               <  4:82954/146‑1 (MQ=255)
                                          ggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAg                                                                                                   <  4:307649/149‑1 (MQ=255)
                                           gTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAgg                                                                                                  >  2:543244/1‑149 (MQ=255)
                                                aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGa                                                                                             <  1:299074/149‑1 (MQ=255)
                                                 aaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGAt                                                                                            <  1:697424/149‑1 (MQ=255)
                                                  aCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATg                                                                                           >  1:520069/1‑149 (MQ=255)
                                                    gTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGc                                                                                         >  3:36824/1‑149 (MQ=255)
                                                       taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTcc                                                                                      >  4:247850/1‑149 (MQ=255)
                                                       taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTcc                                                                                      >  2:487324/1‑149 (MQ=255)
                                                         aTTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAg                                                                                    >  3:698364/1‑149 (MQ=255)
                                                         aTTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAg                                                                                    <  2:616252/149‑1 (MQ=255)
                                                             ccTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTc                                                                                <  4:36824/149‑1 (MQ=255)
                                                             ccTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTc                                                                                <  3:20701/149‑1 (MQ=255)
                                                             ccTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTc                                                                                <  2:112709/149‑1 (MQ=255)
                                                               tGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGc                                                                              <  1:601527/149‑1 (MQ=255)
                                                                 aTGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCAt                                                                            >  1:397790/1‑149 (MQ=255)
                                                                   gcgcTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTa                                                                          <  3:545743/149‑1 (MQ=255)
                                                                    cgcttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAg                                                                         <  1:543244/149‑1 (MQ=255)
                                                                     gcttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGc                                                                        >  4:536453/1‑149 (MQ=255)
                                                                     gcttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGc                                                                        >  4:505052/1‑149 (MQ=255)
                                                                           cttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGAt                                                                  >  4:411573/1‑149 (MQ=255)
                                                                            ttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATg                                                                 >  3:522259/1‑149 (MQ=255)
                                                                             tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATga                                                                <  4:448974/149‑1 (MQ=255)
                                                                             tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATga                                                                <  4:221064/149‑1 (MQ=255)
                                                                             tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATga                                                                <  2:667532/149‑1 (MQ=255)
                                                                               tCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgaga                                                              <  4:317052/149‑1 (MQ=255)
                                                                                 aGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTGGCATTAGCAAAGATGAGatt                                                            >  3:622415/1‑149 (MQ=255)
                                                                                  ggCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGatta                                                           <  1:320844/149‑1 (MQ=255)
                                                                                   gCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGattat                                                          >  2:356105/1‑149 (MQ=255)
                                                                                   gCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGattat                                                          >  2:503925/1‑149 (MQ=255)
                                                                                   gCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGattat                                                          >  1:563717/1‑149 (MQ=255)
                                                                                     cTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTa                                                                          <  1:708355/131‑1 (MQ=255)
                                                                                     cTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTa                                                                          >  2:708355/1‑131 (MQ=255)
                                                                                     cTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTAttt                                                        >  2:698774/1‑149 (MQ=255)
                                                                                        cATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGc                                                     <  3:247850/149‑1 (MQ=255)
                                                                                          tAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCt                                                   <  4:698364/149‑1 (MQ=255)
                                                                                          tAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCt                                                   >  3:71085/1‑149 (MQ=255)
                                                                                            aGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTa                                                                          <  1:292511/124‑1 (MQ=255)
                                                                                            aGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTa                                                                          >  2:292511/1‑124 (MQ=255)
                                                                                            aGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                 >  1:343501/1‑149 (MQ=255)
                                                                                            aGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                 <  4:306559/149‑1 (MQ=255)
                                                                                             ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGc                                                     >  1:705434/1‑144 (MQ=255)
                                                                                             ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGc                                                     <  2:705434/144‑1 (MQ=255)
                                                                                             ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgtg                                                >  3:200693/1‑149 (MQ=255)
                                                                                             ggCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgtg                                                <  1:743993/149‑1 (MQ=255)
                                                                                               cATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTgtggt                                              <  2:640730/149‑1 (MQ=255)
                                                                                                aTTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTgtggtg                                             >  3:217990/1‑149 (MQ=255)
                                                                                                  tGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCa                                           <  3:445977/149‑1 (MQ=255)
                                                                                                       cTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGtttt                                      <  4:267507/149‑1 (MQ=255)
                                                                                                       cTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGtttt                                      <  4:357216/149‑1 (MQ=255)
                                                                                                         gATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGc                                    <  3:536453/149‑1 (MQ=255)
                                                                                                           ttGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtg                                  >  2:448242/1‑149 (MQ=255)
                                                                                                           ttGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtg                                  >  1:341890/1‑149 (MQ=255)
                                                                                                              aaTGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                               >  4:745253/1‑149 (MQ=255)
                                                                                                              aaTGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                               <  4:50093/149‑1 (MQ=255)
                                                                                                                  cTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                           <  1:254317/149‑1 (MQ=255)
                                                                                                                  cTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                           <  4:585446/149‑1 (MQ=255)
                                                                                                                     gagaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                        >  4:719654/1‑149 (MQ=255)
                                                                                                                     gagaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                        >  2:765849/1‑149 (MQ=255)
                                                                                                                     gagaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                        <  1:503925/149‑1 (MQ=255)
                                                                                                                      agaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAac                       >  4:22573/1‑149 (MQ=255)
                                                                                                                        aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                 >  3:574307/1‑121 (MQ=255)
                                                                                                                        aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                 <  4:574307/121‑1 (MQ=255)
                                                                                                                          tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGt                   >  2:662279/1‑149 (MQ=255)
                                                                                                                          tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGt                   <  1:487324/149‑1 (MQ=255)
                                                                                                                            tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTaa                 <  2:397790/149‑1 (MQ=255)
                                                                                                                              tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGt               <  2:520069/149‑1 (MQ=255)
                                                                                                                               atgCCGGATGCGGCGTAAACGCCTTATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTg              >  2:43407/3‑149 (MQ=255)
                                                                                                                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGAt            <  3:548187/149‑1 (MQ=255)
                                                                                                                                 gCCGGATAAGGCGTTTACGCCGCATCCGGCATAACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACGCGTAAGTGAt            >  4:649431/1‑149 (MQ=255)
                                                                                                                                  ccGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATg           <  4:252136/149‑1 (MQ=255)
                                                                                                                                        aaGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgcgc     >  3:459476/1‑149 (MQ=255)
                                                                                                                                          ggCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgcgcgc   <  1:673956/149‑1 (MQ=255)
                                                                                                                                           gCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGcgcgcg  >  4:147929/1‑149 (MQ=255)
                                                                                                                                           gCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGcgcgcg  >  1:407550/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCG  >  NC_000913/4373137‑4373424

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