Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 7.9% intergenic (+131/+14) rraB → / ← yjgM ribonuclease E inhibitor protein B/putative acetyltransferase YjgM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,0200AT 7.9% 252.8 / 8.5 88intergenic (+131/+14)rraB/yjgMribonuclease E inhibitor protein B/putative acetyltransferase YjgM
Reads supporting (aligned to +/- strand):  ref base A (44/36);  new base T (2/5);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGC  >  NC_000913/4478871‑4479167
                                                                                                                                                     |                                                                                                                                                   
tGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                                                                     <  1:1250987/150‑1 (MQ=255)
 gACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAacac                                                                                                                                                    <  1:1056716/150‑1 (MQ=255)
    ggAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTcgcg                                                                                                                                                 <  2:136808/150‑1 (MQ=255)
    ggAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                                                                                 <  2:359119/150‑1 (MQ=255)
      aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTTATGTCAAATCAGGCATAAACAAAAAGTACGTTTTTAATGTTCTAATGCCCGCGTGTTCAACTCGCggg                                                                                                                                               <  2:848123/150‑1 (MQ=255)
      aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                                                                               >  1:91832/1‑150 (MQ=255)
       gTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGc                                                                                                                                              <  1:551299/150‑1 (MQ=255)
         tCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAt                                                                                                                                            <  2:721622/150‑1 (MQ=255)
             cACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttcga                                                                                                                                        <  1:640334/150‑4 (MQ=255)
             cACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttt                                                                                                                                        >  1:157117/1‑150 (MQ=255)
              aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                                                                       >  1:709938/1‑150 (MQ=255)
              aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                                                                       <  2:275028/150‑1 (MQ=255)
               cTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                                                                                                                      <  2:1382296/150‑1 (MQ=255)
                 aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTag                                                                                                                                    >  2:1025953/1‑150 (MQ=255)
                  aTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                                                                                                                   >  2:1261771/1‑150 (MQ=255)
                    tAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATTAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                                                                                                                                 >  1:937571/1‑150 (MQ=255)
                     aaGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                                                                >  2:740998/1‑150 (MQ=255)
                     aaGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                                                                >  2:830691/1‑150 (MQ=255)
                       gCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcg                                                                                                                              <  1:1069255/150‑1 (MQ=255)
                       gCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcg                                                                                                                              >  1:1272797/1‑150 (MQ=255)
                           aGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAATCATGTCACATCCGGCATAAACAAAACGTACGTTGTAAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCa                                                                                                                          <  2:1413569/150‑1 (MQ=255)
                           aGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAGTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGGTGTTAATGTTCGAATGCCCGCGAGGTGAACACGCGGGCATTTTTTTAGAGTTCGCGCa                                                                                                                          <  2:1407911/150‑1 (MQ=255)
                            gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCaa                                                                                                                         <  2:727026/150‑1 (MQ=255)
                             cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAAc                                                                                                                        >  1:1524120/1‑150 (MQ=255)
                              gcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                                                                       <  1:219041/150‑1 (MQ=255)
                              gcATCAGGCCAAGGTAGTTCCTCCAACCTCGTCAATTTGTAGACCGGATAATGCAATCATCACACAACCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                                                                       >  2:423230/1‑150 (MQ=255)
                               cATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                                                                      >  1:396191/1‑150 (MQ=255)
                                    ggCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                                                                 >  2:286058/1‑150 (MQ=255)
                                        cAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                                                                             >  1:1386874/1‑150 (MQ=255)
                                          ggTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                                                                           <  1:1213170/150‑1 (MQ=255)
                                          ggTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac                                                                                                           <  1:440564/150‑1 (MQ=255)
                                                tCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGa                                                                                                     >  1:735446/1‑150 (MQ=255)
                                                  cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                                                                   <  2:91832/150‑1 (MQ=255)
                                                   tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                                                                                                                      <  2:1542279/114‑1 (MQ=255)
                                                   tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                                                                                                                      >  1:1542279/1‑114 (MQ=255)
                                                           ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                          <  2:157117/150‑1 (MQ=255)
                                                           ggTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                          <  2:903451/150‑1 (MQ=255)
                                                              gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCac                                                                                       >  2:1385163/1‑149 (MQ=255)
                                                              gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                                                                                       >  2:1379423/1‑150 (MQ=255)
                                                              gAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTAGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAg                                                                                       >  1:517167/1‑150 (MQ=255)
                                                                 tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCa                                                                                    <  2:707560/150‑1 (MQ=255)
                                                                    gTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                                                 >  2:895271/1‑150 (MQ=255)
                                                                     tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                                                                <  1:740998/150‑1 (MQ=255)
                                                                       gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGCAATTTTTTTAGAGTTAGCGCAACATCCGCACTTCACAAAAGAAATGGCCCGTGAAGCAAAGCGCa                                                                              >  1:23232/1‑150 (MQ=255)
                                                                        cccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                                             <  2:1186222/150‑1 (MQ=255)
                                                                            gATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                                         >  2:1542933/1‑150 (MQ=255)
                                                                            gATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                                         <  1:895271/150‑1 (MQ=255)
                                                                             aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                                                                        >  1:957464/1‑150 (MQ=255)
                                                                                 ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                                                    >  2:1094744/1‑150 (MQ=255)
                                                                                   cATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                                                  >  1:1320180/1‑150 (MQ=255)
                                                                                    aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                                                                 <  2:1579724/150‑1 (MQ=255)
                                                                                      tCATGTCACATCCGGCATAAACAAATCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACATTCGACTTGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCa                                                               <  1:906962/150‑1 (MQ=255)
                                                                                      tCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCa                                                               <  2:709938/150‑1 (MQ=255)
                                                                                       cATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa                                                              <  1:1261771/150‑1 (MQ=255)
                                                                                       cATGTCACATCCGGCATAAACAAAACGTACGTTGTTAAAGTTAGAATGCCCGCGAGTTGACCACGAGGGCATTTTTTTAGAGTTCGCGCAACATCCGCCCTTCACTATCGACATGGCCCGAGCAGCGCAGCGCCTAGTCGAATTGCTCaa                                                              >  2:1601967/1‑150 (MQ=255)
                                                                                         tGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                            <  2:1320180/150‑1 (MQ=255)
                                                                                         tGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                            >  2:1449994/1‑150 (MQ=255)
                                                                                             acaTCCGGCATAAACAATACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTGCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGACAAGCGCATAGTCGATATGCTCAAAGCCCa                                                        <  2:432621/150‑1 (MQ=255)
                                                                                                 ccGCCATAAACAAAACGTACGATCTTAATCTTCGAATGCCCGCGAGTTGCACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAATTCAAAATCTACATCGCCCGTGCAGCCAAGCGCATAGTCGATAAGCTCAAAGCCCAAATg                                                    >  2:1178584/1‑150 (MQ=255)
                                                                                                   ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                                                  >  1:1611186/1‑150 (MQ=255)
                                                                                                   ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                                                  >  1:1351434/1‑150 (MQ=255)
                                                                                                       tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATa                                              >  1:601328/1‑150 (MQ=255)
                                                                                                       tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATa                                              <  2:396191/150‑1 (MQ=255)
                                                                                                       tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATa                                              >  1:1322082/1‑150 (MQ=255)
                                                                                                        aaaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaa                                             >  2:1020236/1‑150 (MQ=255)
                                                                                                         aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                                            >  2:1214933/1‑150 (MQ=255)
                                                                                                            aaaaCGTACGTTATTAAAGTTCTAAAGCCCGCGAGTTAACCACGCGGGCATTTTTTTAGAAATAGAGCAACATCCGCACTTCACAATCGACATGGAAAGAGCAGCCAACCCCATAGTCGATATGCTCAAAGACCAAATGCTCATAAAgcg                                         >  2:697710/1‑150 (MQ=255)
                                                                                                             aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgcgc                                        >  1:932331/1‑150 (MQ=255)
                                                                                                                 gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg                                    <  1:423230/150‑1 (MQ=255)
                                                                                                                 gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg                                    <  1:1035617/150‑1 (MQ=255)
                                                                                                                  tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATgg                                   <  1:1025953/150‑1 (MQ=255)
                                                                                                                    cgtTGATAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCt                                 >  2:58433/1‑150 (MQ=255)
                                                                                                                      ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTc                               <  1:1094744/150‑1 (MQ=255)
                                                                                                                      ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTc                               <  2:1322082/150‑1 (MQ=255)
                                                                                                                         ttAATTTTCGAATGGCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCTCGCAACATCCGCACTTCACAATAGACTTGTCCCGTGCAGCCAAGCGCTTATTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGTCTTCCtt                            <  2:1429718/150‑1 (MQ=255)
                                                                                                                         ttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCtt                            <  1:940313/150‑1 (MQ=255)
                                                                                                                          tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt                           >  2:72792/1‑150 (MQ=255)
                                                                                                                              gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGCCATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAACGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaa                       >  1:600025/1‑150 (MQ=255)
                                                                                                                              gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaa                       <  1:1170330/150‑1 (MQ=255)
                                                                                                                                 cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                    <  2:1386874/150‑1 (MQ=255)
                                                                                                                                    aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCgg                 <  1:1304629/150‑1 (MQ=255)
                                                                                                                                      gCCCGCGAGTTGAACCCGCGGGCAGTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGCCATGGCCCGTGCAGCCAAGCGCATAGTCGAAATGCTCAAAGACCAAATGCTCATAAAGCGCAATGGGTTCCTTTAAAAAAGCGGTc               >  2:1617928/1‑150 (MQ=255)
                                                                                                                                      gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCCCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTc               >  1:1246070/1‑150 (MQ=255)
                                                                                                                                       cccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg              >  1:138305/1‑150 (MQ=255)
                                                                                                                                       cccGCGTGTTCAACTCGCGGGCATTCTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg              >  1:1564140/1‑150 (MQ=255)
                                                                                                                                       cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg              >  1:1374602/1‑150 (MQ=255)
                                                                                                                                       cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg              >  2:483497/1‑150 (MQ=255)
                                                                                                                                         cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGtt            >  2:1034006/1‑150 (MQ=255)
                                                                                                                                             aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCa        <  1:570927/150‑1 (MQ=255)
                                                                                                                                              gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAg       >  1:230111/1‑150 (MQ=255)
                                                                                                                                                   acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATagc  <  1:286058/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGC  >  NC_000913/4478871‑4479167

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