Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,603,724 A→C 12.8% T236P (ACC→CCC)  lsrC → Autoinducer‑2 ABC transporter membrane subunit LsrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,7240AC12.8% 184.6 / 19.4 77T236P (ACC→CCC) lsrCAutoinducer‑2 ABC transporter membrane subunit LsrC
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base C (5/5);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTC  >  NC_000913/1603576‑1603871
                                                                                                                                                    |                                                                                                                                                   
ttATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGAc                                                                                                                                                    <  1:27106/150‑1 (MQ=255)
 tATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGAcc                                                                                                                                                   <  1:1378973/150‑1 (MQ=255)
       aCGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTAcc                                                                                                                                             >  2:11960/1‑150 (MQ=255)
            cGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                        >  2:1116257/1‑150 (MQ=255)
             gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTg                                                                                                                                       <  2:261882/150‑1 (MQ=255)
              aTAATTTACTGGTCGCTCGTCAACTGGGCGTTCGTACTGTAGCCATTCGCATTGTGGCTTTTTCTTTTAACGGCTGCATGGCGGCACTGGCGTGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTgg                                                                                                                                      <  2:1525868/150‑1 (MQ=255)
              aTAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTgg                                                                                                                                      <  1:587781/150‑1 (MQ=255)
               tAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggt                                                                                                                                     <  1:497816/150‑2 (MQ=255)
               tAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                     >  2:1115526/1‑150 (MQ=255)
                aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggtt                                                                                                                                    <  1:742974/150‑3 (MQ=255)
                aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggtt                                                                                                                                    <  2:225905/150‑3 (MQ=255)
                 aTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttg                                                                                                                                   >  1:1309826/1‑147 (MQ=255)
                   ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttggg                                                                                                                                 <  2:63407/150‑6 (MQ=255)
                   ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                 <  1:25678/150‑1 (MQ=255)
                   ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                 <  2:1046597/150‑1 (MQ=255)
                    tACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttgggg                                                                                                                                <  1:50983/150‑7 (MQ=255)
                      cAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGaaa                                                                                                                              >  1:780927/1‑150 (MQ=255)
                       aGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAg                                                                                                                             <  1:1114899/150‑1 (MQ=255)
                        gggCGCTCGTGAACTGGGCGTTCGTACTGAATGCATTCTCATTTTGGCATTTTCGTTGAACGGCTGCATGGTGTGACTTTCTGGAATTGTGTTTGCTTCTCAGATTGGTTTTATCCGGAACCTGAGCGGTACCTGGCTGGTGATGAAAGc                                                                                                                            <  2:101702/150‑1 (MQ=255)
                        gggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGc                                                                                                                            <  1:623512/150‑1 (MQ=255)
                         ggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCa                                                                                                                           <  2:1163113/150‑1 (MQ=255)
                               cGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                     >  2:929287/1‑150 (MQ=255)
                                gTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAg                                                                                                                    >  1:608391/1‑150 (MQ=255)
                                 tCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGc                                                                                                                   >  1:1363192/1‑150 (MQ=255)
                                  cAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGcc                                                                                                                  >  2:23606/1‑150 (MQ=255)
                                  cAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGcc                                                                                                                  >  2:33630/1‑150 (MQ=255)
                                       gggCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGt                                                                                                             >  1:936258/1‑150 (MQ=255)
                                          cGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCt                                                                                                          >  1:475229/1‑150 (MQ=255)
                                           gTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTg                                                                                                         <  1:23606/150‑1 (MQ=255)
                                              cGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggc                                                                                                      <  2:1463634/150‑1 (MQ=255)
                                                  cTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                  >  1:567010/1‑150 (MQ=255)
                                                      aGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGAGCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAg                                                                                              <  1:3166/150‑1 (MQ=255)
                                                                  tgtgGCATTTTCGTTGAACGTCTGCATGGCGGCTCTGGCGGGATTTGTGTTTGCTTGGCAGATTGGTTTTTTCCCCATCCAGACCGGTACCGGGCTGGTGATGAAAGCTATTTCTGGCTGCTTGCTGGGCGGCTTTAGTTTGGTCgttgg                                                                                  <  2:446728/150‑1 (MQ=255)
                                                                  tgtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtgg                                                                                  >  1:272400/1‑150 (MQ=255)
                                                                  ggtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtgg                                                                                  <  1:627567/149‑1 (MQ=255)
                                                                   gtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtggt                                                                                 >  2:636829/1‑150 (MQ=255)
                                                                   gtgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCggtggt                                                                                 <  2:584756/150‑1 (MQ=255)
                                                                       cATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCAAACCAGCCCGGTACAGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCg                                                                             >  1:1149793/1‑150 (MQ=255)
                                                                           ttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGc                                                                         >  1:81711/1‑150 (MQ=255)
                                                                           ttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGc                                                                         <  2:475229/150‑1 (MQ=255)
                                                                                  aaCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCAtt                                                                  <  2:272400/150‑1 (MQ=255)
                                                                                    cGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTgg                                                                <  1:1335710/150‑1 (MQ=255)
                                                                                         gCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCgg                                                           <  1:1115526/150‑1 (MQ=255)
                                                                                          cATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGt                                                          <  2:1024248/150‑1 (MQ=255)
                                                                                            tGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                            >  1:1337767/1‑114 (MQ=255)
                                                                                            tGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                            <  2:1337767/114‑1 (MQ=255)
                                                                                                   tcTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGgcgc                                                 >  2:1188464/2‑150 (MQ=255)
                                                                                                      ggCGGGATTTGTGTTTGCTTCGCAGATTGGGTTTATCCCCAACCAGACCGGTAGCGGGCTGGAGATGAATGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGGGCGGTACTCGGCGCATg                                              <  2:1509628/150‑1 (MQ=255)
                                                                                                      ggCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATg                                              >  2:1405455/1‑150 (MQ=255)
                                                                                                      ggCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATg                                              >  1:276006/1‑150 (MQ=255)
                                                                                                      ggCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCACCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATg                                              >  2:1024440/1‑150 (MQ=255)
                                                                                                      ggCGGGAATTGTGTTTGCTTCGCAGATTCGTGTTATCCCCAACCAGACCGGCACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCAGTAGTTTGCTCGGTGGTTCCGGGGCGATCCTTGGTGCGGTACTCGGCGCATg                                              <  2:156226/150‑1 (MQ=255)
                                                                                                         gggAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGtt                                           <  1:929287/150‑1 (MQ=255)
                                                                                                             aTTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTg                                       <  1:1116257/150‑1 (MQ=255)
                                                                                                                gtgtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                    <  2:936258/150‑1 (MQ=255)
                                                                                                                gtgtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                    >  2:1152788/1‑150 (MQ=255)
                                                                                                                  gtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTCGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCa                                  <  2:780927/150‑1 (MQ=255)
                                                                                                                    ttGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGa                                >  2:847558/1‑150 (MQ=255)
                                                                                                                     tGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGag                               >  2:577835/1‑149 (MQ=255)
                                                                                                                            cAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGAGGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGGGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGc                        <  1:1518364/150‑1 (MQ=255)
                                                                                                                             aGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCg                       <  1:11960/150‑1 (MQ=255)
                                                                                                                                ttGGTTTTATCCCCAACCAGCCCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGAGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGAGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGAGTAc                    >  2:227263/1‑150 (MQ=255)
                                                                                                                                     tttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTg               >  1:14478/1‑150 (MQ=255)
                                                                                                                                        aTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTg            >  2:1467673/1‑150 (MQ=255)
                                                                                                                                         tCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGt           >  1:1541583/1‑150 (MQ=255)
                                                                                                                                          ccccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGtt          >  1:1291774/1‑150 (MQ=255)
                                                                                                                                           cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  1:1435032/1‑150 (MQ=255)
                                                                                                                                           cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  1:1429295/1‑150 (MQ=255)
                                                                                                                                           cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  1:831743/1‑150 (MQ=255)
                                                                                                                                           cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  1:860612/1‑150 (MQ=255)
                                                                                                                                           catctccAGCCCGGTACCGGTCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg         >  1:87134/6‑150 (MQ=255)
                                                                                                                                            atctccAGCCCGGTACCGGTCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTAGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgc        >  2:343229/5‑150 (MQ=255)
                                                                                                                                             cAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgcg       >  2:49125/1‑150 (MQ=255)
                                                                                                                                             cAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgcg       <  2:81711/150‑1 (MQ=255)
                                                                                                                                               aCCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCa     >  2:854369/1‑150 (MQ=255)
                                                                                                                                                 cAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCAtt   >  1:61924/1‑150 (MQ=255)
                                                                                                                                                  aGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTc  >  1:994124/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTC  >  NC_000913/1603576‑1603871

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