Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,641 G→T 6.2% intergenic (‑88/+20) ybcI ← / ← ybcJ DUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,6410GT 6.2% 299.4 / 8.3 105intergenic (‑88/+20)ybcI/ybcJDUF457 family inner membrane protein/ribosome‑associated protein; putative RNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (45/53);  new base T (5/2);  total (50/55)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 109 total aligned reads displayed.

TTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTGC  >  NC_000913/556492‑556790
                                                                                                                                                     |                                                                                                                                                     
ttCGCACCCAACCCTAAACCAATCCAGAGGGGAACGGCGGCGCGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTg                                                                                                                                                       <  1:62502/150‑1 (MQ=255)
ttCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTg                                                                                                                                                       <  2:304442/150‑1 (MQ=255)
 tCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCTTAGAGGGCAAAAGTATGAAAAGAGAACCGGCCTGTTGa                                                                                                                                                      >  2:1689670/1‑150 (MQ=255)
  cGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGAt                                                                                                                                                     >  2:991566/1‑150 (MQ=255)
  cGAACACAACCCTAAACCAATGCATAGGGGAAAGGCGGCTTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATTCTATAAGCATGCATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGGTGAt                                                                                                                                                     >  2:1612763/1‑150 (MQ=255)
    aaaccaacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCTATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATAc                                                                                                                                                   <  1:1275518/147‑1 (MQ=255)
      cccaacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAg                                                                                                                                                 >  2:531484/1‑150 (MQ=255)
      cccaacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAg                                                                                                                                                 >  2:436832/1‑150 (MQ=255)
        caacccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCAATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGc                                                                                                                                               >  2:33243/1‑150 (MQ=255)
         aacccTACACCCATGCATCGGGGACCGGCGGCGTGCGTAATACCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGAGGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGcc                                                                                                                                              <  2:704652/150‑1 (MQ=255)
           cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCgg                                                                                                                                            >  2:414832/1‑150 (MQ=255)
           cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTATCAACAGGCCgg                                                                                                                                            <  1:1658812/150‑1 (MQ=255)
           cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGGCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCgg                                                                                                                                            <  1:964629/150‑1 (MQ=255)
           cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGGTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGCGGGCAAAAGTCTGCAAAGAGACCCGGCCTGGTGATACAGGCCgg                                                                                                                                            >  1:321725/1‑150 (MQ=255)
           cccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGGTGGCATGATGGGCGCTTCGGCAAAATGTCGATGCAATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCgg                                                                                                                                            <  1:767205/150‑1 (MQ=255)
            ccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCggg                                                                                                                                           <  1:446204/150‑1 (MQ=255)
            ccTAAACCAAAGCATAGGGGAACGGCGGCGTGCGTAATAACGGGTGGCAGGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCAGGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCggg                                                                                                                                           <  1:333242/150‑1 (MQ=255)
             cTAATCCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTAGCAGCATGGATGAACGGGGCGGAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGa                                                                                                                                          <  1:773939/150‑1 (MQ=255)
             cTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGa                                                                                                                                          >  1:385299/1‑150 (MQ=255)
                 aCCAATGCATAGGGGAACGGCGGCGTGCGTAATATCGGTTGGCATGAGGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGATCCGGCCTGTTGATACAGGCCGGGAAAgg                                                                                                                                      <  2:1345067/150‑1 (MQ=255)
                     aTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATc                                                                                                                                  <  2:610547/150‑1 (MQ=255)
                      tGCATCGGGGACCGGCGGCGTGCGTAATATCGGTTGGCATGAAGGCCGCTTCGGCAAAATGTCGATCCTATCAGCATGGATGAACGGTGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCa                                                                                                                                 <  2:1577169/150‑1 (MQ=255)
                       tcATAGTGGCCCGGCGGCGTGAGTAATAACGGTTGGCATGATGGTCGCTTCCGCAAAATCTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAACGTCTGCACAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAg                                                                                                                                <  1:163767/149‑1 (MQ=255)
                       gCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAg                                                                                                                                <  2:373436/150‑1 (MQ=255)
                        cATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAgg                                                                                                                               >  2:439086/1‑150 (MQ=255)
                        cATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCGATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAgg                                                                                                                               <  1:970823/150‑1 (MQ=255)
                          tAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGca                                                                                                                             <  1:1273716/150‑1 (MQ=255)
                          tAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGca                                                                                                                             <  2:1563131/150‑1 (MQ=255)
                          tAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGca                                                                                                                             >  2:1205720/1‑150 (MQ=255)
                           aGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGcaa                                                                                                                            >  2:1748035/1‑150 (MQ=255)
                                     cggcgTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTAc                                                                                                                  >  1:469340/1‑150 (MQ=255)
                                      ggcgTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACg                                                                                                                 <  2:1525419/150‑1 (MQ=255)
                                       gcgtgcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGGATGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGc                                                                                                                <  2:183918/150‑1 (MQ=255)
                                        cgtgcgtAATACCGGATGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCt                                                                                                               <  2:307489/150‑1 (MQ=255)
                                           gcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACACCCTGTACGCtgtg                                                                                                            >  2:206334/1‑150 (MQ=255)
                                           gcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                                                                            <  2:712040/150‑1 (MQ=255)
                                           gcgtAATAACGGTCGGCATGATGGGCGCTTCGGCAAAATGTCGATGCTAGCAGCATGGATGAACGGGTCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCtgtg                                                                                                            <  2:928908/150‑1 (MQ=255)
                                              taataaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCCGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGAcc                                                                                                         >  2:121667/1‑150 (MQ=255)
                                                 taaCGGTTGGCATGATGGGCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACAGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGc                                                                                                      <  1:1788869/150‑1 (MQ=255)
                                                  aaCGTTTGGCATTAGGGTCGCTTCGGCAAAATGTCGCTGCTATCAGCATGGAGGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACCGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCa                                                                                                     <  2:303412/150‑1 (MQ=255)
                                                  aaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCa                                                                                                     <  2:856970/150‑1 (MQ=255)
                                                   aCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCaa                                                                                                    >  2:1253351/1‑150 (MQ=255)
                                                        tGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACta                                                                                               >  2:969946/1‑149 (MQ=255)
                                                        tGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTc                                                                                               <  1:1193091/150‑1 (MQ=255)
                                                          gCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCaa                                                                                             >  2:1753061/1‑149 (MQ=255)
                                                          gCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCAc                                                                                             >  2:1416842/1‑150 (MQ=255)
                                                          gCATGATGGTCGCTTCGGCAAAAGGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCTGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACACCCTGTACGCTGTGACCTGCAAAACTCAc                                                                                             >  2:245811/1‑150 (MQ=255)
                                                                tgGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCtn                                                                                       <  2:449036/150‑2 (MQ=255)
                                                                     gCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCgg                                                                                  >  1:947718/1‑150 (MQ=255)
                                                                         cGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgac                                                                              >  1:1031228/1‑150 (MQ=255)
                                                                             aaaaTGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGAAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGCCCGGCGACGAtt                                                                          >  1:333573/1‑150 (MQ=255)
                                                                                 cgTCGATGCTATCAGCCTGGATGCCCGGGGCGTAGAGGGCAAAAGTCTCACAAGAGAACCGCCCTGTTGATACAGGCCGGGAAAGGGATCAGGCCACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGc                                                                      <  2:1518189/149‑1 (MQ=255)
                                                                                   tCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAg                                                                    <  2:679885/150‑1 (MQ=255)
                                                                                        gCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGCCCTGCAAAACTCACTGTCAGACCGGCGCCGATTTTGCAGCGttt                                                               >  2:1155399/1‑150 (MQ=255)
                                                                                        gCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGttt                                                               >  1:588132/1‑150 (MQ=255)
                                                                                          tATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCGGTACGCTGTGACATGCAAAACGCAATGTCTGACCGGAGACGATTTTGCAGCGttttc                                                             >  2:1090525/1‑148 (MQ=255)
                                                                                          tATCAGCAGTGCGGAACGGGGCGTAGAGGGCAAAAGTCGGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTgc                                                             <  1:1509877/150‑1 (MQ=255)
                                                                                            tCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTgcgc                                                           >  1:682622/1‑150 (MQ=255)
                                                                                               gCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGtt                                                        >  2:619604/1‑150 (MQ=255)
                                                                                               gCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGtt                                                        >  1:1757041/1‑150 (MQ=255)
                                                                                                cATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGttt                                                       <  2:385299/150‑1 (MQ=255)
                                                                                                cATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGttt                                                       <  1:616206/150‑1 (MQ=255)
                                                                                                cATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGttt                                                       >  2:1251288/1‑150 (MQ=255)
                                                                                                 aTGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTc                                                      <  2:1740297/150‑1 (MQ=255)
                                                                                                  tGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGATAACCGCCCTGTTGATACAGTCCGGGAAAGGGATCAGGCAACAACATGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCa                                                     <  2:861258/150‑1 (MQ=255)
                                                                                                   ggATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCaa                                                    >  2:485262/1‑150 (MQ=255)
                                                                                                    gATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAc                                                   >  2:1594013/1‑150 (MQ=255)
                                                                                                     aTGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAAcc                                                  <  1:1748035/150‑1 (MQ=255)
                                                                                                            gggCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGt                                           >  2:1005258/1‑150 (MQ=255)
                                                                                                                gTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGAc                                       <  1:1253351/150‑1 (MQ=255)
                                                                                                                  cgagGGCACAAGTCTGACAAGAGCACCGGCCTGTGGATACAGGCCGGGAAAGGGCTCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGGCTGACCGGCGACGATGTTGCAGCGGGTGCGCGTTTCAACCGCACCGTCGACtt                                     <  1:859723/149‑1 (MQ=255)
                                                                                                                   gcgggCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGTGATCCGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACttt                                    <  2:1380097/148‑1 (MQ=255)
                                                                                                                   gagGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACAGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACAGTCGACttt                                    >  2:63406/1‑150 (MQ=255)
                                                                                                                    agGGCAAAAGTCAGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCCCTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCATCCGCACCGGCGACTTTc                                   >  1:311619/1‑150 (MQ=255)
                                                                                                                      ggCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCAc                                 <  1:991566/150‑1 (MQ=255)
                                                                                                                       gCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAACCTCACTGTCGGCCCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTcccc                                >  2:1503777/1‑150 (MQ=255)
                                                                                                                        cAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGCCCGGCGACGATTTTGCTGCGTTTGCGCGGTTCCAGCGCCCCGTCGGCTTTCACCt                               >  1:728749/1‑150 (MQ=255)
                                                                                                                        cAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGCCCGGCGACGATTTTGCAGCGTTTGCGCGTTTCATCCGCACCGTCGACTTTCACCt                               >  1:61409/1‑150 (MQ=255)
                                                                                                                               cTGATAAGAGATCCGGCCTGTTTATACAGGACGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTc                        <  2:1213792/150‑1 (MQ=255)
                                                                                                                               cTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTc                        <  2:1791655/150‑1 (MQ=255)
                                                                                                                                tGAAAAGAGAACCGGGCTGTTGCTACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                       <  1:43020/150‑1 (MQ=255)
                                                                                                                                tGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGAGTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCg                       >  2:1547256/1‑150 (MQ=255)
                                                                                                                                 gAAAAGAGAACCGGCCTGTTGATTCAGGCCGGGAAAGGGATCAGGCCACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTGGGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCgg                      <  1:273032/150‑1 (MQ=255)
                                                                                                                                 gAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCgg                      <  2:82566/150‑1 (MQ=255)
                                                                                                                                 gAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCgg                      >  2:1557366/1‑150 (MQ=255)
                                                                                                                                  aaaaGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGc                     <  1:1205720/150‑1 (MQ=255)
                                                                                                                                   aaaGTGAACCGGCCTGTTGATACAGGCCGGGAAAGTGATCAGGCAACAACGTGTACGCTGTGACCTGCAAAACTCACTGTCTGAGCGGCGACGATTTTGCAGCGTTTGCGCGTTTCATCCGCACCGTCGACTTTCACCTGGCCTTCGGCa                    <  2:128041/150‑1 (MQ=255)
                                                                                                                                     agagAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAACCTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCCCCGTCGACTTTCACCTGGCCTTCGGCact                  >  1:663803/1‑148 (MQ=255)
                                                                                                                                         tcccGGCCTGTATCAACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcgcg              >  2:1049929/3‑150 (MQ=255)
                                                                                                                                         aaCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATcgcg              <  2:669702/150‑1 (MQ=255)
                                                                                                                                          cccGGCCTGTATCAACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGa             >  1:566438/2‑150 (MQ=255)
                                                                                                                                          cccGGCCTGTATCAACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGa             >  2:1423441/2‑150 (MQ=255)
                                                                                                                                           ccGGCCTGTATCAACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACAGTCGACTTTCACCTGGCCTTCGGCAATCGCGAt            >  2:26419/1‑150 (MQ=255)
                                                                                                                                               ccTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTc        <  1:1663812/150‑1 (MQ=255)
                                                                                                                                                cTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCg       <  1:1155399/150‑1 (MQ=255)
                                                                                                                                                  gTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGcc     <  1:619604/150‑1 (MQ=255)
                                                                                                                                                    tGATACAGGCCGGGAAAGGGATCAGGCAACATCCTGTACGCTGTGACCTTCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTg   <  1:33243/150‑1 (MQ=255)
                                                                                                                                                    tGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGCCGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTg   >  1:241716/1‑150 (MQ=255)
                                                                                                                                                     gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTgc  >  2:568273/1‑150 (MQ=255)
                                                                                                                                                     gATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTgc  >  1:1317679/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGTCGACTTTCACCTGGCCTTCGGCAATCGCGATTTTCGCCTGC  >  NC_000913/556492‑556790

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