Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 975,701 A→C 6.0% D27A (GAC→GCC)  mukE → chromosome partitioning protein MukE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913975,7010AC 6.0% 251.1 / 7.6 83D27A (GAC→GCC) mukEchromosome partitioning protein MukE
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base C (2/3);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTACG  >  NC_000913/975554‑975848
                                                                                                                                                   |                                                                                                                                                   
aGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACa                                                                                                                                                   >  1:736227/1‑150 (MQ=255)
    gATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgcgc                                                                                                                                               >  1:1366570/1‑150 (MQ=255)
      tttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCt                                                                                                                                             >  2:1607267/1‑150 (MQ=255)
       ttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAGCGccgg                                                                                                                                            <  1:633067/150‑3 (MQ=255)
       ttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCtt                                                                                                                                            <  2:736227/150‑1 (MQ=255)
        tCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAGCGccgga                                                                                                                                           <  1:1344591/150‑4 (MQ=255)
         cACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                                                                                                                                          >  2:384205/1‑150 (MQ=255)
           ccGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCAGGCGCTGTCCAGCGccggaaat                                                                                                                                        <  2:212977/150‑7 (MQ=255)
              gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCa                                                                                                                                     <  1:477792/150‑1 (MQ=255)
              gACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCa                                                                                                                                     <  2:1005573/150‑1 (MQ=255)
               aCTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAg                                                                                                                                    <  2:897175/150‑1 (MQ=255)
               aCTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAg                                                                                                                                    <  1:1431114/150‑1 (MQ=255)
                  gCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGAc                                                                                                                                 >  2:1007606/1‑150 (MQ=255)
                   ccAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                                                                                                                                >  2:1003835/1‑150 (MQ=255)
                    cAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                               >  2:1562894/1‑150 (MQ=255)
                     aGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGcc                                                                                                                              >  1:1460407/1‑150 (MQ=255)
                            tGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTgg                                                                                                                       <  1:700089/150‑1 (MQ=255)
                            tGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTgg                                                                                                                       <  2:1457071/150‑1 (MQ=255)
                                      ttAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGaa                                                                                                             <  2:1143971/150‑1 (MQ=255)
                                        aaTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACt                                                                                                           >  2:1068582/1‑150 (MQ=255)
                                        aaTGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACt                                                                                                           >  2:395095/1‑150 (MQ=255)
                                          tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTgg                                                                                                         >  1:1028711/1‑150 (MQ=255)
                                           gATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGa                                                                                                        <  2:1432674/150‑1 (MQ=255)
                                            aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               <  2:638653/142‑1 (MQ=255)
                                            aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               <  2:644160/142‑1 (MQ=255)
                                            aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               >  1:644160/1‑142 (MQ=255)
                                            aTTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               >  1:638653/1‑142 (MQ=255)
                                                cGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATc                                                                                                   >  2:1392752/1‑150 (MQ=255)
                                                      cAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCAt                                                                                             >  1:206459/1‑150 (MQ=255)
                                                               ggCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATgg                                                                                    >  1:38839/1‑150 (MQ=255)
                                                                 cgcATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGAt                                                                                  >  2:602880/1‑150 (MQ=255)
                                                                 cgcATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGAt                                                                                  >  2:1578333/1‑150 (MQ=255)
                                                                     tGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTc                                                                              <  1:386869/150‑1 (MQ=255)
                                                                     tGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTc                                                                              <  2:1429344/150‑1 (MQ=255)
                                                                        cATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               <  1:290244/114‑1 (MQ=255)
                                                                        cATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgacg                                                                                                               >  2:290244/1‑114 (MQ=255)
                                                                        cATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAgg                                                                           <  1:153305/150‑1 (MQ=255)
                                                                        cATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAgg                                                                           >  2:665482/1‑150 (MQ=255)
                                                                            gACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAgg                                                                           >  2:159174/1‑146 (MQ=255)
                                                                            gACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAgg                                                                           <  1:159174/146‑1 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAg                                                                            <  1:963383/143‑1 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAg                                                                            >  2:963383/1‑143 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     >  1:432283/1‑150 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     >  2:658245/1‑150 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     <  2:1460407/150‑1 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     >  1:1039455/1‑150 (MQ=255)
                                                                              cAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATc                                                                     <  1:1025632/150‑1 (MQ=255)
                                                                                    ttGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAag                                                               >  1:1382343/1‑150 (MQ=255)
                                                                                    ttGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAag                                                               <  2:481949/150‑1 (MQ=255)
                                                                                          aaGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTTCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGtttt                                                         <  2:969941/150‑1 (MQ=255)
                                                                                                tGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgc                                                   <  1:222444/150‑1 (MQ=255)
                                                                                                tGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgc                                                   >  1:966725/1‑150 (MQ=255)
                                                                                                       tAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATaa                                            <  1:272821/150‑1 (MQ=255)
                                                                                                         aGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACg                                          >  1:931333/1‑150 (MQ=255)
                                                                                                          gCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAata                                                                       <  2:1260267/120‑1 (MQ=255)
                                                                                                          gCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAata                                                                       >  1:1260267/1‑120 (MQ=255)
                                                                                                           cTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGtt                                        >  2:1495433/1‑150 (MQ=255)
                                                                                                            tGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTg                                       <  2:1039455/150‑1 (MQ=255)
                                                                                                            tGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTg                                       <  2:206459/150‑1 (MQ=255)
                                                                                                             ggCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGa                                      <  2:839863/150‑1 (MQ=255)
                                                                                                             ggCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGa                                      >  2:685303/1‑150 (MQ=255)
                                                                                                                gcAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCt                                   <  1:679443/150‑1 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTg                                       >  2:513996/1‑137 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTg                                       <  1:513996/137‑1 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                          >  2:908371/1‑150 (MQ=255)
                                                                                                                         ggCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCgcgc                          >  1:1074175/1‑150 (MQ=255)
                                                                                                                           cGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCAc                        >  1:923152/1‑150 (MQ=255)
                                                                                                                           cGAATCCATTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCAc                        >  2:1008222/1‑150 (MQ=255)
                                                                                                                            gAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCAcc                       <  1:1607267/150‑1 (MQ=255)
                                                                                                                                ccGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGaa                   >  1:1117788/1‑150 (MQ=255)
                                                                                                                                      tttCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttc             >  2:1566155/1‑150 (MQ=255)
                                                                                                                                       ttCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttct            >  1:1302562/1‑150 (MQ=255)
                                                                                                                                        taaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttctt           >  2:1360597/4‑150 (MQ=255)
                                                                                                                                        taaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttctt           >  1:1316208/4‑150 (MQ=255)
                                                                                                                                         ccGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCTCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttcttc          <  1:1003835/150‑1 (MQ=255)
                                                                                                                                         ccGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttcttc          >  2:583443/1‑150 (MQ=255)
                                                                                                                                          cGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTtcttct         >  2:658556/1‑150 (MQ=255)
                                                                                                                                            gcgcTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTAt       <  2:1028711/150‑1 (MQ=255)
                                                                                                                                             cgcTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTAtt      <  1:1068582/150‑1 (MQ=255)
                                                                                                                                               cTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTa    <  2:256179/150‑1 (MQ=255)
                                                                                                                                                tGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTAc   <  2:462342/150‑1 (MQ=255)
                                                                                                                                                tGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTAc   >  2:1117813/1‑150 (MQ=255)
                                                                                                                                                 ggACAGCGGCTTACGCTCAGCACGCCATATTGGCCGCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGAATTACGCGCGTTATCACGTTGACCTTATTCGCGCACCACCAGCGTTCTTCTATTTACg  <  2:1359406/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTACG  >  NC_000913/975554‑975848

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