Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,443,592 A→C 7.7% M169R (ATG→AGG) ‡ secY ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,5920AC 7.7% 259.9 / 12.6 78M169R (ATG→AGG) ‡secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (27/45);  new base C (2/4);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GAGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGG  >  NC_000913/3443444‑3443739
                                                                                                                                                    |                                                                                                                                                   
gAGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACa                                                                                                                                                     <  4:149740/149‑1 (MQ=255)
 aGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACAt                                                                                                                                                    <  3:443902/149‑1 (MQ=255)
  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:59898/149‑2 (MQ=255)
          aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  3:626337/149‑2 (MQ=255)
          aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  1:544839/149‑2 (MQ=255)
          aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtgt                                                                                                                                           >  2:52882/1‑149 (MQ=255)
          aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCtgt                                                                                                                                           <  4:59052/149‑1 (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                                                                                                                      <  3:5013/140‑1 (MQ=255)
                                        tGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGt                                                                                                                      >  4:5013/1‑140 (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                                                                                                                         <  3:3495/115‑1 (MQ=255)
                                                              ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGc                                                                                                                         >  4:3495/1‑115 (MQ=255)
                                                              ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAcc                                                                                       >  1:112956/1‑149 (MQ=255)
                                                              ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAcc                                                                                       <  3:16312/149‑1 (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                                                             >  1:14821/1‑149 (MQ=255)
                                                                                        aTTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCgg                                                             >  4:421217/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:566838/149‑1 (MQ=255)
                                                                                                                 gTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGAc                                    <  1:52882/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                                                    <  4:391484/117‑1 (MQ=255)
                                                                                                                                 gCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGt                                                    >  3:391484/1‑117 (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)
                                                                                                                                                    |                                                                                                                                                   
GAGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGG  >  NC_000913/3443444‑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: