Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 975,700 G→T 5.8% 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.8% 302.5 / 5.6 98D27Y (GAC→TAC) ‡mukEchromosome condensin MukBEF, MukE localization factor
Reads supporting (aligned to +/- strand):  ref base G (39/53);  new base T (4/2);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

AGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTA  >  NC_000913/975554‑975846
                                                                                                                                                  |                                                                                                                                                  
aGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACa                                                                                                                                                 >  1:1701331/1‑150 (MQ=255)
  cAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGTCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgc                                                                                                                                               >  2:324755/1‑150 (MQ=255)
  cAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAgc                                                                                                                                               <  1:162491/150‑1 (MQ=255)
  cAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgc                                                                                                                                               >  2:988473/1‑150 (MQ=255)
    gATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgcgc                                                                                                                                             <  1:159548/150‑1 (MQ=255)
     aTTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGcc                                                                                                                                            <  2:818762/150‑1 (MQ=255)
     aTTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGcc                                                                                                                                            <  1:1284354/150‑1 (MQ=255)
      tttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCt                                                                                                                                           <  2:250232/150‑1 (MQ=255)
       ttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCtt                                                                                                                                          <  1:1374973/150‑1 (MQ=255)
        tCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAGCGccgga                                                                                                                                         <  1:697544/150‑4 (MQ=255)
        tCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGAACCGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTa                                                                                                                                         <  1:695335/150‑1 (MQ=255)
         cACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGATAAGATGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                                                                                                                                        >  2:1140482/1‑150 (MQ=255)
         cACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCAGGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                                                                                                                                        <  1:1705758/150‑1 (MQ=255)
          cccGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACg                                                                                                                                       <  1:850821/149‑1 (MQ=255)
               aCTGACAGCGAAATGGCAGCAGATTAATTATTACGGAGCCAAGGGACAGGCGCAGTTCATCGACAAATATTGACCAAGTGATGCCGGTTAAGCTGGCGCAGGCTCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAg                                                                                                                                  <  1:121221/150‑1 (MQ=255)
                    cAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                             <  2:1612057/150‑1 (MQ=255)
                    cAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                             >  1:162255/1‑150 (MQ=255)
                     aGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGcc                                                                                                                            <  1:538561/150‑1 (MQ=255)
                        gAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGCCGCCata                                                                                                                         >  2:1627974/1‑150 (MQ=255)
                         aaaTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCatat                                                                                                                        >  2:1609783/1‑150 (MQ=255)
                                aGCCGATTAATGATTACGGAGCCTAGGTAGAGGCGCATGTCATCGACAATTTTTGATCAAGTGATGCGGGTGAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTc                                                                                                                 <  2:204612/150‑1 (MQ=255)
                                aGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTc                                                                                                                 >  1:1442374/1‑150 (MQ=255)
                                 tcctaTTAATGATTACGGATCCAAGGAACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCTGGACGCCATATTGGCCTcg                                                                                                                <  2:1329962/146‑1 (MQ=255)
                                 gCCGATTACTGATTACGGAGCCACGGTACAGGCGCATGTCATCGACAAATATTGAACAAGGGAGTCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcg                                                                                                                <  2:825835/150‑1 (MQ=255)
                                 gCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGCTGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcg                                                                                                                <  2:1731077/150‑1 (MQ=255)
                                 gCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcg                                                                                                                <  2:1097980/150‑1 (MQ=255)
                                 gCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcg                                                                                                                >  2:218594/1‑150 (MQ=255)
                                    gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                             <  1:373074/150‑1 (MQ=255)
                                          tGATTCCGGAGCCAAGGACCAGGCGCATGTCATCGACAAATATTGAACAAGTTATGCCGGTTAAGCGGGCGCAGGCGCTGGCGAATCCGTTATATCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTgg                                                                                                       <  1:830912/150‑1 (MQ=255)
                                          tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTgg                                                                                                       >  2:1614860/1‑150 (MQ=255)
                                           gATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGa                                                                                                      <  1:1093808/150‑1 (MQ=255)
                                            aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGAt                                                                                                     >  1:915819/1‑150 (MQ=255)
                                               aCGGAGCCAAGGTACAGGCGCATGTCATCGACAAAAATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTTGACGAACTGGATAAt                                                                                                  <  2:227410/150‑1 (MQ=255)
                                                  gAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCAt                                                                                               <  2:1664874/150‑1 (MQ=255)
                                                   aGCCAAGGTACATGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATg                                                                                              <  2:1635419/150‑1 (MQ=255)
                                                     ccAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCa                                                                                            >  2:871380/1‑150 (MQ=255)
                                                          gTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCGGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCt                                                                                       <  1:51748/150‑1 (MQ=255)
                                                            aCAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGa                                                                                     >  1:978186/1‑150 (MQ=255)
                                                            aCAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTCGATAATCATGCATTCCTGa                                                                                     >  1:1490871/1‑150 (MQ=255)
                                                            aCAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCAGATTGGCCTCGACGACCTGGATTATAATGCATTCCTGa                                                                                     >  2:1153297/1‑150 (MQ=255)
                                                                    aTGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGAtttt                                                                             <  1:1595637/150‑1 (MQ=255)
                                                                     tGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTc                                                                            <  2:1251882/150‑1 (MQ=255)
                                                                       tCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTGCGGCGCTGGACAGCGCCGTTCGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAg                                                                          <  1:589088/150‑1 (MQ=255)
                                                                       tCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAg                                                                          <  2:641231/150‑1 (MQ=255)
                                                                        cATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTACTGATGGATTTTCAgg                                                                         >  2:1550832/1‑150 (MQ=255)
                                                                          tCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGaa                                                                       <  1:598073/150‑1 (MQ=255)
                                                                            gACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAata                                                                     <  1:1451622/150‑1 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGTCTGCGCCTTACGTTCAGGACGCCATATTGGCCTGGACGAACTGGATAATCATGCATTCCTGATGGAGTTTCAGGAatata                                                                   >  1:1069476/1‑149 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                   >  2:61200/1‑150 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                   >  1:1291884/1‑150 (MQ=255)
                                                                               aaaTATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCt                                                                  >  2:391788/1‑150 (MQ=255)
                                                                                 atatTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTgg                                                                >  1:1228377/1‑150 (MQ=255)
                                                                                 atatTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTgg                                                                <  2:1672103/150‑1 (MQ=255)
                                                                                 atatTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTAAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTgg                                                                <  2:964026/150‑1 (MQ=255)
                                                                                  tatTCAACAAGTGATGCCGGTTAAGCGGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTTGa                                                               <  2:1533135/150‑1 (MQ=255)
                                                                                     tGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAagc                                                            >  2:1137884/1‑149 (MQ=255)
                                                                                     tGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAaga                                                            >  1:1821/1‑150 (MQ=255)
                                                                                     tGAACAAGTGATACCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTAATTGATATTCAGGAATATCAGGAaga                                                            >  2:886589/1‑150 (MQ=255)
                                                                                      gAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAagag                                                           <  1:1609783/150‑1 (MQ=255)
                                                                                           aGTGACGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTa                                                      <  1:665675/150‑1 (MQ=255)
                                                                                             tGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGTATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcg                                                    <  1:623992/150‑1 (MQ=255)
                                                                                               aTGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcg                                                  >  1:650418/1‑150 (MQ=255)
                                                                                                 gCCGGTTAAGCTGGCGCAGGCGCTGGAGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATCTTGGCCGCGACGACCTGGATAATCATGCATTCGTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgcg                                                >  1:99740/1‑150 (MQ=255)
                                                                                                  ccGGTTACGCGGGCGCAGGCGGTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGt                                               <  1:207199/150‑1 (MQ=255)
                                                                                                   cGGTTAAGCTGGCGCAGGCGCTGGCGAAACCGTTATTTCCGGCGCTGGCCAGCGCCTTACGGTCAGGACGCACTATTGGCGTCGACGAACTCGATAATCATGCATTCCTGATGGAGTTTCAGGAATATCTGGAAGAGTTTTTCGCGCCtt                                              >  1:126044/1‑150 (MQ=255)
                                                                                                     gTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTtat                                            <  1:841136/150‑1 (MQ=255)
                                                                                                        aaGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAAc                                         <  1:218594/150‑1 (MQ=255)
                                                                                                        aaGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAAc                                         >  2:1114550/1‑150 (MQ=255)
                                                                                                        aaGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAAc                                         <  1:408921/150‑1 (MQ=255)
                                                                                                         aGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACg                                        <  1:1627974/150‑1 (MQ=255)
                                                                                                            tGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCGCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTg                                     <  1:324755/150‑1 (MQ=255)
                                                                                                             ggCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTACTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGa                                    >  2:852600/1‑150 (MQ=255)
                                                                                                               cgcAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGc                                  <  2:162255/150‑1 (MQ=255)
                                                                                                                gcAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACATCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCAAGCATTCCTGATGGATTTTCATGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCt                                 <  1:485458/150‑1 (MQ=255)
                                                                                                                gcAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCAGGAAGAGTTTTACGCGCGTTATAACGTTGAGCt                                 >  2:77220/1‑150 (MQ=255)
                                                                                                                 cAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCtt                                <  2:632250/150‑1 (MQ=255)
                                                                                                                  aGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCATGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTa                               <  1:1715778/150‑1 (MQ=255)
                                                                                                                  aGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTa                               <  2:1434328/150‑1 (MQ=255)
                                                                                                                  aGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTa                               <  2:1440936/150‑1 (MQ=255)
                                                                                                                      gcTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTcg                           >  2:1327643/1‑150 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                        >  2:96628/1‑150 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                        >  2:748071/1‑150 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                        >  2:589072/1‑150 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                        >  2:1141090/1‑150 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                        <  1:576791/150‑1 (MQ=255)
                                                                                                                          gCGAATCCGTTAATTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCa                       <  2:1477841/150‑1 (MQ=255)
                                                                                                                                 cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAg                >  1:1265189/1‑150 (MQ=255)
                                                                                                                                 cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAg                <  2:1701331/150‑1 (MQ=255)
                                                                                                                                  gTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAgg               >  2:982439/1‑150 (MQ=255)
                                                                                                                                       gtaaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGATCTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATTACGTTGAGCTTATTCGCGCACCAGAAGGGttgt          >  1:130828/5‑148 (MQ=255)
                                                                                                                                       gtaaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATTATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTAGCGCACCAGAAGGGttct          >  1:671665/5‑150 (MQ=255)
                                                                                                                                        taaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttctt         >  2:439170/4‑150 (MQ=255)
                                                                                                                                        tCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCa                    <  2:942381/139‑1 (MQ=255)
                                                                                                                                        tCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCa                    >  1:942381/1‑139 (MQ=255)
                                                                                                                                         ccGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttcttc        <  1:1614860/150‑1 (MQ=255)
                                                                                                                                           ggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTa      >  1:1113435/1‑150 (MQ=255)
                                                                                                                                            gcgcTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTAt     <  2:173983/150‑1 (MQ=255)
                                                                                                                                               cTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTa  >  1:495940/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTA  >  NC_000913/975554‑975846

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: