Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 975,700 G→T 5.5% D27Y (GAC→TAC) ‡ mukE → chromosome condensin MukBEF, MukE localization factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913975,7000GT 5.5% 363.7 / 7.8 118D27Y (GAC→TAC) ‡mukEchromosome condensin MukBEF, MukE localization factor
Reads supporting (aligned to +/- strand):  ref base G (49/62);  new base T (4/3);  total (53/65)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 121 total aligned reads displayed.

GCAAGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTACG  >  NC_000913/975551‑975848
                                                                                                                                                     |                                                                                                                                                    
gcAAGCAGATTTCACCTGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTgg                                                                                                                                                      <  2:1426037/150‑1 (MQ=255)
gcAAGCAGAATTCACCGGACTGCCAGCGAAATGGCAGCCGATTAAGGATTACTTTGCAAATGTACAGGCGCATGTCATCGACAAATATTGAACATGTGATGCCGGTTAATCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTgg                                                                                                                                                      <  2:1225311/150‑1 (MQ=255)
  aaGCAGAATTCACCCGACTGCCCGCGAAATGGCAGCCGATTAATGATTACGGAGCCACGGTACAGGCGCATGACATCGACAAATATTGAACAAGTGATGCCGGTTAAGCGGGCGCAGGCGCTGGCGAATCCGTTATTACCGGCGCTGGAc                                                                                                                                                    <  2:78731/150‑1 (MQ=255)
   aGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACa                                                                                                                                                   >  2:193532/1‑150 (MQ=255)
   aGCAGATTTCACCGGACTGCCAGCGAAATAGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTACGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGAcc                                                                                                                                                   >  1:1598390/1‑149 (MQ=255)
     cAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgc                                                                                                                                                 >  2:1401159/1‑150 (MQ=255)
     cAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgc                                                                                                                                                 <  1:732549/150‑1 (MQ=255)
         tttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCt                                                                                                                                             >  1:1924556/1‑150 (MQ=255)
          ttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAGCGccgg                                                                                                                                            <  2:747495/150‑3 (MQ=255)
          ttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAGCGccgg                                                                                                                                            <  2:442271/150‑3 (MQ=255)
          ttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAACGccgg                                                                                                                                            <  1:1560988/150‑3 (MQ=255)
           tCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGGCTAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAAGCCATTATTTCCGGCGCTGGACAGCGCCtgg                                                                                                                                           >  1:342774/1‑148 (MQ=255)
            cACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                                                                                                                                          >  1:143492/1‑150 (MQ=255)
                 gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCa                                                                                                                                     <  1:685964/150‑1 (MQ=255)
                 gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCa                                                                                                                                     <  1:1389598/150‑1 (MQ=255)
                    tGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGa                                                                                                                                  >  2:1384437/1‑150 (MQ=255)
                     gCCAGCGAAATGGCAGCCGATTAATGATTAGGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGAc                                                                                                                                 <  1:65147/150‑1 (MQ=255)
                     gCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGAc                                                                                                                                 <  1:615492/150‑1 (MQ=255)
                      ccAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                                                                                                                                >  2:347200/1‑150 (MQ=255)
                      ccAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                                                                                                                                >  2:1480364/1‑150 (MQ=255)
                       cAGCGAAATGGCAGCCGATTAATGATTCCGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                               <  2:159033/150‑1 (MQ=255)
                       cAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                               >  1:642760/1‑150 (MQ=255)
                        aGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCTTTACGTTCAGGACGcc                                                                                                                              <  2:100686/150‑1 (MQ=255)
                         gCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACACGTGATGCCGGTTAAGCGGGCGCCGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCGTACGTTCAGGACGCCa                                                                                                                             <  2:854946/150‑1 (MQ=255)
                         gCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCa                                                                                                                             <  2:854943/150‑1 (MQ=255)
                         gCGAAATGGCAGCAGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCAGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCa                                                                                                                             >  2:1804127/1‑150 (MQ=255)
                              aTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTg                                                                                                                        <  2:375152/150‑1 (MQ=255)
                                   aGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAATTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGCCGCCATATTTGCCTc                                                                                                                   >  1:355515/1‑150 (MQ=255)
                                    gCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcg                                                                                                                  >  2:1476424/1‑150 (MQ=255)
                                       gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGACCAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               <  2:1230143/150‑1 (MQ=255)
                                       gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               <  2:30773/150‑1 (MQ=255)
                                        aTTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCGGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacga                                                                                                              <  2:1248995/150‑1 (MQ=255)
                                           aaTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACt                                                                                                           <  1:525315/150‑1 (MQ=255)
                                           aaTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACt                                                                                                           >  2:897071/1‑150 (MQ=255)
                                            aTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGTACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTg                                                                                                          <  1:94497/150‑1 (MQ=255)
                                            aTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTg                                                                                                          >  1:1878556/1‑150 (MQ=255)
                                             tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCTCAGGCGCTGGCGAATCCGTTATTTCCGTCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGACTCGACGAACTgg                                                                                                         >  1:1237855/1‑150 (MQ=255)
                                             tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGAGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTgg                                                                                                         >  1:1971154/1‑150 (MQ=255)
                                               aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGAt                                                                                                       >  2:45554/1‑150 (MQ=255)
                                               aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTAAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGAt                                                                                                       <  2:1094257/150‑1 (MQ=255)
                                                ttACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATa                                                                                                      >  1:927601/1‑150 (MQ=255)
                                                  aCGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAAt                                                                                                    >  1:1525124/1‑150 (MQ=255)
                                                  aCGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAAt                                                                                                    <  2:1678442/150‑1 (MQ=255)
                                                   cGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATc                                                                                                   <  2:1742765/150‑1 (MQ=255)
                                                   cGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATc                                                                                                   <  2:753509/150‑1 (MQ=255)
                                                     gAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCAt                                                                                                 <  2:143492/150‑1 (MQ=255)
                                                      aGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATg                                                                                                <  1:1601222/150‑1 (MQ=255)
                                                         cAAGGTACAGGCGCTGGGCATCGACAAATATTGAACTAGTGAGGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCAt                                                                                             <  1:1380966/150‑1 (MQ=255)
                                                           aGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTc                                                                                           >  1:551540/1‑150 (MQ=255)
                                                             gTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCt                                                                                         <  2:279806/150‑1 (MQ=255)
                                                               aCAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGa                                                                                       >  1:1844658/1‑150 (MQ=255)
                                                                cAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGAt                                                                                      >  1:1352480/1‑150 (MQ=255)
                                                                cAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGAt                                                                                      <  1:385860/150‑1 (MQ=255)
                                                                 aGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATg                                                                                     <  2:1782431/150‑1 (MQ=255)
                                                                 aGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATg                                                                                     <  2:1932768/150‑1 (MQ=255)
                                                                 aGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATg                                                                                     <  1:1136607/150‑1 (MQ=255)
                                                                    cgcATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGGAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGAt                                                                                  <  2:759082/150‑1 (MQ=255)
                                                                              cGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAat                                                                        >  1:466343/1‑150 (MQ=255)
                                                                              cGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAat                                                                        <  1:896967/150‑1 (MQ=255)
                                                                              cGACAAATAATGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATAATGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAat                                                                        >  2:487448/1‑150 (MQ=255)
                                                                               gACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAata                                                                       <  1:59317/150‑1 (MQ=255)
                                                                                 cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     <  1:1861135/150‑1 (MQ=255)
                                                                                 cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     >  2:481829/1‑150 (MQ=255)
                                                                                 cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     >  2:901483/1‑150 (MQ=255)
                                                                                   aaTATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTg                                                                   <  1:1990207/150‑1 (MQ=255)
                                                                                    atatTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTgg                                                                  <  1:860611/150‑1 (MQ=255)
                                                                                       ttGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATATTCATGCATTCCTGATGGATTTTCAGGAATATCTGGAag                                                               >  2:195791/1‑150 (MQ=255)
                                                                                         gAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAagag                                                             >  2:726409/1‑150 (MQ=255)
                                                                                         gAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAagag                                                             >  2:868195/1‑150 (MQ=255)
                                                                                         gAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCAGATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAagag                                                             >  1:478188/1‑150 (MQ=255)
                                                                                         gAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACACCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAagag                                                             >  1:1947445/1‑150 (MQ=255)
                                                                                          aaCAAGTGCTGCCGGTTAAGCTGGCGCAGGCGCGGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGt                                                            <  1:1259721/150‑1 (MQ=255)
                                                                                               gTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAc                                                       <  2:1271017/150‑1 (MQ=255)
                                                                                                  aTGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcg                                                    <  2:909572/150‑1 (MQ=255)
                                                                                                   tGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTTCGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgc                                                   >  1:337307/1‑150 (MQ=255)
                                                                                                   tGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgc                                                   >  2:1166662/1‑150 (MQ=255)
                                                                                                    gCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgcg                                                  <  1:1112191/150‑1 (MQ=255)
                                                                                                      cGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGTATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGATGAGTTTTACGCGCGtt                                                <  1:233869/150‑1 (MQ=255)
                                                                                                           aaGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAAc                                           <  1:1716429/150‑1 (MQ=255)
                                                                                                             tcTGGCGCAGGCGCAGGCGAAACCGTTATTTCAGGCGCTGGACAGCGCCTTACTTTCACGACGCCATATTGGCCTCCACGAACTGGATAATCATGCATTCCTGATGGATATTCAGGAATATCTGGAAGAGATTTACGCGAGTAATAACGt                                         >  2:739406/2‑150 (MQ=255)
                                                                                                             gCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGt                                         <  1:347200/150‑1 (MQ=255)
                                                                                                               tGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTg                                       <  2:1774812/150‑1 (MQ=255)
                                                                                                                   gcAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCt                                   <  1:1267627/150‑1 (MQ=255)
                                                                                                                        cgcTGGCGAATCCGTTATTTCCGGCGCTTGACAGCGCCTTACGTTCAGGACGCCAAAATGGCCTCGACGAACGGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTc                              >  2:845707/1‑150 (MQ=255)
                                                                                                                            ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTAgcgc                          >  2:434753/1‑150 (MQ=255)
                                                                                                                             gCGAATCCGTTATTTCGGGCGCTGGACAGCGCCTTAGGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTACTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCa                         <  2:1237855/150‑1 (MQ=255)
                                                                                                                                aaTCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCa                      >  1:822647/1‑150 (MQ=255)
                                                                                                                                aaTCCGTTATTTCCGGCGCTGGACAGCGCATTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCa                      <  1:433600/150‑1 (MQ=255)
                                                                                                                                    cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAg                  <  2:953362/150‑1 (MQ=255)
                                                                                                                                    cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAg                  <  2:1143276/150‑1 (MQ=255)
                                                                                                                                      ttATTTCCTGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAggg                <  1:1401159/150‑1 (MQ=255)
                                                                                                                                        acgtaaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGtt              >  2:1194683/7‑150 (MQ=255)
                                                                                                                                          ttCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttct            <  1:726409/150‑1 (MQ=255)
                                                                                                                                           taaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttctt           >  1:236620/4‑150 (MQ=255)
                                                                                                                                           tCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttctt           >  2:1718789/1‑150 (MQ=255)
                                                                                                                                            ccGGCGCTGGACAGCGCCTTACGTTCAGGTCGCCATATTGGCCTCGACGAACTGGATAATCATTCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttcttc          <  1:45554/150‑1 (MQ=255)
                                                                                                                                             aggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCGTATTCGCGCACCAGAAGGGTgcttct         >  2:1213247/2‑150 (MQ=255)
                                                                                                                                              ggCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGTTTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTa        <  2:1860535/150‑1 (MQ=255)
                                                                                                                                                 gcTGGACAGCGCCTTACGATCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTAttt     <  1:730185/150‑1 (MQ=255)
                                                                                                                                                    ggACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTACg  >  2:1147103/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAAGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTACG  >  NC_000913/975551‑975848

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