Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,443,594 G→T 8.0% T168T (ACC→ACA secY ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,5940GT 8.0% 254.3 / 7.2 74T168T (ACC→ACAsecYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (26/42);  new base T (2/4);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGG  >  NC_000913/3443446‑3443739
                                                                                                                                                    |                                                                                                                                                 
ggAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATg                                                                                                                                                   <  1:608974/149‑1 (MQ=255)
ggAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCACAACCACATCAGGAACATg                                                                                                                                                   <  1:323412/149‑1 (MQ=255)
 gAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATgg                                                                                                                                                  >  3:343585/1‑149 (MQ=255)
  aaGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGt                                                                                                                                                 <  3:485517/149‑1 (MQ=255)
   aGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGtt                                                                                                                                                >  1:532533/1‑149 (MQ=255)
        aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  1:544839/149‑2 (MQ=255)
        aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  3:626337/149‑2 (MQ=255)
        aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  1:59898/149‑2 (MQ=255)
        aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtgt                                                                                                                                           <  4:59052/149‑1 (MQ=255)
        aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtgt                                                                                                                                           >  2:52882/1‑149 (MQ=255)
                 tGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGCTCCTGTGACCAGACt                                                                                                                                  <  1:433830/149‑1 (MQ=255)
                   aCGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTa                                                                                                                                <  3:240681/149‑1 (MQ=255)
                   aCGCGCTTGCTCGAGAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTa                                                                                                                                <  4:343585/149‑1 (MQ=255)
                      cgcTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTacaa                                                                                                                             <  3:270752/149‑1 (MQ=255)
                           gCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCg                                                                                                                        <  4:585699/149‑1 (MQ=255)
                             tCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGt                                                                                                                      <  4:475941/149‑1 (MQ=255)
                              cGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTg                                                                                                                     <  3:197506/149‑1 (MQ=255)
                               gATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTacaa                                                                                                                             >  4:681970/1‑140 (MQ=255)
                               gATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTacaa                                                                                                                             <  3:681970/140‑1 (MQ=255)
                               gATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGa                                                                                                                    >  3:637194/1‑149 (MQ=255)
                                aTAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGaa                                                                                                                   <  2:376556/149‑1 (MQ=255)
                                   gTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTa                                                                                                                <  3:224949/149‑1 (MQ=255)
                                    tATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAAGGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAg                                                                                                               >  3:555114/1‑149 (MQ=255)
                                      tGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGt                                                                                                                      >  4:5013/1‑140 (MQ=255)
                                      tGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGt                                                                                                                      <  3:5013/140‑1 (MQ=255)
                                      tGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGaa                                                                                                             <  4:342307/149‑1 (MQ=255)
                                      tGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGaa                                                                                                             <  1:576902/149‑1 (MQ=255)
                                             tggctggcGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAg                                                                                                      >  4:324131/1‑149 (MQ=255)
                                                  ggcGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCgg                                                                                                 <  2:532533/149‑1 (MQ=255)
                                                  ggcGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCgg                                                                                                 <  2:53344/149‑1 (MQ=255)
                                                         gTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATc                                                                                          <  3:669761/149‑1 (MQ=255)
                                                          tCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCa                                                                                         <  1:590373/149‑1 (MQ=255)
                                                          tCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCa                                                                                         <  4:577433/149‑1 (MQ=255)
                                                           cccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAc                                                                                        <  2:435237/149‑1 (MQ=255)
                                                            ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGc                                                                                                                         >  4:3495/1‑115 (MQ=255)
                                                            ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGc                                                                                                                         <  3:3495/115‑1 (MQ=255)
                                                            ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAcc                                                                                       <  3:16312/149‑1 (MQ=255)
                                                            ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAcc                                                                                       >  1:112956/1‑149 (MQ=255)
                                                            ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAcc                                                                                       >  1:420077/1‑149 (MQ=255)
                                                             cgcgACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCa                                                                                      <  4:273664/149‑1 (MQ=255)
                                                                   aaTACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCtt                                                                                <  3:204252/149‑1 (MQ=255)
                                                                   aaTACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCtt                                                                                <  4:600346/149‑1 (MQ=255)
                                                                         ggagAAGATAATGATTGAAAGACCGTTGCCGATACCTAGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGAACAGACTTACAACAGCAGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATAc                                                                          <  2:548070/146‑1 (MQ=255)
                                                                           cGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCa                                                                        <  3:224474/149‑1 (MQ=255)
                                                                                      aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCgg                                                             >  4:421217/1‑149 (MQ=255)
                                                                                      aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCgg                                                             >  1:14821/1‑149 (MQ=255)
                                                                                       atcagaTACCGTTCCCGATACCCCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGCTTCCTGTGACCAGACTTACAACAGCGGTGAAATAGAATACAAAGCCCGGGCTAATCACCAGTCCTTGCATACCACGCATATTCGGc                                                            >  4:323461/6‑149 (MQ=255)
                                                                                          aaaTACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGa                                                         >  1:625090/1‑149 (MQ=255)
                                                                                              aCCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCgg                                                     <  1:684094/149‑1 (MQ=255)
                                                                                                cGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTa                                                   >  1:91337/1‑149 (MQ=255)
                                                                                                 gTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAg                                                  >  2:572289/1‑149 (MQ=255)
                                                                                                        aTACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATAcc                                           <  3:386760/149‑1 (MQ=255)
                                                                                                           ccTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGAt                                        <  2:420077/149‑1 (MQ=255)
                                                                                                            cTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATc                                       >  2:614606/1‑149 (MQ=255)
                                                                                                               gTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGAc                                    <  1:52882/149‑1 (MQ=255)
                                                                                                               gTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGAc                                    <  1:566838/149‑1 (MQ=255)
                                                                                                               gTTCAGTAATCGGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGAc                                    >  3:393943/1‑149 (MQ=255)
                                                                                                                 tCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTg                                  <  4:637194/149‑1 (MQ=255)
                                                                                                                    gTAATCTGTTCGCCCACCCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGaa                               <  4:555114/149‑1 (MQ=255)
                                                                                                                     tAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGc                                                                     <  1:466602/110‑1 (MQ=255)
                                                                                                                     tAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGc                                                                     >  2:466602/1‑110 (MQ=255)
                                                                                                                        tCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATAtt                           >  2:239633/1‑149 (MQ=255)
                                                                                                                            ttCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCa                       >  1:494799/1‑149 (MQ=255)
                                                                                                                               gCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGt                                                    >  3:391484/1‑117 (MQ=255)
                                                                                                                               gCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGt                                                    <  4:391484/117‑1 (MQ=255)
                                                                                                                               gCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGACAGCa                    >  2:197429/1‑149 (MQ=255)
                                                                                                                                  cAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCa                 <  4:393943/149‑1 (MQ=255)
                                                                                                                                       acaTCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTa            <  2:574728/149‑1 (MQ=255)
                                                                                                                                        tcacAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATAcc                                                                         >  3:606045/4‑87 (MQ=255)
                                                                                                                                        tcacAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATAcc                                                                         <  4:606045/84‑1 (MQ=255)
                                                                                                                                         cacAGGAACCATGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTAcc          >  4:2012/3‑149 (MQ=255)
                                                                                                                                         aTCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGcc                        <  2:141104/135‑1 (MQ=255)
                                                                                                                                         aTCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGcc                        >  1:141104/1‑135 (MQ=255)
                                                                                                                                          tCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCg                                                                                                  <  4:156450/60‑1 (MQ=255)
                                                                                                                                          tCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCg                                                                                                  >  3:156450/1‑60 (MQ=255)
                                                                                                                                                aCATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCgg   <  3:324131/149‑1 (MQ=255)
                                                                                                                                                 cATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCggg  >  2:513485/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGG  >  NC_000913/3443446‑3443739

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: