Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,786,902 G→A 5.8% S22N (AGC→AAC)  envC → murein hydrolase activator EnvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,9020GA 5.8% 491.0 / 10.8 154S22N (AGC→AAC) envCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base G (78/66);  new base A (2/7);  total (81/73)
Fisher's exact test for biased strand distribution p-value = 8.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGG  >  NC_000913/3786753‑3787051
                                                                                                                                                     |                                                                                                                                                     
tATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGTATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAg                                                                                                                                                       >  1:212515/1‑150 (MQ=255)
 aTGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGATATTTTTTATTTTTTCAATGAGGACCATGATCTACACCAGCGTTCTTAgc                                                                                                                                                      <  2:75842/150‑1 (MQ=255)
   ggAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAgcgc                                                                                                                                                    <  2:359869/150‑1 (MQ=255)
     aaaTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCTATCGTTCTTAGCGCtt                                                                                                                                                  <  2:350023/150‑2 (MQ=255)
     aaaTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTg                                                                                                                                                  >  2:521989/1‑150 (MQ=255)
       aTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGc                                                                                                                                                <  1:1196485/150‑1 (MQ=255)
         cccGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGt                                                                                                                                              <  1:533997/150‑1 (MQ=255)
         cccGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGt                                                                                                                                              <  1:587565/150‑1 (MQ=255)
         cccGCAAGAGATGACTGGTAAGCAGCTGTTCATCGTGGAATAATCACTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCTAAGCTTTGCAATCAGGCACATCATCTACGCCAGCGTTCTTATCGCTGGCGt                                                                                                                                              >  1:914716/1‑150 (MQ=255)
          ccGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGTa                                                                                                                                             <  1:529154/150‑1 (MQ=255)
           cccAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGCGAAAGGCGATTAATAAAATGACACGGGCCGTGAAACCGCGCAGCTATGCAATCAGGCCCATCATCTACGCCAGCGTTCATAGCGCTcgcgaaa                                                                                                                                            >  2:170572/3‑147 (MQ=255)
           cGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAt                                                                                                                                            >  2:279231/1‑150 (MQ=255)
           cGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGTag                                                                                                                                            <  2:287404/150‑2 (MQ=255)
              aaGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAttgt                                                                                                                                         <  2:1394148/150‑1 (MQ=255)
               agagATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAttgttg                                                                                                                                       >  2:549214/1‑150 (MQ=255)
                gagaTGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAttgttg                                                                                                                                       <  2:202514/150‑1 (MQ=255)
                gagaTGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAttgttg                                                                                                                                       >  2:1193964/1‑150 (MQ=255)
                  gaTGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTtgtg                                                                                                                                     <  2:604500/150‑1 (MQ=255)
                          gTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCtttt                                                                                                                             <  2:1042472/150‑1 (MQ=255)
                               ccGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGcc                                                                                                                        >  2:1335422/1‑150 (MQ=255)
                               ccGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGcc                                                                                                                        <  1:1409629/150‑1 (MQ=255)
                               ccGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGcc                                                                                                                        <  1:626271/150‑1 (MQ=255)
                                     ttCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCAcgcg                                                                                                                  <  1:15323/150‑1 (MQ=255)
                                       cATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGa                                                                                                                >  2:471713/1‑150 (MQ=255)
                                       cATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGa                                                                                                                <  2:974778/150‑1 (MQ=255)
                                         tCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATg                                                                                                              <  1:279231/150‑1 (MQ=255)
                                          cGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGa                                                                                                             >  1:161829/1‑150 (MQ=255)
                                           gTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAg                                                                                                            >  1:135927/1‑150 (MQ=255)
                                           gTGGAATAATCCCTCCCCATGAGGGGAAAGGCAATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAg                                                                                                            >  2:1377849/1‑150 (MQ=255)
                                               aataatCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTg                                                                                                        >  1:1146245/1‑150 (MQ=255)
                                                  aatCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGAcc                                                                                                     <  1:427208/150‑1 (MQ=255)
                                                   atCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCa                                                                                                    >  1:230913/1‑150 (MQ=255)
                                                      cctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACt                                                                                                 >  2:326240/1‑150 (MQ=255)
                                                      cctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACt                                                                                                 >  1:452570/1‑150 (MQ=255)
                                                       ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGACTTTTCCGCCCAGGCGGATGAGCGTGACCAACTc                                                                                                <  2:546556/150‑1 (MQ=255)
                                                            cATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATc                                                                                           <  1:146650/150‑1 (MQ=255)
                                                            cATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTCGCGTATTGTTGTGCGCCTTTTCCTCCCACGCGGATCAGCGTGACCAACTCAAATc                                                                                           >  2:1021352/1‑150 (MQ=255)
                                                             aTGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCt                                                                                          <  2:510731/150‑1 (MQ=255)
                                                              tGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTa                                                                                         >  1:679678/1‑150 (MQ=255)
                                                               gAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTAt                                                                                        >  2:150984/1‑150 (MQ=255)
                                                                   ggAAAGGCGATTAGTACCATGAAACGGGCCGTGAACCCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAATTCAAATCTATTCAg                                                                                    >  2:1254349/1‑150 (MQ=255)
                                                                   ggAAAGGCGATTAATACCATGCCACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAg                                                                                    >  1:1035501/1‑150 (MQ=255)
                                                                   ggAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAg                                                                                    <  2:314877/150‑1 (MQ=255)
                                                                    gAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAgg                                                                                   <  1:846983/150‑1 (MQ=255)
                                                                     aaaGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGc                                                                                  >  1:3465/1‑150 (MQ=255)
                                                                     aaaGGCAATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGc                                                                                  <  1:1377849/150‑1 (MQ=255)
                                                                       aGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTg                                                                                >  2:949287/1‑150 (MQ=255)
                                                                        ggCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGa                                                                               >  1:64137/1‑150 (MQ=255)
                                                                           gATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACAt                                                                            <  2:380280/150‑1 (MQ=255)
                                                                              tAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATcgc                                                                         >  2:1069234/1‑150 (MQ=255)
                                                                                aTACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATcgccg                                                                       <  1:471713/150‑1 (MQ=255)
                                                                                 tACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATcgccgc                                                                      <  1:694464/150‑1 (MQ=255)
                                                                                  aCCATGACCCGGGCCGTGAAACCCCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCa                                                                     <  2:156106/150‑1 (MQ=255)
                                                                                   ccATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCaa                                                                    >  1:18474/1‑150 (MQ=255)
                                                                                   ccATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCaa                                                                    >  1:879561/1‑150 (MQ=255)
                                                                                   ccATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCaa                                                                    >  1:13149/1‑150 (MQ=255)
                                                                                      tGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAag                                                                 >  2:1357369/1‑150 (MQ=255)
                                                                                        acacGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAagag                                                               <  1:326240/150‑1 (MQ=255)
                                                                                             ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAgcgcgc                                                          <  2:64137/150‑1 (MQ=255)
                                                                                              gCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGcc                                                                                                                        >  1:1094819/1‑87 (MQ=255)
                                                                                              gCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGcc                                                                                                                        <  2:1089607/87‑1 (MQ=255)
                                                                                                cGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGt                                                       <  1:537546/150‑1 (MQ=255)
                                                                                                 ttGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTa                                                      >  2:594944/2‑150 (MQ=255)
                                                                                                 gTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTa                                                      >  2:851804/1‑150 (MQ=255)
                                                                                                  tGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTAc                                                     >  2:841671/1‑150 (MQ=255)
                                                                                                      aCCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCa                                                 >  2:1204894/1‑150 (MQ=255)
                                                                                                        cgcgCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCaaa                                               <  1:29203/150‑1 (MQ=255)
                                                                                                         gcgcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCaaaa                                              <  2:496649/150‑1 (MQ=255)
                                                                                                         gcgcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCaaaa                                              <  2:872166/150‑1 (MQ=255)
                                                                                                         gcgcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCaaaa                                              >  1:929034/1‑150 (MQ=255)
                                                                                                          ctcAGGTTTGAAATCATGCCCATCATCTACGCCATCGTTATTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGTCGAAAAAGAGCGCGCGGTACGCCAAAAg                                             <  2:914716/148‑1 (MQ=255)
                                                                                                          cgcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAg                                             <  2:822716/150‑1 (MQ=255)
                                                                                                           gcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGc                                            <  2:230913/150‑1 (MQ=255)
                                                                                                            cAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGca                                           <  1:910505/150‑1 (MQ=255)
                                                                                                              ggTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaac                                         <  1:1281660/150‑1 (MQ=255)
                                                                                                               gTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaaca                                        <  1:1052765/150‑1 (MQ=255)
                                                                                                                   gCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGCCCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaacaa                                    >  2:1261819/1‑150 (MQ=255)
                                                                                                                   gCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaacaa                                    >  1:578311/1‑150 (MQ=255)
                                                                                                                   gCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaacaa                                    <  2:451371/150‑1 (MQ=255)
                                                                                                                   gCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCCCTGCCGTATTGTACTGAGCCTTTTCCGCCCACGCGGATCAGACTGACCAACTCAAATCTATTCAGGCTGACATCCACGCAAAAGAGCGCCCGGTACGACAAAAGcaacaacaa                                    >  1:1026716/1‑150 (MQ=255)
                                                                                                                      aTCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAAcgc                                 >  2:484267/1‑150 (MQ=255)
                                                                                                                         aGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCa                              <  1:1323529/150‑1 (MQ=255)
                                                                                                                          ggCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCaa                             >  1:493060/1‑150 (MQ=255)
                                                                                                                              catcatCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCt                         <  2:879561/150‑1 (MQ=255)
                                                                                                                              catcatCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTg                        >  1:98010/1‑150 (MQ=255)
                                                                                                                               atcatcTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTg                        <  1:1069234/150‑1 (MQ=255)
                                                                                                                                   tcTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCg                    >  1:797957/1‑150 (MQ=255)
                                                                                                                                   aatACGCCAGCGCTAAGAACGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGc                       <  1:1342384/145‑1 (MQ=255)
                                                                                                                                   aatACGCCAGCGCTAAGAACGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGc                       >  2:1342384/3‑147 (MQ=255)
                                                                                                                                       cGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAg                                             >  2:1108893/1‑121 (MQ=255)
                                                                                                                                       cGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGcaca                <  2:1322167/150‑1 (MQ=255)
                                                                                                                                       cGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGcaca                >  1:650663/1‑150 (MQ=255)
                                                                                                                                         ccAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAAt              >  2:565448/1‑150 (MQ=255)
                                                                                                                                           aGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTc                     >  1:416173/1‑141 (MQ=255)
                                                                                                                                           aGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTc                     <  2:416173/141‑1 (MQ=255)
                                                                                                                                           aGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTg            <  2:761242/150‑1 (MQ=255)
                                                                                                                                           aGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTg            >  1:1269791/1‑150 (MQ=255)
                                                                                                                                                  ttAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGc     <  2:1046747/150‑1 (MQ=255)
                                                                                                                                                   tAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCa    >  1:508722/1‑150 (MQ=255)
                                                                                                                                                    aGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAg   <  1:1146370/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAAAGCAGG  >  NC_000913/3786753‑3787051

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: