Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,443,593 T→A 7.6% M169L (ATG→TTG) ‡ secY ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,5930TA 7.6% 265.0 / 7.8 77M169L (ATG→TTG) ‡secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (27/44);  new base A (2/4);  total (29/48)
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

AGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGG  >  NC_000913/3443445‑3443739
                                                                                                                                                    |                                                                                                                                                  
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:544839/149‑2 (MQ=255)
         aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  3:626337/149‑2 (MQ=255)
         aGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACCATGTTCCtga                                                                                                                                           <  1:59898/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                                                                                                                             <  3:681970/140‑1 (MQ=255)
                                gATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTacaa                                                                                                                             >  4:681970/1‑140 (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                                                                                         <  4:577433/149‑1 (MQ=255)
                                                           tCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCa                                                                                         <  1:590373/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:420077/1‑149 (MQ=255)
                                                             ccGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCAcc                                                                                       >  1:112956/1‑149 (MQ=255)
                                                              cgcgACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCa                                                                                      <  4:273664/149‑1 (MQ=255)
                                                                    aaTACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCtt                                                                                <  4:600346/149‑1 (MQ=255)
                                                                    aaTACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCtt                                                                                <  3:204252/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: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                                                                     >  2:466602/1‑110 (MQ=255)
                                                                                                                      tAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGc                                                                     <  1:466602/110‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
AGGAAGTGCAGGTCGCCTTGACGCGCTTGCTCGATAGTATGGGCAATGGCTGGCGGGAGTCCCGCGACAATACCGGCGAAGATAATGATTGAAATACCGTTGCCGATACCTCGTTCAGTAATCTGTTCGCCCAACCACATCAGGAACATGGTTCCTGTGACCAGACTTACAACAGCGGTGAAGTAGAATGCAAAGCCCGGGTTAATCACCAGGCCTTGCATACCAGGCATATTCGGCAGACCGGTAGCAATACCGATCGACTGGAATATTGCCAGCACCAGAGTACCGTAGCGGG  >  NC_000913/3443445‑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: